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Abstract 

We review the new possibilities offered by the reaction dynamics of asymmetric 
heavy ion cohisions, using stable and unstable beams. We show that it represents a 
rather unique tool to probe regions of highly Asymmetric Nuclear Matter {ANM) in 
compressed as well as dilute phases, and to test the in-medium isovector interaction 
for high momentum nucleons. The focus is on a detailed study of the symmetry 
term of the nuclear Equation of State {EOS) in regions far away from saturation 
conditions but always under laboratory controlled conditions. 

Thermodynamic properties of ANM are surveyed starting from nonrelativistic 
and relativistic effective interactions. In the relativistic case the role of the isovector 
scalar 5-meson is stressed. The qualitative new features of the liquid-gas phase tran- 
sition, "diffusive" instability and isospin distillation, are discussed. The results of 
ab-initio simulations of n-rich, n-poor, heavy ion collisions, using stochastic isospin 
dependent transport equations, are analysed as a function of beam energy and cen- 
trality. The isospin dynamics plays an important role in all steps of the reaction, 
from prompt nucleon emissions to the final fragments. The isospin diffusion is also 
of large interest, due to the interplay of asymmetry and density gradients. In rela- 
tivistic collisions, the possibility of a direct study of the covariant structure of the 
effective nucleon interaction is shown. Results are discussed for particle production, 
collective flows and iso-transparency. 

Perspectives of further developments of the field, in theory as well as in experi- 
ment, are presented. 
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1 Introduction 



A key question in the pliysics of unstable nuclei is the knowledge of the EOS 
for asymmetric nuclear matter (ANM) away from normal conditions. We re- 
call that the symmetry energy at low densities has important effects on the 
neutron skin structure, while the knowledge in high densities region is cru- 
cial for supernovae dynamics and neutron star properties. The paradox is 
that while we are planning second and third generation facilities for radioac- 
tive beams our basic knowledge of the symmetry term of the EOS is still 
extremely poor. Effective interactions are obviously tuned to symmetry prop- 
erties around normal conditions and any extrapolation can be quite danger- 
ous. Microscopic approaches based on realistic NN interactions, Brueckner or 
variational schemes, or on effective field theories show a rather large variety 
of predictions. As an example, in Fig. 1-1 we collect the isospin dependence of 
some EOS^s which have the same saturation properties for symmetric NM 
(top): SKM* [2,3], SLy230b (SLyA) [4,5,6] and BPAL32 [7,8,9]. 

In Fig. 1-1 (bottom) we report the density dependence of the potential symme- 
try contribution for the three different effective interactions. While all curves 
obviously cross at normal density po, quite large differences are present for 
values, slopes and curvatures in low density and particularly in high density 
regions. 

Moreover even at the relatively well known "crossing point" at normal den- 
sity the various effective forces are presenting controversial predictions for the 
momentum dependence of the fields acting on the nucleons and consequently 
for the splitting of the neutron/proton effective masses, of large interest for 
nuclear structure and dynamics. 

In the recent years under the stimulating perspectives offered from nuclear as- 
trophysics and from the new Radioactive Ion Beam (RIB) facilities a relevant 
activity has started in the field of the isospin degree of freedom in heavy ion 
reactions, see the refs.[14,15,16]. Here we review the field trying to pin down 
the most interesting theory questions and eventually the related key observ- 
ables. We will follow non-relativistic and relativistic approaches to construct 
effective interactions. In general the physics is not dependent on the theoret- 
ical framework, however we will see some genuine pure relativistic effects in 
the dynamics of the isovector part of the EOS. 

In Section 2, we look at the density dependence of the symmetry term around 
saturation in order to relate slope and curvature to physics properties of ex- 
otic nuclei, bulk densities, neutron distributions and monopole frequencies. 
The momentum dependence of the interactions in the isovector channel is 
thoroughly analysed discussing the expected effects on the energy-slope of the 
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Fig. 1-1. EOS for various effective forces. Top: neutron matter (up), symmetric 
matter (down); Bottom: potential symmetry term. 

Lane Potential at normal density and on the symmetry field seen by high 
momentum nucleons. 

Section 3 is devoted to the study of symmetry properties at low density, i.e. 
to the Liquid-Gas Phase Transition in Asymmetric Matter. We start with 
a thermodynamical approach to the mechanical (vs. density variations) and 
chemical (vs. concentration variations) instability regions. We show then that 
the line of maximum instability for a two component system goes along a mix- 
ing of the two directions with a "mixing angle" that depends on microscopic 
properties of the interaction among the components. 

Wc pass then to a dynamical approach of great relevance since it leads to 
the space-time properties of the unstable normal modes which give rise to 
fragments inside the extended spinodal boundary. Characteristic coupled dis- 
persion relations must be solved. We remind that the coupling of the isoscalar 
and isovector collective response is well known for stable modes [17,18] in 
neutron rich exotic nuclei, with transition densities that show a mixed na- 
ture. Here we extend the same results to the unstable responses of interest in 
fragmentation reactions. The mixed nature of the unstable density oscillations 
win naturally lead to the isospin distillation effect, i.e. a different concentration 
(N/Z) between the liquid and the gas phase [19,20,21,22,23]. This is indeed 
the main new feature of the liquid-gas phase transition which is behind the 
cluster formation in the expansion stage of the reaction dynamics. We can 
expect that quantitatively this effect will be dependent on properties of the 
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symmetry energy in very dilute matter, well below the saturation point. 

The early emission of particles in heavy ion collisions is relevant for under- 
standing both the reaction dynamics and the mechanism for particle produc- 
tion. As discussed before, the density behavior of the symmetry potential is 
strictly related to the value of the pressure that one obtains in asymmetric 
nuclear matter. The early reaction dynamics is mainly governed by the pres- 
sure of the excited nuclear matter produced during the initial stage of the 
collision [24,25,26]. Therefore we can expect important symmetry effects on 
transport properties ( fast particle emission, collective flows) in the asymmet- 
ric NM that will be probed by Radioactive Beam collisions at Fermi and 
intermediate energies. In Section 4 we review several observables that are re- 
lated to the early particle emission: two-nucleon correlation functions, light 
cluster formation, transverse and elliptic collective flows and we will focus on 
the difference observed between neutron and proton behavior. Results of colli- 
sion simulations based on microscopic kinetic equation are discussed, starting 
from realistic effective interactions widely used for symmetric systems. Differ- 
ent parametrizations in the momentum dependence of the isovector channel 
are tested, taking care that the symmetry energy, including its density depen- 
dence, will be not modifled. We study in particular the transport effect of the 
sign of the n/p effective mass splitting m* — m* in asymmetric matter at high 
baryon and isospin density. 

In Section 5, we continue the discussion of nuclear reactions in the Fermi en- 
ergy domain. We remind that this is the transition region between a dynamics 
mainly driven by the mean-field, below 15 — 20 AMeV, and the one where 
the nucleon-nucleon collisions play a central role, above 100 AMeV. We can 
have then a very rich variety of dissipativc reaction mechanisms, with related 
different isospin dynamics. One of the new distinctive features is the enhanced 
production of Intermediate Mass Fragments {IMF, 3 < Z < 20). We analyse 
then with great detail the isovector channel effect on the onset of the frag- 
mentation mechanism and on the isotopic content of the produced fragments. 
In the reaction simulations it is essential the use of a Stochastic Transport 
Approach since fluctuations, instabilities and dynamical branchings are very 
important in this energy range. 

As we may see, even from Fig. 1-1, the potential symmetry term for various 
effective interactions shows quite different behaviors in the region around nor- 
mal density and at very low densities, where we enter the spinodal zone and 
the cluster formation initiates. While around po the density dependence be- 
comes steeper when we go from asy-soft (Skyrme-hke) to asy-superstiff (this 
suggests the names) at subsaturation densities it manifests an opposite trend. 
We will see that with the centrality of the collision we can have different sce- 
narios for fragment production, from the growing instabilities of dilute matter 
in central reactions to the cluster formation at the interface between low and 
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normal density regions in semicentral collisions. We can then expect a large 
variety of isospin effects, probing different regions of the symmetry energy 
below saturation. 

The traditional approach to nuclear physics starts from non-relativistic for- 
malisms in which non-nucleonic degrees of freedom arc integrated out giving 
nucleon-nucleon potentials. Then nuclear matter is described as a collection of 
quantum nonrelativistic nuclcons interacting through an istantaneous effective 
potential. Although this approach has had a great success, a more appropriate 
set of degrees of freedom consists of strongly interacting effective hadron fields, 
mesons and baryons. These variables are the most efficient in a wide range 
of densities and temperatures and they arc the degrees of freedom actually 
observed in experiments, in particular in Heavy Ion Collisions at intermediate 
energies. Moreover this framework appears in any case a fundamental "Door- 
way Step" towards a more microscopic understanding of the nuclear matter. 
Rclativistic contributions to the isospin physics for static properties and reac- 
tion dynamics will be discussed in the second part of the report. Sections 6, 7 
and 8. 

The QHD (Quantum-Hadro-Dynamics) effective field model represents a very 
successful attempt to describe, in a fully consistent rclativistic picture, equi- 
librium and dynamical properties of nuclear systems at the hadronic level 
[27,28,29]. In this report we mostly focus our attention on the dynamical 
response of Asymmetric Nuclear Matter (ANM). We present a rclativistic ki- 
netic theory with the aim of a transparent connection between the collective 
and reaction dynamics and the coupling to various channels of the nucleon- 
nucleon interaction. We show that the same isospin physics described in detail 
at the non-relativistic level can be naturally reproduced. Moreover some new 
genuine rclativistic effects will be revealed, due to the covariant structure of 
the effective interactions. 

One of the main points of our discussion is the relevance of the coupling to 
a scalar isovector channel, the effective 5[ao(980)] meson, not considered in 
the usual nuclear structure studies [12,13]. A related feature of interest is the 
dynamical treatment of the Fock terms. We like to note that recently , see the 
conclusions of the refs. [30,31], the 5-field coupling has been reconsidered as 
an interesting improvement of covariant approaches, in the framework of an 
Effective Field Theory as a rclativistic Density Functional Theory, since 
contribution to this channel are mainly coming from correlation effects. One 
of the main tasks of our work is just to try to select the dynamical observables 
more sensitive to it. 

In Section 6, this extension of the QHD model is presented. Equilibrium 
properties are discussed, like the nuclear EOS and the corresponding ther- 
modynamical instabilities. Particular attention is devoted to the expectations 
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for the splitting of the nucleon effective masses in asymmetric matter, with 
a detailed analysis of the relationship between the Dirac masses of the effec- 
tive field approach and the Schrodinger masses of the non-relativistic models. 
A relativistic Dirac-Lane potential is deduced, which shows a structure very 
similar to the Lane potential of the non-relativistic optical model, but now 
in terms of isovector self-energies and coTipling constants. As a general trend 
we have a close parallelism between relativistic and non-relativistic results, as 
it should be since the main physics is the same. We always note this point 
but we also stress several new features coming from a consistent use of a field 
theory approach. 

We follow a Relativistic Mean Field (RMF) approximation that is allow- 
ing more physics transparent results, often even analytical. We will always 
keep a close connection to the more microscopic Dirac-Brueckner-Hartree-Fock 
(DBHF) approaches, in their extension to asymmetric matter, [32,33,34,35] It 
is well known that correlations are naturally leading to a density dependence 
of the coupling constants, see ref. [33] for the DBHF calculations and ref. [36] 
just for the basic Fock correlations. Within the RMF model we can get a clear 
qualitative estimation of the contribution of the various fields to the nuclear 
dynamics. The price to pay is that when we try to get quantitative effects 
we are forced to use different sets of couplings in different baryon density re- 
gions. In this case we use the DBHF results as guidelines. In particular for 
the controversial 5-mcson field, expected to be important at densities above 
saturation, we fix the corresponding coupling from the analysis of refs. [34,35], 
where it actually appears not strongly density dependent in a wide range of 
baryon densities. 

An important outcome of our work is to show that the two effective cou- 
plings, vector and scalar, in the isovector channel are influencing in a different 
way the static (symmetry energy) and dynamic (collective response, reaction 
obscrvables) properties of asymmetric nuclear matter. All that will open new 
possibilities for a phenomenological determination of these fundamental quan- 
tities. 

In Section 7, we discuss a fully relativistic Landau Fermi liquid theory based on 
the Quantum Hadro- Dynamics (QHD) effective field picture of Nuclear Mat- 
ter. From the linearized kinetic equations we get the dispersion relations of 
the propagating collective modes. The relation between static properties and 
the collective response is analysed stressing the different role of the various 
channels present in the effective nuclear interaction. We focus our attention 
on the dynamical effects of the interplay between scalar and vector field con- 
tributions. An interesting "mirror" structure in the form of the dynamical 
response in the isoscalar /isovector degree of freedom is revealed, with a com- 
plete parallelism in the role respectively played by the compressibility and the 
symmetry energy. In particular we study the infiuence of a scalar-isovector 
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channel (coupling to a (5-mcson-like effective field) on the collective response 
of asymmetric nuclear matter (ANM). Interesting contributions are foTind on 
the propagation of isovector-like modes at normal density and on an expected 
smooth transition to isoscalar-hke oscillations at high baryon density. 

Important "chemical" effects on the neutron-proton structure of the normal 
modes are shown. For dilute ANM we have the isospin distillation mechanism 
of the unstable isoscalar-like oscillations, as already shown in Section 3 in a 
non-relativistic frame, while at high baryon density we predict an almost pure 
neutron wave structure of the propagating sounds. 

Results for relativistic Heavy Ion Collisions (HIC) in the AGeV beam energy 
region are presented in Section 8. We recall that intermediate energy HICs 
represent the only way to probe in terrestrial laboratories the in-medium ef- 
fective interactions far from saturation, at high densities as well as at high 
momenta. Within a relativistic transport model it is shown that the isovector- 
scalar 5-meson, which affects the high density behavior of the symmetry term 
and the nucleon effective mass splitting, influences the isospin dynamics. The 
effect is largely enhanced by a relativistic mechanism related to the covariant 
nature of the fields contributing to the isovector channel. The possibility is 
emerging of a direct measurement of the Lorentz structure of the effective 
nuclear interaction in the isovector channel. 

Quantitative calculations are discussed for collective flows, charged pion pro- 
duction and isospin stopping. Asymmetric systems where some data are avail- 
able have been studied. Although the data are mostly of inclusive type (and 
the colliding nuclei not very neutron rich), quite clearly a dependence of some 
observables on charge asymmetry is emerging. Very sensitive quantities ap- 
pear to be the elliptic flow, related to the time scale of the particle emissions, 
and the isospin transparency in central collisions. 

Finally in Section 9 a general outlook for theory and experiment is presented. 
Particular attention is paid to a selection of the expected most sensitive observ- 
ables to the isospin dynamics, in heavy ion reactions from low to relativistic 
energies. 
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2 Symmetry term effects on compressibility, saturation density 
and nucleon mean field 



In asymmetric matter the energy per nucleon, i.e. the equation of state, will 
be a functional of the total (p = p„ + Pp) and isospin (p3 = p„ — pp) densities. 
In the usual parabolic form in terms of the asymmetry parameter / = ps/p = 
{N — Z) / A we can define a symmetry energy -^^(p): 

= + (2-1) 



The symmetry term gets a kinetic contribution directly from the basic Pauli 
correlations and a potential contribution from the properties of the isovector 
part of the effective nuclear interactions in the medium. Since the kinetic 
part can be exactly evaluated we can separate the two contributions, reducing 
the discussion just to a function F{u) of the density u = p/po linked to the 
interaction: 

e.,.^^(p) = ^ + f F(.), (2-2) 



with F{1) — 1, where po is the saturation density and the parameter C is of 
the order C ~ 32MeV to reproduce the 04 term of the Bethe-Weiszacker mass 
formula. The validity of the parabolic form Eq.(2-1) comes directly form the 
isospin structure of the nucleon- nucleon interaction. The effects of many-body 
correlations are not much affecting this behaviour, see ref. [32,37], well verified 
in the nuclear data systematics [38]. Corrections are neglible up to very high 
densities, of the order of 5po- The corrections in the kinetic term are also 
exactly evaluated to be of 4% for any density. 

A traditional expansion to second order around normal density is used [39,40,12] 
esym = -^{P) - a4 + - [-^) + [-^) , (2-3) 

in terms of a slope parameter 

i . 3p.(%=) ^ 1P..M (2-4) 

^ dp ) P=PQ Pn 



simply related to the symmetry pressure Psym = P^desym/dp at and a 
curvature parameter 
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Fig. 2-1. Geometric picture of the lowering of the saturation density in asymmetric 
matter, Eq.(2-6). 

a kind of symmetry compressibility. We remark that our present knowledge 
of these basic properties of the symmetry term around saturation is still very 
poor, see the recent analysis in ref.[30] and refs. therein. In particular we 
note the uncertainty on the symmetry pressure at po, of large importance for 
structure calculations. These points, in connection to the possibility of getting 
a new insight from reaction data, will be largely discussed in this report. 

We first mention some simple considerations on asymmetry effects on equi- 
librium density and compressibility, observables related respectively to bulk 
densities and monopole resonances in medium heavy nuclei [10] . From a linear 
expansion around the value at symmetry, 7 = 0, we get for the variation of 
saturation density 



_9p^ ± 

Knm{I = 0) dp^ 
3poL 



P=Po 



Knm{I = 0) 



P < 0, (2-6) 



where Kj^m{I — 0) and po ^-re respectively compressibility and saturation 
density of symmetric NM. 

The Eq.(2-6) has an intuitive geometrical meaning, which is qualitatively 
shown in Fig.2-1. Asymmetry brings an extra pressure Pgym — P^d^sym/dp 
that can be compensated just moving to the left the saturation point (P = 0) 
of the quantity Pgym/ {dP/ dp) at po ( we recall that, for symmetric matter 
at Po, 9{dP/dp) = K^m{I = 0)). For the compressibility shift we have, after 
some algebra. 



AKnm{I)^9po 



'dp" 



d n 

dp- 



'Sym 



iP) 



P=PO 



Table 1 

Symmetry term at saturation 



F{u) 


L 




[Ksyni ~ 6ij] 


const = 1 


+25MeV 


-25MeV 


-125MeV 


\/u 


+4QMeV 


-61MeV 


-355MeV 


u 


+75MeV 


-25MeV 


-4:75MeV 


u''/{l+u) 


+100MeF 


+50MeV 


-550Mey 



= [K^ym - GL]P < 0, (2-7) 

where we note the interplay between slope and curvature of the symmetry 
term. 

The Eq.(2-6), and the related geometrical interpretation Fig. 2-1, is a particular 
case (for P = 0) of the general variation of the density in asymmetric matter 
corresponding to a fixed pressure P, see the discussion in rcfs. [41,42]. The 
corresponding compressibility shift, Eq.(2-7), is also evaluated in the isobaric 
case, where there is the possibility to have some experimental information 
from Giant Monopole Resonance (GMR) data in charge asymmetric nuclei. 
Actually Eq.(2-7) represents an approximate form where higher order terms 
in density variations are neglected, see refs. [41,42]. The general negative sign 
of the compressibility shift, Eq.(2-7), is a natural consequence of Eq.(2-6), i.e. 
of the fact that the saturation density is decreasing with charge asymmetry. 

In order to have a quantitative idea, we now show explicitly the influence 
on the L, Ksym parameters of a different density dependence in the potential 
part of the symmetry energy around saturation, i.e. of the function F{u) of 
Eq.(2-2): 

L^-eF + -C^Fiu) 

3 2 du ^ ' y-=-i- 

We obtain the very instructive Table 1 for various functional forms F{u),u = 
yo/po, around po- 



The choices in the table of the F{u) behaviors are in fact not arbitrary. They 
reflect the wide spectrum of theory predictions for effective forces in the isovec- 
tor channel, as discussed in detail in the rest of the report. The constant trend 
of the flrst row (around saturation density) is typical of Skyrme-like forces, 
the -/u behavior is obtained in variational approaches with realistic A^A^ inter- 
actions [43,44], the linear dependence in Brueckner-Hartree-Fock (BHF) non 
relativistic calculations [7,8,37] as well as in Relativistic Mean Field (RMF), 
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[10], and Dirac-Brueckner-Hartree-Fock (DBHF), [32], approaches. Finally 
the more repulsive, nearly parabolic, dependence of the fourth row, [45,46], 
can be related to non-relativistic predictions with three-body forces, either 
extended BHF or variational [37,47] as well as to other relativistic DBHF 
estimations, refs. [33,34], or RMF with scalar isovector meson-hke contribu- 
tions [12,13]. 

Fixing the same compressibility at saturation for symmetric matter Kmm{I — 
0), the relativistic effective forces always predict larger shifts in both equihb- 
rium density and compressibility for asymmetric matter. This can be seen also 
in the comparison with Skyrme-like interactions in the finite nuclei calcula- 
tions shown in rcf. [10], even at relatively small charge asymmetries. There are 
therefore good chances of obtaining some direct experimental indications from 
GMR measurements in N ^ Z nuclei. A recent systematic study of the isospin 
dependence of GM 7?'s in Sn isotopes, [48] , seems to reveal a relatively large 
decrease of the Giant Monopole centroid with increasing asymmetry, more in 
agreement with a stiff behavior of the symmetry term around saturation, 
as expected in the relativistic models. We note that the relation between the 
"symmetry compressibility", Eq.(2-5), and the variation of the compressibility 
of asymmetric matter is not trivial: we can have cases where an increase of 
Kgyni actually corresponds to a softening of the Equation of State of asym- 
metric matter. 

It is finally instructive to evaluate the density gradient of the symmetry pres- 
sure as a function of the slope and curvature of the symmetry term: 

d 2 1 



that around normal density gives 
d ^ 10 



-^Psyra = — + C — + --p^ F {u) 



dp 27 



d 1 d^ 



du 2 du'^ 



u=l 



A stiffer symmetry term in general enhances the pressure gradient of asym- 
metric matter. We can expect direct effects on the nucleon emissions in the 
reaction dynamics, fast particles and collective flows. Moreover due to the 
different fields seen by neutrons and protons, we shall observe even specific 
isotopic effects. This point will be analysed in detail in Sect. 4. 



2. 1 The symmetry term of Skyrme forces 



Since in this report we will often show reaction results from non-relativistic 
kinetic equations with Skyrme forces [2] , we will expand a little the discussion 
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Fig. 2-2. Density dependence of the potential symmetry term for various Skyrme 
effective forces, see text. The bottom right panel shows the kinetic contribution 

on the isospin dependence of these widely used effective interactions. In a 
Skyrme-like parametrization the symmetry term has the form: 



^ Esym . . _ epip) C{p) p_ 

esym- ^ 3 + 2 po 



(2-9) 



with the function C(p), in the potential part, given by: 



C{P) 
Po 

+ 
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to(l + 2X0) + 1(1 + 2x3) P" 



(2-10) 



with a > and the usual Skyrme parameters. Wc remark that the second 
term is related to isospin effects on the momentum dependence [4] . In the Fig. 
2-2 we show the density dependence of the potential symmetry term of various 
Skyrme interactions, 5"///, SGII, SKM* (see [3] and refs. therein) and the 
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Tabic 2 



Parameters of Skyrme Forces 



Force 


SIII 


SGII 


SkM* 


SLya 


SLyA 


SLyl 


to{MeVfm^) 


-1128.75 


-2645.0 


-2645.0 


-2490.3 


-2488.91 


-2482.41 


i3(MeF/m3+3°) 


14000.0 


15595.0 


15595.0 


13803.0 


13777.0 


13677.0 


xq 


0.45 


0.09 


0.09 


1.1318 


0.8340 


0.846 




1.0 


0.06044 


0.0 


1.9219 


1.3539 


1.391 


ti{MeVfm^) 


395.0 


340.0 


410.0 


489.53 


486.82 


457.97 


t2{MeVfm^) 


-95.0 


-41.9 


-135.0 


-566.58 


-546.39 


-419.85 


Xl 


0.0 


-0.0588 


0.0 


-0.8426 


-0.3438 


-0.511 




0.0 


1.425 


0.0 


-1.0 


-1.0 


-1.0 


a 


1 


1/6 


1/6 


1/6 


1/6 


1/6 



more recent Skyrme-Lyon forms, SLya and Slyh (or SlyA), see [4,5]. We also 
separately present the local and non-local contributions, first and second term 
of the Eq.(2-10). We clearly see a sharp change from the earlier Skyrme forces 
to the Lyon parametrizations, with almost an inversion of the signs of the 
two contributions. The important repulsive non-local part of the Lyon forces 
leads to a completely different behavior of the neutron matter EOS, of great 
relevance for the neutron star properties. Actually this substantially modified 
parametrization was mainly motivated by a very unpleasant feature in the 
spin channel of the earlier Skyrme forces, the collapse of polarized neutron 
matter, see discussion in [49,4,5,6]. 

In the Table 2 we collect the parameters, of interest for Nuclear Matter prop- 
erties, of the various Skyrme forces used here, where we have also included 
the recent Sly7 force particularly tuned for neutron rich systems [6] . Note the 
transition to X2 = — 1 values in the Skyrme-Lyon Forces just to cure the col- 
lapse of ferromagnetic neutron stars discussed before. In correspondence the 
predictions on the isospin effects on the momentum dependence of the sym- 
metry term arc quite different, see Fig. 2-2. A very important consequence for 
the reaction dynamics is the expected inversion of the sign of the n/p effective 
mass splitting, which will be widely discussed in the next sections. 



2.1.1 Mean field and chemical potentials 

We can derive a general Skyrme-like form for neutron and proton mean field 
potentials [50,51]: 



12 




9€pot{Pq, pq') 
dpg 




o+l 



(2-11) 



Spot being the potential energy density. Here p = p^ + Pp and p^ = pn — Pp are 
respectively isoscalar and isovector densities, and q — n,p, Tg — +1 {q — n), 



In the following we will always compare results obtained with forces that 
have the same saturation properties for symmetric NM [52]. We will refer 
to an "asy — stiff EOS (e.g. hke BPAL32 of Fig.l) when we are con- 
sidering a potential symmetry term linearly increasing with nuclear density 
and to a "asy — so/t" EOS (e.g. like SKM* of Fig.l) when the symmetry 
term shows a saturation and eventually a decrease above normal density. In 
some cases, in order to enhance the dynamical effects, we will consider also 
"asy — super sti f f behaviours, i.e. with a roughly parabolic increase of the 
symmetry term above normal density [9,45,46]. Wc focus our discussion on 
single particle properties since in this case the symmetry contribution will be 
linearly dependent on the asymmetry of the matter. 

In Figs.2-3, 2-4 we report, for an asymmetry (N — Z)/A = 0.2 representative of 
^"^"^Sn, the density dependence of the symmetry contribution to the mean-field 
potential (Fig.2-3, second part of Eq.(2-ll)) and of the chemical potentials 
(Fig. 2-4) for neutrons (top curves) and protons (bottom curves) , for the 
different effective interactions in the isovector channel. 

From Fig.2-3 we note that in regions just off normal density the field "seen" 
by neutrons and protons in the three cases is very different, in particular 
below saturation density. We thus expect important transport effects during 
reactions at intermediate energies (prompt particle emissions, collective flows, 
n/p interfcrometry) : the interacting asymmetric nuclear matter will experience 
compressed and expanding phases before forming fragments around normal 
density. 

A transparent picture of the isospin dynamics can be obtained from the anal- 
ysis of the density dependence of the "local" values of the neutron/proton 
chemical potentials Pq = de{pq, pqr)/dpq, e being the energy density. We re- 
mind that the chemical potentials contain all the contributions to the energy 
per particle, including the isoscalar part and the kinetic symmetry term. In 
non-equilibrium processes the mass flow is determined by the differences in 
the local values of chemical potential and it is directed from the regions of 
higher chemical potential to regions of lower values until equalization. 

From Fig. 2-4 we can already predict the Isospin Distillation effect and even 



-1 (q^p). 
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Fig. 2-3. Symmetry contribution 
to the mean field at / = 0.2 

for neutrons (upper curves) 
and protons (lower curves): 
dashed lines "asy-soft", long 
dashed lines "asy-stifF', solid lines 



Fig. 2-4. Density dependence, for 
/ = 0.2, of neutron (upper curves) 
and proton (lower curves) chemical 
potentials for asy-supcrstiff (solid 
lines) and asy-soft (dashed line) 



EOS. 



" asy-superstifF' 

the differences between the results in the isospin dynamics during fragment 
formation obtained using different symmetry terms, see Sects. 3 and 5. The 
chemical potential for protons in a system having the ^^'^Sn asymmetry, / = 
0.2, is below the corresponding value of symmetric nuclear matter while for 
neutrons is above, indeed /U„ — /Xp = Aesym{p)I. From the density dependence 
in the low densities region we see that when the inhomogenities develop both 
neutrons and protons have the tendency to move in phase from lower to higher 
density regions. This is in qualitative agreement with the rigorous proof of 
the fact that the system is unstable against isoscalar-like fluctuations, see 
next section. Since the variations of the two chemical potentials are different 
(larger for protons) we expect a lower asymmetry in the liquid phase, i.e. in 
the clusters formed through such bulk instability mechanism. Moreover from 
the larger difference in the neutron/proton slopes in dilute matter, around 
l/3poi for the asysoft case (dashed lines of Fig.2-4) we can even expect a 
larger Isospin Distillation mechanism with such symmetry term. 

In the case of a clusterization in presence of a contact between more dilute 
and "normal" density regions in order to understand the isospin dynamics we 
have to look at the density dependence of proton/neutron chemical potentials 
in the region between 0.08 fm^^ and 0.16 fm^^. We see from Fig.2-4 that in 
this range the neutrons have the tendency to move towards more dilute regions 
producing a n-enrichment while the protons will migrate in opposite direction. 
Such mechanism is present in the "neck fragmentation", [53,54,55,56]: the neck 
IMF' s will be always more n— rich compared to the fragments produced in 
the case of bulk fragmentation. This effect, clearly seen in experiments, will 
be discussed in detail in Sect. 5 since naturally it appears very sensitive to the 
stiffness of the symmetry term around saturation density. 
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2.2 Effective masses in neutron-rich matter 



In the Fig. 2-2 we have noticed a dramatic change in the local vs. non-local 
contributions to the symmetry energy going from the old to the new Skyrme 
forces of the Lyon type, explicitly built for asymmetric matter. A related 
interesting effect can be seen on the neutron/proton effective masses. 

In asymmetric matter we consistently have a splitting of the neutron/proton 
effective masses given by: 

m*~^ = + gip + g2pq, (2-12) 

with 

Pq=n,p 2 P ■ 

The gi, g2 coefficients are simply related to the momentum dependent part 
of the Skyrme forces: 



gi^^[tl{2 + Xi)+t2{2 + X2)] 

92 = + 2^2) - ti(l + 2xi)] (2-13) 

This result derives from a general q — structure of the momentum dependent 
part of the Skyrme mean field 

Uq,MD = m{9iP + 92pq)E (2-14) 



where E is the nucleon kinetic energy. 

In the Fig. 2-5 we show the density behavior of in neutron rich matter 
/ = 0.2 for the same effective interactions. From the Eqs.(2-12, 2-13) we see 
that the sign of the g2 univocally assigns the sign of the splitting, i.e. g2 < 
gives larger neutron masses m* > m* while we have the opposite for g2 > 0. 

In the Table 3 we report some results obtained with various Skyrme forces 
for quantities of interest, around saturation, for the present discussion. We 
show also the E — slope of the corresponding Lane Potential, see later, simply 
related to the isospin dependent part of Eq.(2-14). For the effective mass 
parameters of Eq.(2-12) we observe that while the gi coefficients are always 
positive, corresponding to a decrease of the nucleon mass in the medium, the 
isospin dependent part shows different signs. In particular we see that in the 
Lyon forces the g2 values are positive, with neutron effective masses below the 
proton ones for n-rich matter as shown in Fig. 2-5. 
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Fig. 2-5. Density dependence of the neutron/proton effective mass splitting for var- 
ious Skyrme effective forces, see text. The asymmetry is fixed at / = 0.2, not very 
exotic. 

Table 3 

Properties at saturation 



F 01 ■('■(' 


SIII 


SGII 


SkM* 


SLjja 


SLjji 


SLyJ 


gi {10-^){MeV-^fm^) 


+3.85 


+3.31 


+3.53 


+10-5 


+1.67 


+1.70 


g2 (10-^) {MeV-^fm^) 


-3.14 


-2.96 


-3.50 


+5.78 


+2.53 


+2.76 




0.150 


0.1595 


0.1603 


0.160 


0.1595 


0.1581 


a^iMeV) 


28.16 


26.83 


30.03 


31.97 


32.01 


32.01 


C{po){MeV) 


31.72 


29.06 


35.46 


39.40 


39.42 


39.42 


E — slope{LanePot.) 


-0.22 


-0.21 


-0.26 


+0.43 


+0.19 


+0.20 



In general we obtain a splitting of the order of 10 — 15% at normal density po, 
and increasing with baryon density. Unfortunately from the present nuclear 
data we have a very little knowledge of this effect, due to the low asymmetries 
available. This issue will be quite relevant in the study of drip-line nuclei. 

The sign itself of the splitting is very instructive. Passing from Skyrme to 
Skyrme- Lyon we see a dramatic inversion in the sign of the n/p mass split- 
ting. The Lyon forces predict in n-rich systems a neutron effective mass always 
smaller than the proton one. We will come back to this point. Here we just 
note that the same is predicted from microscopic relativistic Dirac-Brueckner 
calculations [34] and in general from the introduction of scalar isovector virtual 



16 



mesons in RMF approaches [12,13]. At variance, non-relativistic Brueckner- 
Hartree-Fock calculations arc leading to opposite conclusions [57,58]. We re- 
mind that a comparison between relativistic effective (Dirac) masses and non- 
relativistic effective masses requires some attention. This point will be care- 
fully discussed later in the effective field theory approach to the in-medium 
interactions. Sect. 6. 



2.2.1 Effective masses and Landau Parameters 

The quasiparticle nucleon energies in a general Skyrme form are given by: 



(27r/i) 



where Uq{pq, Pg') is the local part of the mean field, Eq.(2-ll), and f{p),fqip) 
are the nucleon momentum distributions. The corresponding Landau param- 
eters defined by: 

^''^tJ (^[^^^(^^'^O^/'^b') + FgAP:P'WAp')] (2-16) 



have the form: 



Fqqip,p') = 7V^(|(p - + |(p - p'r + 

FqAp,p')-N,{^ip-pr + §^) (2-17) 



where Nq = is the energy level density of the q-nucleons. 
Prom the expansion: 

F<i,' = F^q' + Fiq,{p-p') (2-18) 

we get, at the Fermi momentum, the explicit Skyrme form of the "local" and 
"non-local" Landau parameters: 



F^q = Nq(p^) [(g, + g^)p\ + ^ 



dpqdpq. 
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Fgq' = -NgipF)giPFqPFq' 



(2-19) 



The (n,p) effective masses have the compact form: 



m* 1 
= 1 + - 
m 3 



(2-20) 



which nicely leads to the (1 -|- \F^) result for symmetric matter [59]. For the 
(n,p) mass splitting in asymmetric matter we have the expression: 



m 



F' + i?^^' 



'PFn' 



J np \ 



PFn 
PFp 



(2-21) 



Since all the F^^, are negative this result has been used to predict a larger m* 
mass in n-rich systems due to the larger neutron Fermi momentum, see the 
recent ref.[60]. However this is actually not generally correct because in fact 
all the F^ Landau parameters also depend on Fermi momenta and the final 
balance will be fixed by the microscopic structure of the effective interaction, 
i.e. by the interplay between the gi and g2 quantities in the case of Skyrme- 
like forces, Eq.(2-19). In particular, we see that when we pass from the "old" 
Skyrme to the "Lyon" parametrizations, i.e. from g2 < to g^ > 0, the 
combination gi + g2 becomes much larger and the term is dominant in 
the n-rich case, see Eq.(2-21), leading to a m* < m* mass splitting, as already 
discussed. 



2.2.2 Energy dependence of the Lane Potential 

We note that the sign of the splitting will directly affect the energy dependence 
of the Lane Potential, i.e. the difference between {n,p) optical potentials on 
charge asymmetric targets, normalized by the target asymmetry [61]. From 
the Eqs.(2-ll, 2-14) we obtain the explicit Skyrme form of the Lane Potential: 

U,ane ^ = C {po) + ^g^E (2-22) 



where C(po) gives the potential part of the 04 parameter in the mass formula, 
see Eq.(2-2). We see that the E — slope has just the sign of the g2 parameter, 
and so we have opposite predictions from the various Skyrme forces analysed 
here. The change in the energy slope is reported in the last row of the Table 



18 



150 



100 



> 

(D 



CN 



1 1 1 1 1 1 1 1 1 

: — m*>m*p 


1 1 1 1 1 1 1 1 1 


_ m,<mp 


SLy7 _ 




_,._^.:::^---'-''''"SLy4: 


' '''X. ■ — ■ — 


- 


1 1 1 1 1 1 1 1 1 


s\\\'--: 
1 1 1 1 1 1 1 1 1 



100 200 300 

E.n (MeV) 



400 



Fig. 2-6. Energy dependence of the Lane Potential for different Skyrme forces. The 

dotted (m* < m*) and solid (m* > m*) curves correspond to the more general 
momentum dependent mean fields BGBD — 1,2 of the ref.[65], see text. 



3. The difference in the energy dependence of the Lane Potential is quite 
dramatic, as we can see from Fig.2-6. 

An important physical consequence of the negative slopes is that the isospin 
effects on the optical potentials tend to disappear at energies just above 
100 MeV (or even change the sign for "old" Skyrme- like forces). Unfortu- 
nately results derived from neutron/proton optical potentials at low energies 
are not conclusive, [61,62,63], since the effects appear of the same order of 
the uncertainty on the determination of the local contribution. More neutron 
data are needed at higher energies, in particular a systematics of the energy 
dependence. 

Moreover we can expect important effects on transport properties ( fast par- 
ticle emission, collective flows) of the dense and asymmetric NM that will 
be reached in Radioactive Beam collisions at intermediate energies. Indeed at 
supra-saturation density the difference in the predictions will be enhanced. 
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2.3 Isospin effects on the momentum dependence of the mean field 



In presenting the symmetry energy results of various effective Skyrme forces 
we have stressed the conflicting predictions on the isospin dependence of the 

effective masses, i.e. on Isospin Momentum Dependent, Iso-MD, effects. This 
will be one of the main questions to address in the reaction dynamics of 
exotic nuclear systems, in particular for the close connection to fundamental 
properties of the nuclear interaction in the medium. 

We will review some results on this direction using non-relativistic and rela- 
tivistic (later in Sects. 6 and 8) microscopic kinetic approaches. Starting from 
reahstic effective interactions widely used for symmetric systems we will test 
very different parametrizations in the momentum dependence of the isovector 
channel, taking care that the symmetry energy, including its density depen- 
dence, will be not modified. We study in particular the transport effect of the 
sign of the n/p effective mass splitting m* —m* in asymmetric matter at high 
baryon and isospin density. 

In a non-relativistic frame we can extend the general form of effective mo- 
mentum dependent interactions first introduced by Bombaci et al. [64,57] for 
astrophysical and heavy ion physics applications. The isovector part can be 
modified in order to get new parametrizations with different n/p effective mass 
splittings while keeping the same symmetry energy at saturation, including a 
very similar overall density dependence, see ref.[65]. 

The energy density as a function of the asymmetry parameter reads: / = 



^{pni Pp) — ^kin + + + ^MD 



(2-23) 



(2n) 



2 



. ^. B 2 B p'^+i /I \ ^ 



SMoiPn: Pp: I) = C—{In + 1p) + ^ J'^^ —I{In - Ip) 
Po 5 

where the integrals 2r(-^) = -j^^ I d?kfT{k)g{k, A) include the momentum 



k-<k> 
A 



; the subscript 



dependent part of the mean field g{k, A) 
T — n,p stands for neutrons and protons respectively; — 0.16 /m~^ is the 



1 -1 
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normal density of nuclear matter. We refer to this parametrization as the 
BGBD — EOS. For symmetric nuclear matter (/ = 0), the energy density 
Eq.(2-23) reduces to the parametrization proposed by Gale, Bertsch and Das 
Gupta [GBD interaction, [66]). In a sense this interaction represents an exten- 
sion of the Gogny non locality since it gives very similar results in the range 
p/po ^ 1-5 and k < 4fm~^, see refs. [67,68] where nice applications in the 
Fermi energy region can be found. 

The parameters A, B, C, a and A take the same values as in [66] {A = 
-lUMeV, B = 203.3 MeV, C = -75MeV, a = |, A = 1.5pP, where pP is 
the Fermi momentum at normal density), and provide a soft EOS for sym- 
metric matter, with a compressibility -R'ata/ ~ 210 MeV. The value of Zi sets 
the strength of the momentum dependence [MD) in the isospin channel; the 
remaining parameters xq and can be set to fix the symmetry energy. 

Taking the functional derivative of the energy density with respect to the 
distribution function f^-, we obtain, apart from the kinetic term Skin, the mean 
field potential for neutrons and protons: 




p/po (a. J.) p/pc (a, u.) 



Fig. 2-7. Potential symmetry energy as a function of density (solid line) for m* < ni* 
{BGBD - 1 set, left) and m* > m* (BGBD - 2 set, right). Dotted lines refer 
to density-dependent contributions, dashed-dotted to momentum dependent ones. 
Small panels indicate the corresponding behaviour of proton and neutron effective 
masses as a function of density, with I = 0.2. 
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^2-24) 



5Po 



|^(3C-4zi)X, + (C + 2zi)T,,| 



where the subscripts in the integrals are t ^ t'\ the upper signs refer to 
neutrons, the lower ones to protons. 

The last BGBD term includes, besides the usual GBD momentum depen- 
dence, an isospin-dcpcndcnt part from which we can get different effective 
masses for protons and neutrons. In fact, the effective mass is defined as: 



and we see that mass splitting is determined not only by different momentum 
dependence of mean field, but also by Fermi momenta of neutrons and protons. 




Fig. 2-8. Mean field potential as a function of momentum at normal density, for an 
asymmetry / = 0.2. Solid lines refer to the case m* > m* {BGBD — 1 set), dashed 
lines to m* < m* {BGBD - 2 set). 
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In the case considered here we have: 



ml 



m 



— < 



1 + T 



!+(¥)' 



(2/3) / p x(2/3) 



1 2 



In order to investigate the effects of mass splitting on non-relativistic heavy 
ion colhsions two sets of parameters {BGBD — 1,2 shown in Table 4) which 
give opposite splitting, but quite similar behaviour of the symmetry energy, 
will be considered: 



force 


mass splitting 




Xo 


3^3 


BGBD-1 


m* < m*p 


28 


1.925 


0.41 


BGBD-2 


ml > m* 


-36.75 


-1.477 


-1.01 



Table 4 

Values of the parameters zi, xq and X3, for opposite mass splitting but giving the 
same Esym{po) = ^3AMeV, which characterize the two used effective forces. 



Figure 2-7 shows the potential part of the symmetry energy as a function 
of density (solid line) for the two choices of mass splitting. We remark the 
very similar overall density dependence. The dashed-dotted and dotted lines 
indicate respectively the contributions from the momentum dependent and 
density-dependent part of the EOS. A change in the relative sign of mass 
splitting is related to opposite behaviours of these two contributions, exactly 
as already noticed in the Introduction for the Skyrme-like forces. The small 
panels on the top left of each graph illustrate the corresponding mass splitting 
as a function of density for an asymmetry I — 0.2 (the ^^"^Au asymmetry). 

We finally discuss the relation between effective mass and momentum depen- 
dence. Prom the definition Eq.(2-25) we see that the effective mass is inversely 
proportional to the slope of mean field at the Fermi momentum: 

• For n-rich systems the n/p mean field potential difference t/„ — Up will 
increase with nucleon momentum in the m* < m* case and decrease in the 
opposite m* > m* case. 

• For nucleons with smaller effective masses the potential will be more re- 
pulsive at momenta higher than the Fermi one, and more attractive at low 
momenta. 

Mean field potentials at normal density as a function of momentum are shown 
in Fig.2-8 for asymmetry / = 0.2. As we can see, the parametrizations de- 
scribed here give rise to opposite behaviors for low and high momentum parti- 
cles. This is a very general feature, due to the meaning itself of effective mass 
Eq.(2-25). As we will see in the reaction dynamics this behavior can give rise 
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to some compensation effects between low and high momenta contributions. 
This will lead to an expected larger isospin — M D sensitivity of more exclu- 
sive measurements, in particular with a tranverse momentum selection of the 
nucleons emitted for a given rapidity. 
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3 Instabilities in Two-component Fluids: the Liquid-Gas Phase 
Transition in Asymmetric Matter 



Since the fragment production represents a relevant dissipation mechanism 
for intermediate energy reactions and since it is seen as a consequence of the 
hquid-gas phase transition in Asymmetric Nuclear Matter (ANM), we are 
devoting one section to a detailed discussion of it. In this way we set some 
guidelines in order to understand the physics behind the transport simulation 
results of Sect. 5, in particular on the possibility of extracting some informa- 
tion on the nuclear effective forces in the isovector channel. A very transpar- 
ent picture of the Isospin Distillation, actually Neutron Distillation, effect is 
emerging. 

One-component systems may become unstable against density fluctuations 
as the result of the strong attraction between constituents. In symmetric bi- 
nary systems, like Symmetric Nuclear Matter (SNM), one may encounter 
two kinds of density fluctuations: i) isoscalar, when the densities of the two 
components oscillate in phase with equal amplitude, ii) isovector when the 
two densities fluctuate still with equal amplitude but out of phase. Then me- 
chanical instability is associated with instability against isoscalar fluctuations 
leading to cluster formation while chemical instability is related to instability 
against isovector fluctuations, of repulsive character, leading to species sepa- 
ration. Turning now to asymmetric binary systems, as the ANM of interest 
here, this direct correspondence between the nature of fluctuations and the 
occurrence of mechanical or chemical instabilities is lost and wc face a more 
complicated scenario, where isoscalar and isovector instabilities are coupled 



An appropriate framework for the study of instabilities is provided by the 
Fermi liquid theory [69], which has been applied, for instance, to symmetric 
binary systems as SNM (the two components being protons and neutrons) 
[70] and the liquid ^He (spin-up and spin-down components) [71,72]. 

3. 1 Thermodynamical study 

Let us flrst discuss the thermodynamical stability of ANM at T = 0. We 
will review here the results obtained in [20] with some more discussion on the 
physics observables. The introduction of flnite temperatures is straightforward. 

The starting point is an extension to the asymmetric case of the formalism 
introduced in [71]. The distribution functions for protons and neutrons are: 



[19,20,22]. 




(3-1) 
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where /Xg arc the corresponding chemical potentials. The nucleon interaction 
is characterized by the Landau parameters: 

^ N y^J!lL_ . N N(T)=[ (3-2) 

where H is the energy density, V is the volume and iV^ is the single-particle 
level density at the Fermi energy. At T = this reduces to 

N,{0) = mpFjiir^K') = 3p,/(2e^,,), 

where pF,q and eF,q are Fermi momentum and Fermi energy of the g-component. 
Thermodynamical stability for T = requires the energy of the system to be 
an absolute minimum for the undistorted distribution functions, so that the 
relation: 

6H - Hp6pp - Hn^pn > (3-3) 



is satisfied when we deform proton and neutron Fermi seas. 

Only monopolar deformations will be taken into account, since we consider 
here momentum independent interactions, so that Fj^l^'' are the only non- 
zero Landau parameters. In fact, for momentum independent interactions, 
all the information on all possible instabihties of the system is obtained just 
considering density variations. However one should keep in mind that in the 
actual dynamical evolution of an unstable system in general one observes 
deformations of the Fermi sphere, hence the direction taken by the system in 
the dynamical evolution is not necessarily the most unstable one defined by 
the thermodynamical analysis. 

Then, up to second order in the variations, the condition Eq.(3-3) becomes 
SH - PpSpp - PnSpn = ^{0,6 pp^ + h6pn + c6pp6pn) > (3-4) 

where 

a^N-\0){l + Fn ; b^N-\0){l + Fn ; 

c = N;\0)Fr + N-\0)F^' = 2N;\0)Fr. (3-5) 
The r.h.s. of Eq. (3-4) is diagonalized by the following transformation: 



u — cos 13 6pp + sin(3 6pm 
V — —sin[3 Spp + cosP Spn, 



(3-6) 
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where the mixing angle 0</3<7r/2is given by 

^ ^ a-b 7V-H0)(1 + F^^) - N-\0){1 + Fr) ' ^ " ^ 



Then Eq. (3-4) takes the form 

5H - iip5pp - iJinSpn = Xu^ + Yv'^ > (3-8) 

where 



X^-{ a + b + sign{c)yj{a-by + c^ ) 



(iVp(0) + iV„(0))- 



1 + ^0. ) (3-9) 



and 



y = -( a + & — sign(c)\j (a — 6)^ -|- ) 



defining the new generahzed Landau parameters Fq^ 



•s.a 



Hence, thanks to the rotation Eq.(3-6), it is possible to separate the total 
variation Eq.(3-3) into two independent contributions, the "normal" modes, 
characterized by the "mixing angle" which depends on the density of states 
and the details of the interaction. 

In the symmetric case, Np ^ = N, F^"" = F^^ and F^^ = F^, Eq.(3-4) 
reduces to 



SH - HpSpp - Pn5pn = ^ (1 + Fq){6Pp + 6pnf 



+ ^^^(l + ^o")(<^Pp- W (3-11) 

where Fq = F^"' + Fq^ and Fq = Fq^ — Fq^ are symmetric and antisymmetric 
(or isoscalar and isovector) Landau parameters, and the usual Pomeranchuk 
stability conditions for pure isoscalar /isovector fluctuations are recovered [51]. 

In the general case u- and f- variations, Eq.(3-6), can be interpreted as new 
independent isoscalar-like and isovector-like directions appropriate for asym- 
metric systems and F^^ and F^^, defined by Eqs. (3-9,3-10), can be considered 
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as generalized symmetric and antisymmetric Landau parameters. Now, the 
instability of the system can be studied completely just looking at the curva- 
tures along the u and v directions. If the system is stable against these two 
directions, it cannot be unstable in any other directions. On the contrary, if 
the system in unstable against u or v direction, or both, it can be unstable 
also in other directions of the {Spn, Spp) plane. 

Thus the thermodynamical stability requires X > and F > 0. Equivalently, 
the following conditions have to be fulfilled: 

1 + Fo^>0 and 1 + > 0, (3-12) 



They represent Pomeranchuk stability conditions extended to asymmetric bi- 
nary systems. 

As one can intuitively expect, the new stability conditions, Eq.(3-12), are 
equivalent to mechanical and chemical stability of a thermodynamical state, 
[73], i.e. 



where P is the pressure and y the proton fraction. In fact, mechanical and 
chemical stability are very general conditions, deduced by requiring that the 
curvatures of thermodynamical potentials, such as the free energy (or the 
entropy) with respect to the extensive variables are positive (negative). In the 
case discussed here, it can be proven that [20]: 



XY = N-\0)N-\0)[{1 + FDil + - F^^Fn 



(1 -y)p2 yopj 



T,y 



dy 



(3-14) 



T,P 



and: 



oc X{\/tcos[5 -\- -j^sinfiY + Y{Vtsinf3 — -j^cosfiY 



with t^—^. (3-15) 
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Moreover, from Eq.3-15 (first line), it is possible to see that {dP/dp)T,y is 
proportional to the variation of the energy along the direction dpn/Spp = 1/t = 
(1 — y)/y = Pn/ Pp- In fact, from Eq.(3-4), along the direction 6pp = t 6pn, one 
has: 5H — Pp6pp — Pn^Pn = \t {at + b/t + c)5p^. Thus the sign of {dP/dp)T,y 
gives information about the stability (or the instability) of the system against 
variations that preserve the initial proton to neutron density ratio. However, 
as we have seen, this is not the isoscalar-like direction of the normal mode, 
that is given by u. 



3.1.1 Isoscalar-like and isovector-like instabilities 

From Eq.s (3-8,3-9,3-10,3-12), one can define as isoscalar-like instability the 
case when the state is unstable against m— fluctuations, i.e. when 1 -|- Fq^ < 
(or X < 0). The name isoscalar-like comes from the fact that since the mixing 
angle is in the interval < /3 < 7r/2 the neutron and proton oscillations are in 
phase in the normal-mode m— direction. Analogously we deal with isovector- 
like instability when the system is unstable against fluctuations i.e. when 
l + F^g < (or F < 0). 

Let us consider first the case when X < 0. Since, as seen before, the me- 
chanichal instability {{dP/dp)T,y < 0) is not along the normal-mode -u-direction, 
when the isoscalar-like instability starts to appear (X = 0), it cannot be a me- 
chanical instability, so it corresponds to a chemical instability. On the contrary, 
when the isoscalar-like instability becomes stronger, we will have mechanical 
instabilities and the chemical instability will in turn disappear. However the 
nature of the unstable mode has not changed. So the distinction between me- 
chanical and chemical instabilities is purely semantic. This was shown for the 
first time in ref.[20] for nuclear matter and confirmed in the case of finite 
nuclear systems in ref. [21]. 

It is also interesting to observe that when {dP/dp)T,y changes the sign, passing 
through zero, the quantity (^)^p changes also the sign, but passing through 

infinity. In fact, from Eq.(3-14), one sees that the product (^) (f^) oc 
XY is a finite negative number, since X <0 and Y > 0. 

The case of symmetric nuclear matter is easily recovered. In fact, as expected, 
now the isoscalar instability, {X < 0,Y > 0), appears as mechanical instabihty 
and the isovector instability, (X > 0,y < 0), as chemical instability. Indeed 
we have t — 1, a — b, (3 — tt/A and so X and (f^)^ are proportional (see 
Eq.(3-15)) in this case. 

A more quantitative analysis can be performed in the case of asymmetric 
nuclear matter considering that the quantities a and b remain positive. In this 
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way one can study the effect of the interaction between the two fluids, given by 
c, on the instabihties of the mixture. If c < 0, i.e. for an attractive interaction 
between the two components, from Eq.(3-10) one sees that the system is stable 
against isovector-hke fluctuations but it becomes isoscalar unstable if c < 
—2\fab (see Eq.(3-9)). However thermodynamically this instabihty against 
isoscalar-like fluctuations will show up as a chemical instability if {—ta—b/t) < 
c < —2\/ab or as a mechanical instability if c < (—ta — b/t) < — 2\/a6 (see 
Eq.(3-15)). So the distinction between the two kinds of instability (mechanical 
and chemical) is not really relevant since the nature of density fluctuations is 
essentially the same, i.e. isoscalar-like. If c > 0, i.e. when the interaction 
between the components is repulsive, the thermodynamical state is always 
stable against isoscalar-like fluctuation, but can be isovector unstable if c > 
2\/ab. Since the system is mechanically stable (a, 6, c > 0, see Eq. (3-15)), the 
isovector instability is now always associated with chemical instability. Such 
situation will lead to a component separation of the liquid mixture. Following 
this line a complete analysis of the instabilities of any binary system can be 
performed, in connection to signs, strengths and density dependence of the 
interactions. 



3.1.2 Nuclear matter instabilities 



We show now quantitative calculations for asymmetric nuclear matter which 
illustrate the previous general discussion on instabilities. Let us consider a 
potential energy density of Skyrme type, [50,51], 



H (n ^ _ -4 (Pn + P,f . B (p„ + p,)"+^ 

+ (C, _ C2{-^T) (3-16) 
Po Po 

where po = 0.16 fm^'^ is the nuclear saturation density. The values of the 
parameters A = -356.8 MeV, B = 303.9 MeV, a = 1/6, Ci = 125 MeV, = 
93.5 MeV are adjusted to reproduce the saturation properties of symmetric 
nuclear matter and the symmetry energy coefficient. 



We focus on the low density region, where phase transitions of the liquid-gas 
type are expected to happen, in agreement with the experimental evidences 
of multifragmentation [74,75]. Since a, 6 > and c < 0, we deal only with 
instability against isoscalar-like fluctuations, as for symmetric nuclear matter. 
In Fig. 3-1 the circles represent the spinodal line corresponding to isoscalar- 
like instability, as deflned above, for three values of the proton fraction. For 
asymmetric matter, y < 0.5, under this border one encounters either chemical 
instability, in the region between the two lines, or mechanical instability, under 
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Fig. 3-1. Spinodal line corresponding to isoscalar-like instability of asymmetric nu- 
clear matter (circles) and mechanical instability (crosses) for three proton fractions: 
y = 0.5 (a), y = 0.25 (b), y = 0.1 (c). 

the inner line 



crosses). The latter is defined by the set of values (p, T) for 
= 0. We observe that the line defining chemical instability 

is more robust against the variation of the proton fraction in comparison to 
that defining mechanical instability: reducing the proton fraction it becomes 
energetically less and less favoured to break in clusters with the same initial 
asymmetry. 



However, we stress again the unique nature of the isoscalar-like instability. The 
change from the chemical to the mechanical character along this border line 
is not very meaningful and does not affect the properties of the system. Later 
we will show a quantitative realistic case of the clusterization of a dilute and 
heated asymmetric nuclear matter in a box, Subsect. 3.3. As seen before, in 
1.50 




0.00 



0.00 0.02 0.04 0.06 0.08 0.10 

p(fm-^) 



Fig. 3-2. Density dependence of the function see text, for three proton fractions, 
y = 0.5 (solid), y = 0.4 (open circles), and y = 0.1 (full circles) at T = IMeV. 

ANM isoscalar-like fluctuations are always associated with a chemical effect 
(change of concentration), which is responsible for the Isospin Distillation in 



31 



phase transitions [74] . Indeed the variation of the asymmetry (7 = 1 — 2y) — 

(Pn- Pp)/{pn + Pp) is 

SI = Sppi^^^^M - 1] ^ 6p,[x - 1] (3-17) 

where Iq is the initial asymmetry. The defined function x{p) is reported in 
Fig. 3-2 for different asymmetries. For y = 0.5, SI = 0, but for an initial 
proton fraction y < 0.5, one finds x < 1. Therefore 5/ < if 6pp > 0, thus 
when the density increases (liquid phase), the asymmetry decreases. 

The effect of the isospin distillation can be connected to the strength of the 
symmetry energy and, in particular, to the derivative of the symmetry energy 
coefficient with respect to p. This is easy to demonstrate if we adopt the sim- 
plest form for the symmetry energy: Egym = ^C-^P = \Csym{p)I'^- According 
to Eq.(3-7), the ratio 8pp/5pm can be expressed as: 5pp/5pn — tg(5' — l/tg(3, 
with tg{2P') = —tg{2(5) = —c/{a — h). The difference (a — h) can be written 
as: {a-h) = {N'^ - A^'^^) + (N-^F^^ - A^-^Fq""). For the simple interac- 
tion considered, the second term of the sum vanishes and (a — b) is equal to 
the difference of the inverse of proton and neutron single-particle level den- 
sity, which is positive in n-rich matter. The quantity c is the derivative of the 
proton potential Up — ^ (see Eq.(3-2)) with respect to the neutron density. 
This is equal to the sum of a negative term coming from the isoscalar part of 
the considered effective interaction and a negative term {—C/po), that comes 
from the symmetry energy. Hence the term — c is positive and increases when 
C increases. In this case tg2l3' becomes larger and consequently Spp/6pn in- 
creases leading to a larger distillation effect. One can easily notice that C/po 
coincides with the derivative of Csym with respect to p, hence the distillation 
effect is related to the derivative of the symmetry energy coefficient. 

Within this simple model even the effect of the Coulomb repulsion on the 
distillation effect can be easily understood. Taking into account the Coulomb 
interaction, the term {N~^Fq^ — N~^Fq^) does not vanish anymore, but it 
is equal to the derivative of the Coulomb potential with respect to the pro- 
ton density, that is a positive term. Hence the difference (a — b) increases 
and consequently the distillation effect decreases when including the Coulomb 
interaction, as intuitively expected. 

Investigations on instabilities in nuclear matter can be extended also to the 
high density region. Fig.3-3 shows the density dependence of the generahzed 
Landau parameters Eqs.(3-9), (3-10). For the interaction considered here the 
system exhibits another instability at high density, around 1.5/m^"^ (far away 
from the validity of this NM model), where the quantity c is positive and 
1 + F^g < (or y < 0). From Eq.(3-14) one finds that this is again a chemical 
instability. However now it results from isovector-like fluctuations, in contrast 
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Fig. 3-3. Density dependence of the generalized Landau parameters for two proton 
fractions, y = 0.4 (a) and y = 0.1 (b) (symmetric, solid, and antisymmetric, dashed) 
at T = IMeV. 

to the low density instability. The reason is the change in the character of 
the interaction between the two components. Since the interaction becomes 
repulsive the nuclear phase can become unstable against proton-neutron sep- 
aration. We also notice in Fig. 3-3 that the generalized Landau parameters 
display a discontinuity where the quantity c changes the sign. Other effective 
forces, with more repulsive symmetry terms, will not show this high density 
chemical instability [50] , that actually could be of interest for other many body 
systems. 



3.2 Dynamical analysis 



The dynamical behaviour of a two-fluid system can be described, at the semi- 
classical level, by considering two Vlasov equations, for neutrons and protons 
in the nuclear matter case [76,77,50,51], coupled through the self-consistent 
nuclear field : 

dfq{r,p,t) p dfg dUqirJ) dfg 

P ■ ^ Z ^ = 0' q = n,p. (3-18) 

at m or or op 



For simplicity effective mass corrections are neglected. In fact in the low den- 
sity region, of interest for our analysis of spinodal instabilities, effective mass 
corrections should not be large. 

Ug{r, t) is the self-consistent mean field potential in a Skyrme-hke form [50,51] : 
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U.= ^-^ = A(P-]+B(P-]"\c((^]r, 



^Pq \PoJ \PoJ \P0 

2 dp po 



where 

HMPn, Pp)-W- + Z^^ + ^- + ^(Vp)' - ^(Vpa)' (3-20) 



2 Po a + 2 p^+^ 2 Po 2^ ^' 2 



is the potential energy density (see Eq.3-16), where also surface terms are 
included; p = Pn + Pp and Pz = Pn — Pp are respectively the total (isoscalar) 
and the relative (isovector) density; Tg — +1 {q — n), -1 {q — p)- 

The value of the parameter D = 130 MeV-fm^ is adjusted to reproduce the 
surface energy coefficient in the Bethe-Weizsacker mass formula agurf = 18.6 
MeV. The value = 40 MeV-fm^ ~ D/3 is chosen according to Ref. [78], 
and is also close to the value — 34 MeV-fm^ given by the SKM* interaction 
[3]. 

Let us now discuss the linear response analysis to the Vlasov Eqs. (3-18), cor- 
responding to a semiclassical RPA approach. For a small amplitude pertur- 
bation of the distribution functions fq{r,p,t) , periodic in time, Sfq{r,p,t) ~ 
exp{—iu!t), Eqs. (3-18) can be linearized leading to the following form: 

.,^Sf, + ^'-^-'Jffl-'-^?f = 0, (3-21) 
m or or op or op 



where the superscript (0) labels stationary values and 6Uq is the dynamical 
component of the mean field potential. The unperturbed distribution function 
f^°^ is a Fermi distribution at finite temperature : 

/f (4) = M ^ w^^i ■ (3-22) 

" " exp (e^ - pq)/T + 1 



Since we are dealing with nuclear matter, Vrf/q°^ = in Eq. (3-21) and 6fq oc 
exp(— ia;t -|- ikr). Following the standard Landau procedure [72,50], one can 
derive from Eqs. (3-21) the following system of two equations for neutron and 
proton density perturbations : 

[1 + FrXn]6pn + K''Xn]6pp = , (3-23) 
[FrXp]Spn + [1 + FS''Xp]5pp = , (3-24) 
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where: 

Xq{^, k) 



kv a/f 

Ng{T)J (27rn)3u; + i0-kv de$ ' 



2 dp 



1 /■2 c 
?2^ 



(3-25) 



is the long- wavelength hmit of the Lindhard function [72], y — p/m and 



(3-26) 



are the usual zero-order Landau parameters, as already introduced in Eq.(3- 
2), where now the /c-dependence is due to the presence of space derivatives in 
the potentials (see Eq.(3-19)). For the particular choice of potentials given by 
Eq.(3-19), the Landau parameters are expressed as: 



Po Po Po 



12 



dCp^ 
dp Po' 



+ + 



dp^ 2po 



.(3-27) 



Multiplying the first equation by A'^ ^Xp ^^'^ the second one by ^Xn, we 
are led to define the following functions: 



a{k,u;) = N^\l + Fl'\,)xn ; bik,u) = N-\l + FrXn)Xp ; 

c{k, u) = {N-'Fr + N-'F-^)xnXv = 2N-'FrXnXp, (3-28) 

in some analogy with Eqs.(3-5) and we obtain the following system of equa- 
tions: 



a6pp + c/2 6pn = 0; 

c/2 dpp + bdpn = (3-29) 

The system can be diagonalized with eigenvalues and Aj solutions of the 
equation: 

(a-A,,i)(6-A,,,)-cV4 = . 

Formally wc obtain for Xg i the same expressions as given in Eqs. (3-9,3-10) 
for X and Y, but now a, b and c depend on to. The unstable solutions for u 
are obtained by solving the equations: A^ = (for isoscalar-hke fluctuations), 
Aj = (for isovector-like fluctuations). This problem is completely equivalent 
to solve the equation: c^(a;, k) — 4a{u!, k)b{uj, k), i.e. the dispersion relation 

(1 + FrXn){l + FE'Xv) - F^FrXnXp = , (3-30) 
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that is also obtained directly by imposing the determinant of the system of 
Eqs.(3-23), (3-24) equal to zero. 

The dispersion relation is quadratic in cu and one finds two independent solu- 
tions (isoscalar-like and isovector-like solutions): uj^ and cuf. Then the structure 
of the eigenmodes can be determined and one finds: 

Spp/Spn = -26(0;^, k)/c{ujs, k), 

for the isoscalar-like modes and 

Spp/Spn = -2b{ui, k)/c{uji, k), 

for isovector-like oscillations. However, it is important to notice that the cor- 
responding angles are not equal to the angle /? determined in the thermo- 
dynamical analysis, Eq.(3-7), because of the cu dependence in a, b and c. They 
only coincide with (3 when cu — (and thus Xn,p — !)> i-e. at the border of the 
unstable region. 

3.2.1 Illustrative results for liquid-gas phase transitions 

The dispersion relation, Eq.(3-30), have been solved for various choices of the 
initial density, temperature and asymmetry of nuclear matter. Fig. 3-4 reports 
the growth rate F = Im u}{k) as a function of the wave vector k, for three 
situations inside the spinodal region. Results are shown for symmetric (/ = 0) 
and asymmetric (/ = 0.5) nuclear matter. The growth rate has a maximum 
Fo = 0.01 -7- 0.03 c/fm corresponding to a wave-vector value around ko ~ 
0.5 1 fm~^ and becomes equal to zero at /c ~ 1.5/co, due to /c-dependence of 
the Landau parameters, as discussed above. One can see also that instabilities 
arc reduced when increasing the temperature, an effect also present in the 
symmetric N = Z case [79,80,23]. 

At larger initial asymmetry the development of the spinodal instabilities is 
slower, decrease of the maximum of the growth rate. One should expect also 
an increase of the size of the produced fragments, decrease of the wave number 
corresponding to the maximum growth rate. From the long dashed curves of 
Fig. 3-4 we can predict the asymmetry effects to be more pronounced at higher 
temperature, when in fact the system is closer to the boundary of the spinodal 
region. 

3.2.2 Coulomb effects on instabilities 

The influence of Coulomb effects on the growth rates can be easily investigated 
within the formalism outlined above, as done in [81]. It suffices to add to the 
energy density, Eq.(3-20), the Coulomb energy density, that can be calculated 
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Fig. 3-4. Growth rates of instabilities as a function of the wave vector, as calculated 
from the dispersion relation Eq.(3-30), for three situations inside the spinodal region. 
Lines are labelled with the asymmetry value /. The insert shows the asymmetry of 
the perturbation dpi/Sps, as a function of the asymmetry / of the initially uniform 
system, for the most unstable mode, in the case p = 0.4po, T = 5 MeV. 

in the Hartree-Fock approximation, with the Fock term evaluated in the local 
density approximation: 

H'-' (r) = :f p,(r) / dr'iM^^ - 5 (i) . (3-31) 




This only modifies the Fq^ Landau parameter, implemented by adding the 
term: 

^-^y-^V' (3-32) 



Fig. 3-5 shows a comparison between the results obtained with and without 
the Coulomb interaction, for density p = 0.4po and T = 0. The Coulomb 
force causes an overall decrease of growth rates. This decrease is almost in- 
dependent of the asymmetry. Moreover, as discussed in [81], it depends only 
slightly on the temperature. It is also observed that, when the Coulomb force 
is included, the wave vector k must exceed a certain value kmin in order to 
observe instabilities. The two effects (the decrease of the growth rate and the 
appearance of kmin) are due to the competition between the Coulomb and the 
nuclear forces. In fact, the Coulomb forces push the protons towards regions 
of lower density, the nuclear forces instead push the neutrons in this direction 
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Fig. 3-5. Growth rates of the unstable modes for T=0, p = 0.4 pQ and three different 
values of the asymmetry parameter / (from top to bottom / = 0,0.3,0.6). Results 
including the Coulomb interaction (full line) and without the Coulomb interaction 
(dashed) are shown. Taken from [81] 

(neutron distillation). 



3.2.3 Isospin distillation 

A better understanding of liquid-gas phase transitions in a two-component 
system can be achieved by studying the chemical composition of the growing 
mode. This is shown in the insert of Fig. 3-4, where the asymmetry of the per- 
turbation Ipt = {6pn — Spp)/ (ySpn + Spp) = Spi/Sps, as obtained by solving the 
system of Eqs.(3-23), (3-24), is diplayed as a function of the asymmetry of the 
initially uniform system / = (p^°) — p^p^)/{Pn^ + Pp"-*), for the dominant (most 
unstable) mode (for the case p — O.A po, T — 5 MeV). Without any chemical 
processes it should be Ipt — I. However one obtains Ipt < 0.5 /. This means 
that a growing mode produces more symmetric high-density regions (liquid 
phase) and less symmetric low-density regions (gas phase). Hence, during the 
fragmentation, a collective diffusion of protons from low-density regions to 
high-density regions takes place. This is the isospin fractionation (distillation) 
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effect already discussed in the context of liquid-gas phase transitions in two- 
component systems. This is essentially due to the increasing behaviour of the 
symmetry energy per nucleon with density, in the density region considered 
here. The distillation effect is represented in Fig. 3-6 , where the ratio 5pn/5pp, 
as obtained by solving the dispersion relation, is displayed (direction and mag- 
nitude of the arrows) for some points of the (pp, pn) plane inside the spinodal 
region. One observes that the distillation effect is more pronounced at large 
asymmetry (i.e. smaller proton fraction y). 

We recall that the Isospin Distillation has been revealed even in n-rich finite 
nuclei, performing quantal RPA calculations [21] (see also [23] for a thorough 
discussion on instabilities in finite systems). 




Fig. 3-6. Direction of 6pn/Spp in several points of the {pp,Pn) plane. 



As discussed before, the size of the effect is related to the stiffness of the sym- 
metry energy, as illustrated in Table 5, where the ratio Spn/^Pp, as obtained 
using a stiff or a soft (in brackets) parameterization of the symmetry energy, is 
reported. As expected, the neutron distillation effect (reduction of the asym- 
metry in the formed clusters) is systematically larger at lower densities in the 
asy-soft case while the opposite is seen already at p = 0.4po- This is nice evi- 
dence of the dependence on the slope of the symmetry term in the low density 
region (see Fig.l, Sect.l). 
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I = 0.2 {N/Z = 1.5) 


/ = 0.8 {N/Z = 9.0) 


0.2 


1.23 (1.2) 


3.15 (2.61) 


0.3 


1.19 (1.186) 


2.44 (2.4) 


0.4 


1.15 (1.18) 


1.99 (2.32) 



Table 5 

Isospin content of the clusters, density variations Spn/Spp, formed in dilute matter 
at different initial densities for two initial asymmetries, N/Z = 1.5 and 9.0. The 
values in brackets are obtained with a softer symmetry term (asy-soft EOS) 

3.3 Simulation results: heated nuclear matter in a box 

The previous analytical study is restricted to the onset of Fragmentation, 
and related Isospin Distillation, in Nuclear Matter, in a linearized approach. 
Numerical calculations have been performed in order to study all stages of the 
fragment formation process [51,85]. We report on the results of ref.[51] where 
the same effective Skyrme interactions have been used. 

In the numerical approach the dynamical response of nuclear matter is studied 
in a cubic box of size L imposing periodic boundary conditions. The Landau- 
Vlasov dynamics is simulated following a phase-space test particle method, 
using gaussian wave packets [86,87,88]. The dynamics of nucleon-nucleon col- 
lisions is included by solving the Boltzmann-Nordheim collision integral using 
a Monte-Carlo method [87]. The width of the gaussians is chosen in order to 
correctly reproduce the surface energy value in finite systems. In this way a 
cut-off appears in the short wavelength unstable modes, preventing the for- 
mation of too small, unphysical, clusters [79] . The calculations are performed 
using 80 gaussians per nucleon and the number of nucleons inside the box is 
fixed in order to reach the initial uniform density value. An initial temperature 
is introduced by distributing the test particle momenta according to a Fermi 
distribution. 

We have followed the space-time evolution of test-particles in a cubic box with 
side L — 24/m for three values of the initial asymmetry / = 0, 0.25 and 0.5, 
at initial density p*^°^ = 0.06/m~^ ~ 0.4po and temperature T = 5 MeV. 
The initial density perturbation is created automatically due to the random 
choice of test-particle positions. Results for the initial asymmetries / = and 
/ = 0.5, are reported in Fig. 3-7, (a) and (b) respectively. The figure shows 
density distributions in the plane z — 0, which contains the center of the 
box, at three time steps t = 0, 100 and 200 fm/c, corresponding respectively 
to initial conditions, intermediate and final stages of the cluster formation. 
Clearly, a growth of the small initial density perturbations takes place, this 
time non-linear and hard particle collision effects are included. 
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(b) 




Fig. 3-7. Time evolution {t = 0.0, 100.0 and 200.0 fm/c) of the density p{x,y) in 
the plane z = 0. Initial conditions: p = 0.06/m~^, T = 5MeV and asymmetries 
/ = 0.0 (a); / = 0.5 (6). Upper panels: contour plots; Lower panels: corresponding 
two-dimensional surfaces. Densities in fm~^. 



We can compare the dynamical evolution with the analytical predictions of the 
previous sections. To do this, two variables are constructed: the total density 
variance (see [82]) 



a=< (p-p(°)r> 



all 



(3-33) 



and the correlation function between proton and neutron density perturba- 
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tions, normalized to the neutron density variance, 



< {P p - Pf^){pn - P'^n'^) >all 



^ KFp- Hp jyHn- Hn J -^all looA\ 
iipn - — (5)—- ■ i'3-'34j 



In Eqs. (3-33), (3-34) < ... >aii denotes the average over aU test particles, 
while < ... >n denotes the average over neutrons only. The densities p, p„ 
and pp are calculated in the position of the test particle considered by taking 
contributions from gaussians of all test particles. For a dominant plane-wave 
perturbation we have the limit: 

a oc exp(2n), Rp^ = ^ . (3-35) 



Fig. 3-8 shows the evolution of a (a) and of the (tcst-particle) perturbation 
asymmetry Ipt = (1 — -Rp„)/(1 + Rpn) (b) for the same initial conditions dis- 
cussed above, i.e. T = 5 MeV , = 0.4po and asymmetries / = 0.0, 0.25, 0.5. 
A general feature is the clear linear increase of ln{a) in the time interval 




Fig. 3-8. Time evolution of the density variance (a) and of the perturbation asymme- 
try (b) (see text) for the box simulations at initial density p = T = 5MeV 
and asymmetries I = 0.0, 0.25, 0.5. The straight lines in (a) show linear logarithmic 
fits in the initial stage of the cluster formation. 

50 < t < 150 fm/c. During the first 50 fm/c the system is quickly "self- 
organizing" , selecting the most unstable normal mode. Afterwards the variance 
(Eq. (3-33)) increases exponentially with a time scale given by T = Im uj{k), 
finally it saturates. In correspondence (sec Fig.3-8b), the perturbation asym- 
metry Ipt reveals also a quick saturation at t ~ 50 fm/c. At times before 
50 fm/c the proton and neutron density perturbations are not correlated 
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{Ipt{t = 0) ~ 1.0 since nicely Rpn{t = 0) ~ 0.0) , but at t > 50 fm/c the 
correlation of plane-wave type {Spp/Spn — const > 0) develops. 

We notice that the time scales necessary to reach the asymmetry value char- 
acteristic of the most important growing modes, which are quite short in 
our calculations, generally depend on the structure of the initial noise put 
in the neutron and proton densities. In our calculations all modes are nearly 
equally excited. This causes the quick appearance of the features associated 
with the dominant mode. In agreement with analytical calculations, the in- 
stability grows slower in the case of larger asymmetry. 

For an initial asymmetry / = 0.5, the extracted values of growth time F ~ 0.01 
c/fm and perturbation asymmetry Ipt ~ 0.24 (see Fig.3-8), and of wavelength 
A ~ 12 fm (from the distance between the density distribution maxima in 
Fig. 3- 7b), are in good agreement with the analytical results presented before. 

The Spinodal Decomposition Mechanism, SD, leads to a fast formation of the 
liquid (high density) and gaseous (low density) phases in the matter. Indeed 
this dynamical mechanism of clustering will roughly end when the variance 
(Eq. (3-33)) saturates [83], i.e. around 250 fm/c in the asymmetric cases (see 
Fig.3-8a). We also discuss the "chemistry" of the liquid phase formation. In 



t=0 fm/o 



t=50 fm/o 



1=100 tm/o t=150 fm/c 



t=250 fm/c 
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Fig. 3-9. Time evolution of neutron (thick lines) and proton (thin lines) abundances 
(a) and of asymmetry (b) in different density bins. The calculation refers to the case 
of T = 5MeV, with initial average density p = 0.06/m~^ and asymmetry I = 0.5 
(see the first panel of the (b) plots) . 

Fig. 3-9 we report the time evolution of neutron (thick histogram in Fig. 3- 
9a) and proton (thin histogram in Fig.3-9a) abundances and of asymmetry 
(Fig.3-9b) in various density bins. The dashed lines respectively shows the 
initial uniform density value p ~ 0.4po (Fig.3-9a) and the initial asymmetry 
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/ = 0.5 (Fig.3-9b). The drive to higher density regions is clearly different for 

neutrons and protons: at the end of the dynamical clustering mechanism we 
have very different asymmetries in the liquid and gas phases (see the panel at 
250/m/c in Fig.3-9b). 

It was shown in Refs. [78,19,84], on the basis of thermodynamics, that the two 
phases should have different asymmetries, namely, Igas > Iiiquid, and actually 
a pure neutron gas was predicted at zero temperature if the initial global 
asymmetry is large enough (/ > 0.4) [78]. Here we are studying this chemical 
effect in a non-equilibrium clustering process, on very short time scales, and 
we confirm the predictions of a linear response approach discussed before. 

We can directly check the important result on the unique nature of the most 
unstable mode, independent of whether we start from a mechanical or from 
a chemical instability region. The isospin distillation dynamics presented in 
the Fig. 3-9 refers to the initial conditions of T = 5MeV, average density 
p = 0.06/m~^ and asymmetry / = 0.5, i.e. we start from a point well inside 
the mechanical instabihty region of the used EOS, see Fig.3-l(b). We can 
repeat the calculation at the same temperature and initial asymmetry, but 
starting from an initial average density p = 0.09/m~^, i.e. inside the chemical 
instability region of Fig. 3- 1(b). The results for the Isospin Distillation Dynam- 
ics are shown in the Fig. 3-10. The trend is the same as in the previous Fig. 3-9. 

t=0 fm/o t=100 fm/e t=200 fm/o t=250 fm/c t=300 fm/c 




0.0 0.1 0.0 0.1 0.0 0.1 0.0 0.1 0.0 0.1 



p (fm-») 




Fig. 3-10. Same calculation as in Fig.3-9 but with initial average density 
p = 0.09/m~^, inside the chemical instability region. 

This nicely shows the uniqueness of the unstable modes in the spinodal insta- 
bility region, as discussed in detail in the previous subsection. Such result is 
due to gross properties of the n/p interaction, thus it should be not dependent 
on the use of a particular effective force. This has been clearly shown recently 
in the linear response frame, [22], and in full transport simulations, [85]. 
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As intuitively expected, and as confirmed by the RPA analysis (see [51]), 
the isospin distillation effect becomes more important when increasing the 
initial asymmetry of the system. At the same time, the instability growth 
rates become smaller for the more asymmetric systems, see Fig. 3-4. 

Moreover it is possible to observe a rather smooth and continuous transition 
from the trend observed at p = (mechanical unstable region) to the 

trend observed at p = 0.09/m~^ (chemical unstable region), thus indicating 
that there is no qualitative change between the two kinds of instabilities. In 
fact they actually correspond to the same mechanism, the amplification of 
isoscalar-like fiuctuations, with a significant chemical component (change of 
the concentration). 

The conclusion is that the fast spinodal decomposition mechanism in neutron- 
rich matter will dynamically form more symmetric fragments surrounded by 
a less symmetric gas. Some recent experimental observations from fragmen- 
tation reactions with neutron rich nuclei at the Fermi energies seem to be in 
agreement with this result on the fragment isotopic content : nearly symmet- 
ric Intermediate Mass Fragments (IMF) have been detected in connection 
to very neutron-rich light ions [74], [75]. This will be extensively discussed in 
Sect. 5. 



45 



4 Symmetry term effects on fast nucleon emission and collective 
flows 

4.1 Pre- equilibrium dynamics 

The early reaction dynamics is mainly governed by the pressure of the excited 
nuclear matter formed in the initial stage of the collision, [24,25,26]. We recall 
that a stiff symmetry energy gives a larger gradient of pressure than the 
asy — soft shown in Sect. 2, Eq.(8). This should result in an overall 

faster emission of particles. Moreover, the soft behaviour of the symmetry 
energy leads (at low density) to a larger repulsion of neutrons, with respect to 
the stiff case (see Fig.2-3 of Sect. 2). Thus neutrons are expected to be emitted 
at earlier times, with respect to protons, in the asy — soft case. The isotopic 
content of the pre-equilibrium emission is also sensitive to the symmetry part 
of the EOS. For the same reasons explained above more neutrons are emitted 
in the soft case. 

These ideas have lead to investigate the effect of the symmetry energy on 
the nucleon emission. Of course we do not have direct access to the particle 
emission time experimentally, but two-particle correlation functions, through 
final-state interactions and quantum statistics effects, have been shown to be 
a sensitive probe to the temporal and spatial distribution of emission sources 
during the reaction dynamics of heavy ion collisions at incident energies rang- 
ing from intermediate [89,90,91] to RHIC energy [92]. 

To discuss isospin effects on pre-equilibrium emission, we will review the results 

obtained by Chen et al., [93,94], using an isospin dependent transport code, 
IBUU . This is a one-body dynamics approach without fluctuating terms [149] 
but the results discussed here can be reliable since mainly due to the average 
nuclear dynamics. 

We will consider a typical reaction with radioactive beams: ^^Ca -|- ^^Ca aiE = 
80 MeV/nucleon [95], central collisions. This study has been performed for two 
kinds of density dependence, labelled as "soft" and "stiff", that correspond 
to 7 = 0.5 and 7 = 2.0 respectively, in the parametrization of the symmetry 
energy (total, kinetic -|- potential) 

-Esym(p) = -Esym(po) " u' ■ 

It is important to note that, in this case, "soft" refers to a slope around normal 
density L = 52 MeV, which is close to the common value, while in the stiff 
case the slope is L = 210 MeV, which is much larger than what has been 
labelled as "super-stiff" (L ~ 100 MeV) in the previous Sect. 2. 
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4.1.1 Average emission times 



For the two symmetry energy parametrizations we report in Fig. 4-1 {soft: 
squares, stiff: triangles) the average emission time of protons and neutrons 
as a function of their momenta. It is seen that the average emission time of 
nucleons with a given momentum is carher for the stiff symmetry energy than 
for the soft one. This is expected because of the overall larger pressure gra- 
dient in the stiff case. We also notice that the relative emission sequence of 
neutrons and protons is, however, determined by the difference in their sym- 
metry potentials. In other words if we average over neutron and proton the 
emission time is faster for the stiff case, but, on top of this there is a differ- 
ence in the neutron-proton emission that is driven by the symmetry potential 
which is more attractive for protons causing an emission delay with respect 
to neutrons. 

The neutron-proton difference is larger for the soft case because the isovector 
mean field at low densities, where most nucleons are emitted, is larger for the 
soft symmetry energy than for the stiff symmetry potential. In such a reasoning 
we have disregarded the effect of Coulomb interaction which reduces the delay 
of protons with respect to neutrons. For ^^Ca-|-^^Ca the effect is quite small, 
but it grows with the size of the colliding system modifying the p — n difference 
in the emission time, \tp — Therefore for heavier systems one can expect a 
reduction of the \tp — t„| for the soft case and an increase of the difference for 
the stiff case, because the stronger Coulomb effect will shift down the proton 
average emission time. We will discuss the consequences on the correlation 
functions and the cluster formation in the corresponding subsections. 

We notice the relation between larger momenta and earlier emission. For the 
reaction considered here, it is shown that most of the particles with momem- 
tum P > 250 MeV/c are emitted within t < 70 fm/c. This allows to connect 
the momenta of the emitted particles with the time evolution of the reaction, 
giving hints also about the average density at which particles are emitted. 
The crossing point between the neutron-proton emission time is connected to 
the crossing in the symmetry potential, including Coulomb, at density around 
O.Spo- We can say, for example, that particles with momenta of the order of 
300 MeV show features of the interaction already at density around and below 
the saturation one. 

Using a coalescence model also the emission time of light clusters, such as d, t, 
^He has been studied [96]. The authors have found that light cluster produc- 
tion occurs mainly from the pre-equilibrium stage and has a shorter duration 
than nucleon emission. We notice that this is different from the cluster forma- 
tion at lower energies, where the multifragmentation process is dominant, see 
the previous Sect. 3. A stiffer symmetry energy causes an earlier emission of 
nucleons that reflects the earlier formation of clusters. The average emission 
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Fig. 4-1. Average emission times of protons and neutrons as functions of their mo- 
menta for different symmetry energies, for ^^Ca+^^Ca at 80 MeV/nucleon b=0 fm, 
see text (from [94]). 

time decreases with increasing kinetic energy as for the nucleons, moreover 
the heavier is the cluster and earher it is formed. This has a simple explana- 
tion: the probability to coalesce decreases with the average density, moreover 
we can have a cluster break-up in the expansion phase. The dependence is 
stronger for clusters with more nucleons. 

Therefore the study of cluster isobars has two advantages: they are easier 
to detect with respect to neutrons and they carry information from higher 
density. However the influence of the symmetry energy on the average emission 
time is reduced with respect to the nucleon one. This leads to a reduced 
sensitivity of observables like the correlation function of light clusters, see 
next subsections. 

In summary a simple analysis of the average emission time nicely shows 
how the particle emission is regulated by the symmetry potential. However, 
as mentioned above, the particle emission time cannot be directly accessed 
by experiments, but the nucleon emission function can be extracted from 
two-particle correlation functions, see, e.g., Refs. [89,90,91,92] for reviews. In 
most studies, only the two-proton correlation function has been measured 
[97,98,99,100,24]. In particular, the correlation function of two nonidentical 
particles has been found to depend on their relative space-time distributions at 
freeze out and thus provides a useful tool for measuring the emission sequence, 
time delay, and separation between the emission sources for different particles 
[101,102,103,104,105,106]. Recently, data on two-neutron and neutron-proton 
correlation functions have also become available. The neutron-proton corre- 
lation function is especially useful as it is free of correlations due to wave- 
function anti-symmetrization and Coulomb interactions. Indeed, Ghetti et al. 
have deduced from measured neutron-proton correlation function the time 
sequence of neutron and proton emissions [107,108]. 
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4-2 Nucleon-nucleon correlation functions 



In the standard Koonin-Pratt formalism [109,110,111], the two-particle corre- 
lation function is obtained by convoluting the emission function g(p,x), i.e., 
the probability for emitting a particle with momentum p from the space-time 
point X = (r, t), with the relative wave function of the two particles, i.e., 

r(r a) - /^^ ■'■i^^^^'-2^(P/2-Ak(P/2..r,) |o(q.r)|^ 
^ '"^^ Jd^x,g{I>/2,x,)Jd^X29{P/2,X2) ' ^"^ 



In the above, P(= pi -|- P2) and q(=|(pi — P2)) are, respectively, the total 
and relative momenta of the particle pair; and 0(q, r) is the relative two- 
particle wave function with r being the relative position, i.e., r = (r2— ri) — 
|(vi -|- V2)(t2 — ^i)- This approach has been very useful in studying effects of 
nuclear equation of state and nucleon-nucleon cross sections on the reaction 
dynamics of intermediate energy heavy-ion collisions [90]. The extension to 
the study of symmetry energy from a theoretical point of view is instead quite 
recent [93,94]. 

Shown in Fig. 4-2 are two-nucleon correlation functions gated on total mo- 
mentum P of nucleon pairs from central collisions of ^^Ca + ^^Ca at = 80 
MeV/nucleon. The left and right panels are for P < 300 MeV/c and P > 
500 MeV/c, respectively. Both neutron-neutron (upper panels) and neutron- 
proton (lower panels) correlation functions peak at g MeV/c. The proton- 
proton correlation function (middle panel) is, however, peaked at about g = 20 
MeV/c due to the strong final-state s-wave attraction, hut is suppressed at 
g = as a result of Coulomb repulsion and anti-symmetrization of two-proton 
wave function. These general features are consistent with those observed in 
experimental data from heavy ion collisions [107]. 

Since emission times of low-momentum nucleons do not change much with 
the different -Esym(p) used in IBUU model as shown in Fig. 4-1, two-nucleon 
correlation functions are not much affected by the stiffness of the symmetry 
energy. On the other hand, the emission times of high-momentum nucleons, 
which are dominated by those with momenta near 250 MeV/c, differ appre- 
ciably for the two symmetry energies considered here. Correlation functions 
of high-momentum nucleon pairs thus show an important dependence on nu- 
clear symmetry energy. Gating on nucleon pairs with high total momentum 
allows one to select those nucleons that have short average spatial separations 
at emission and thus exhibit enhanced correlations. The strength of the cor- 
relation function is stronger for the stiff symmetry energy than for the soft 
symmetry energy: about 30% and 20% for neutron-proton pairs and neutron- 
neutron pairs at low relative momentum q = 5 MeV/c, respectively, and 20% 
for proton- proton pairs at g = 20 MeV/c. 
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Fig. 4-2. Two- nucleoli correlation functions gated on total momentum of nucleon 
pairs using the soft (filled squares) or stiff (filled triangles) symmetry energy. Left 
panels are for P < 300 MeV/c while right panels are for P > 500 MeV/c (from 
[93]). 



The neutron-proton correlation function thus exhibits the highest sensitivity 
to variations of Esym{p), however one should heed that proton-proton corre- 
lation function is much easier to measure thanks to the proton charge and 
to the larger relative momentum at which the effect is expected. As shown 
in Fig. 4-1 and discussed earlier, the emission sequence of neutrons and pro- 
tons is sensitive to Esy^{p). With a stiff i?sym(p), high momentum neutrons 
and protons are emitted almost simultaneously, and they are thus temporally 
strongly correlated, leading to a larger neutron-proton correlation function. 
At variance in the case of soft Esym{p) we have a delay in the proton emission, 
larger attractive field, which reduces the temporal correlation with neutrons. 

Furthermore, both neutrons and protons are emitted earlier with stiff Esym{p), 
so they have smaller spatial separation. This demonstrates that correlation 
functions of nucleon pairs with high total momentum can indeed reveal sen- 
sitively the effect of nuclear symmetry energy on the temporal and spatial 
distributions of emitted nucleons. However, one must be aware that such an 
effect is delicate and can vanish as a function of impact parameter, beam 
energy, colliding system. 

The impact parameter dependence has been studied again for the ^^Ca + ^'^Ca 
at = 80 MeV/nucleon. It is seen that the symmetry energy effect is stronger 
in central and semi-central collisions and becomes weaker in semi-peripheral 
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Fig. 4-3. Dependence of nucleon-nucleon correlation functions on the incident energy 
for hight total momentum (P > 500 MeV) neutron-neutron (a) and neutron-proton 
(c) pairs with relative momentum q = 5 MeV/c and proton-proton (b) pairs with 
relative momentum g = 20 MeV/c by using Ko = 380 MeV and free JV-JV cross 
sections with the soft (filled squares) or stiff (open circles) symmetry energy, (from 
[94]). 

and peripheral collisions. This reflects the decrease in nuclear compression 
with the consequent change in the reaction mechanism. 

The dependence of nucleon-nucleon correlation functions on the incident en- 
ergy of heavy ion collisions is shown in Fig. 4-3. It is seen that the values of 
nucleon-nucleon correlation functions increase with increasing incident energy. 
This is understandable since nuclear compression increases with increasing in- 
cident energy, and nucleons are also emitted earlier, leading to a smaller source 
size. The difference due to the symmetry energy does not change significantly 
with the incident energy at least in the range 50 — 150 MeV/nucleon. 

A dependence on the mass of the colliding system has been already studied 
in Ref. [94] for central collisions. The authors have looked at ^^^Sn -|- ^^^Sn 
at = 80 MeV/nucleon. It is seen that, even if the isospin asymmetry is 
similar to ^^Ca -|- "^^Ca system, the effect is reduced by about a factor two 
for pp and np, while nn is less affected. This can be understood with the 
increased Coulomb potential which is expected to reduce the difference in 
Cpn- Considering the emission times in ^^Ca -|- ^^Ca as a reference, for heavier 
systems the Coulomb potential increases the correlation for the soft case but 
decreases it for the stiff case, where for ^^Ca -|- "^^Ca neutrons and protons 
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have almost equal emission times. 

In conclusion the symmetry term effect on nucleon nucleon correlation finc- 
tions shows up only if particular conditions of impact parameter and colliding 
systems are matched. An experimental plan to pursue the investigation of 
such an observable needs a careful choice of the system and the possibility 
to have a good centrality selection. Futhermore it is important to choose the 
various colliding systems in such a way that the isospin asymmetry is as much 
different as possibile while the total mass is unchanged [112,113]. Moreover 
this kind of comparative analysis is in general more reliable than the study 
of absolute values of correlation functions, which themselves are not perfectly 
reproduced hy BUU simulations already for isospin symmetric systems, see 
ref.[24]. 



4.2.1 Isospin Momentum Dependence 

During the preparation of this report, phenomenological potentials extracted 
from the Gogny effective interaction [114] or extensions of the Bertsch-Gale- 
Das Gupta (BGBD) [65,115] have been considered to improve IBUU and 
BNV transport codes, to include the effect of a momentum dependence in the 
isovector channel, namely a splitting of the effective masses of neutrons and 
protons, see the discussion in Sect. 2. A detailed study of particle emission with 
a momentum dependent mean field was firstly done in [67,68], but without 
including isospin effects. Generally the momentum dependence causes a faster 
and richer pre-equilibrium emission. This reflects in an increasing of the two 
nucleon correlation function and of the yield of the clusters formed. 

A first analysis of the effect of the momentum dependence on the particle 
emission has shown a weakening of the role of the symmetry energy leading 
to a weaker dependence of the correlation function on the symmetry energy 
[116]. In fact such a study has employed an effective interaction with a splitting 
of the effective masses as m*^{p) > m*{p), which leads to a reduction of the 
effective action of the symmetry energy at higher density. Indeed the difference 
of the neutron-proton potentials decreases with momentum [65,114], see Fig.2- 
8 of Sect. 2. However, if a momentum dependence in the isovector channel 
corresponding to a splitting in the opposite direction is considered, (p) < 
m*{p), a strong effect on the correlation function can be envisaged, as we will 
see in the results concerning collective flows discussed later in the Section. 

Since this is a stimulating open problem, directly related to the slope of the 
Lane Potential of Sect. 2, a systematic study of the beam energy dependence 
of the (n, p) correlation functions appears to be very promising. A related in- 
teresting observable would also be just the isotopic content of the free nucleon 
emissions and/or of light isobars, see following. 
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4-2.2 Light Clusters 



Another observable of interest for the Esym{p) search is the yield and the en- 
ergy spectrum of light clusters that differ in the isospin content. As mentioned 
above, the hght clusters (isobars) offer the possibility of having particles with 
different isospin that are charged and therefore can be detected in a much 
simpler way with respect to neutrons. Moreover, as already discussed, their 
production time is shorter and they represent good probes of the higher den- 
sity phase. 

Isospin effects on cluster production and isotopic ratios in heavy ion collisions 
have been previously studied using either the lattice gas model [119] or a 
hybrid oilBUU and statistical multifragmentation model [120]. These studies 
are, however, at lower energies, where multifragmentation processes, as a result 
of possible liquid-gas phase transition, play an important role, see Sect. 3. 

Except deuterons, both tritons and ^He are only a small fraction of the total 
multiplicity in heavy ion collisions at beam energy around 100 AMeV. In this 
energy range the reaction dynamics can be described in terms of a relatively 
large compression followed by a faster expansion and vaporization. In such 
a case it is quite suitable to apply a coalescence model for determining the 
cluster production. However a simple coalescence model can only be applied 
if the role of correlations play a negligible role and the binding energy of the 
formed clusters can be discarded. The last issue is expected to affect more 
the absolute value of the yields and can be overtaken looking at the ratios 
of particle yields or spectra. In fact the t/^He ratio which is relevant for the 
isospin physics is expected to be independent on the binding energy effect, 
that is quite similar for the two clusters. 

Within a coalescence model, using the freeze-out distributions obtained in 
IBUU, the behavior of d, t, ^He production has been studied again for the 
^^Ca + ^^Ca system in central collisions, ref.[116]. The feed-down from heavier 
fragments is not taken into account, however especially at large kinetic ener- 
gies, its contribution is less relevant. It is found that the yield of light clusters 
is quite sensitive to the symmetry potential. As it can be inferred from the 
result on the particle emission the stiff symmetry energy leads to a larger yield 
of clusters, in particular the yield of '^He can be a factor two larger than in 
the soft case. The stiff symmetry energy induces a stronger pressure gradient 
causing an earlier emission of neutrons and protons. Furthermore in the stiff 
case neutrons and protons are more correlated in time, see Fig.4-1 and 4-2, 
enhancing the probability of ^He formation in phase space. If one looks at the 
t/^He ratio, the effect of the symmetry energy is reduced because both yields 
are modified so that the increase in the ratio is just about 15%. 

In Ref. [117] the ratio of neutron to proton as a function of their kinetic energies 



53 



was found to be quite sensitive to the symmetry energy. This appears also in 
the hght cluster production as recently found by Chen et al. [118,116]. In 
particular the t/^He ratio shows a quite different qualitative behaviour: while 
it increases with the kinetic energy for the soft symmetry energy, it decreases 
for the stiff one. This is due to the fact that, using a soft symmetry energy, 
a larger number of neutrons is emitted during the early stage of the collision, 
thus enhancing the production of tritons with high kinetic energy. 

Since clusters are expected to be good probes of the high density region reached 
during the reaction dynamics all the previous features will be very sensitive 
to the Momentum Dependence (MD) of the effective forces in the isovector 
channel. For the same ^^Ca + ^^Ca system in ref. [118] it has been shown that 
symmetry energy effects on the t,^ He yields are largely reduced when a Gogny 
effective interaction is used. This is not surprising since the isospin-MD is in 
the m* > m* direction, which implies a quenching of symmetry effects for high 
momentum nucleons. Of course an enhancement is expected in the case of an 
opposite effective nucleon mass splitting. Similar information can be obtained 
from the beam energy dependence of the yields. 



4-3 Collective Flows 

In the search for the EOS of the nuclear matter in various conditions of density 
and temperature a chief role is played by the study of the collective flow, which 
is a motion characterized by space-momentum correlation of dynamical origin. 

The build up of sideward and elliptic flow is realized around the higher density 
stage of the reaction and thus it is a powerful tool for the search of the high den- 
sity behaviour of the symmetry energy. It represents a mean of investigation 
very general, giving information on the dynamical response of excited nuclear 
matter in heavy-ion collisions, from the Fermi energies up to the ultrarelativis- 
tic ones, in the search for a phase transition to QGP, [121,122,123,26,124,125]. 
The collective motion can be characterized in several ways that pin down 
different kind of space-momentum correlation that can be generated by the 
dynamics. The kind of collective flows that have been suggested and employed 
to get information on the equation of state can be divided into three categories: 
radial, sideward and elliptic. 

The sideward and elliptic flow have been and are currently useful tools for 
the study of the compressibility of symmetric nuclear matter. In the search 

for the density behaviour of the symmetry energy similar concepts can be 
exploited but highlightening the difference between neutrons and protons or 
light clusters with different isospin. We will define now the different types of 
collective flow and we will discuss the current status of the effects expected due 
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to different Egymip), and related momentum dependence. We will sec that first 
experimental results with stable beam already show hint of the effect of the 
symmetry energy. Thus future, more exclusive, experiments with radioactive 
beams should be able to set stringent constraints on the density dependence 
of the symmetry energy far from ground state nuclear matter. 



4.3.1 Definitions 

Sideward flow is a deflection of forwards and backwards moving particles, 
within the reaction plane [126]. It is formed because for the compressed and 
excited matter it is easier to get out on one side of the beam axis than on the 
other. The sideward flow is often represented in terms of the average in-plane 
component of the transverse momentum at a given rapidity < Px{y) >'■ 



The particular case in which the slope of the transverse flow is vanishining in 
a region around midrapidity is referred as balance energy. It comes out from 
a balance between the attraction of the mean fleld and the repulsion of the 
two-body collisions. 

For the isospin effect the sum over the particles in Eq.4-2 is separated into 
protons and neutrons. In Ref.[127] also the neutron-proton differential flow 

F^^{y) has been suggested as very useful probe of the isovector part of the 
EOS since it appears rather insensitive to the isoscalar potential and to the in 
medium nuclear cross section and, as we will discuss, it combines the isospin 
distillation effects with the direct dynamical flow effect. The deflnition of the 
differential flow Fpn{y) is 

I Niy) 

Fpn{y) = E Pxi^i (4-3) 



where N{y) is the total number of free nucleons at the rapidity y, is the 
transverse momentum of particle i in the reaction plane, and Tj is +1 and -1 
for protons and neutrons, respectively. 

The flow observables can be seen respectively as the flrst and second coeffi- 
cients from the Fourier expansion of the azimuthal distribution [128]: 

dN 

-TT{y,Pt) = 1 + V^1COS(0) + 2^2005(20) 
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where pt — ^Jp^ + Py is the transverse momentum and y the rapidity along 
beam direction. 

The transverse flow can be also expressed as: 

It provides information on the azimuthal anisotropy of the transverse nucleon 
emission and has been used to study the EOS and cross section sensitivity of 
the balance energy [127]. 

The second coefficient of the expansion defines the elliptic flow V2 that can be 
expressed as 

V2{y.Pt)^{^^) 
Pt 

It measures the competition between in-plane and out-of-plane emissions. The 
sign of V2 indicates the azimuthal anisotropy of emission: particles can be 
preferentially emitted either in the reaction plane (V2 > 0) or out-of-plane 
{squeeze - out, V2 < 0) [128,129]. 

The pt-dependence of V2, which has been recently investigated by various 
groups [124,129,130,131], is very sensitive to the high density behavior of the 
EOS since highly energetic particles {pt > 0.5) originate from the initial 
compressed and out-of-equilibrium phase of the collision, sec e.g. ref.[131]. 
Also at high energy it is allowing to get insight of the partonic stage and 
hadronization mechanism in ultrarelativistic heavy- ion collisions [125]. 



4.3.2 Collective Flows at the Fermi Energies 

The isospin dependence of the transverse collective flow near the balance en- 
ergy was flrst pointed out in Ref.[132], where it is shown that the reactions 
involving neutron-rich nuclei have a significant stronger attractive flow and 
consequently a higher balance energy. A discussion of the effect of a different 
density dependent symmetry energy started only some year later [113,46,25]. 
It was also stressed that the effects are expected to be more evident for 
semi-peripheral collisions (see also Fig. 4-5) around the balance energy. Clear 
signals of the mean field momentum dependence (MD) were also revealed 
([113,133,129]). In the range of the Fermi energies it is indeed essential an 
appropriate treatment of the momentum dependence of the mean field [123]. 

Further studies that look at more n-rich systems investigating the effect of 
different density behaviour of the symmetry energy and based on a more 
exclusive analysis (like pt dependence), have shown that isospin effects on 
collective flows can be found also at higher energies up to beam energies of 
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lAGeV [134,135,65]. 



The search for the density dependence of Egymip) has also to disentangle ef- 
fects coming from the isospin dependence of the nucleon-nucleon cross sections 
[136,137] . In fact little as been done till now to address quantitatevly this prob- 
lem, the best strategy is of course to look for obscrvables that mainly depend 
on the collision rate. A possibility is offered by the pre-equilibrium gamma 
radiation, as mentioned in Ref.[67]. Another possibility is the measurement of 
the radial flow for system with the same mass but different isospin asymmetry. 
Some observables have been found to be weakly sensitive to an overall change 
of the absolute value of the pp, nn, pp cross sections. We will pay particular 
attention at these features. 



4-3.3 Isospin effects around the Balance Energy 

In Ref.[113] the investigation of the isospin effect on transverse flows at Fermi 
energies has been carried out, including also a momentum dependence of the 
mean field. The main structure of this interaction has been already introduced 
in Sect. 2. It is a GBD — MD form in a Skyrme-like effective mean field, 
[66,138]. For a large dynamical range, in particular at the Fermi energies, it 
gives results quite similar to non-local forces of Gogny type [139], that well 
describe collective flows in symmetric systems with realistic compressibilities 
around K ~ 220MeV, [138,140]. For the local part a general Skyrme form is 
used, with two choices for the symmetry term: i) asysoft with the SIII-C (p) 
parametrization and ii) asystiff with the choice C{p) — const — 32MeV (see 
Sect. 2). In this calculation isospin effects on the momentum dependence are 
not included. 

The microscopic transport BNV equations are solved following a test particle 
evolution on a lattice [113,67,141]. Isospin effects on nucleon cross sections and 
Pauli blocking are consistently evaluated. Medium effects on the A^A^ cross 
sections are included following the prescriptions of ref.[136]. A dependence on 
aNN variations is also analysed. 

We present {p, n) transverse flow results for the collisions ^^Fe Fe (n-rich) 
and ^^Ni Ni (n-poor) in the energy range 55 — 105 AMeV where nice 
data arc available for Z = 1,2,3 flows [142,143,144] and some other calcu- 
lations are existing for the proton flows [14,144]. Certainly the asymmetry 
(N — Z)/A of ^^Fe is quite small, nonetheless isospin effects have been found 
both experimentally [142,143,144] and theoretically to be sizeable. In particu- 
lar the comparison of iso-transport simulations with experimental results seem 
to show that a stiff behaviour of Egymip) can give a good description of the 
difference in the proton transverse flows between ^^Fe Fe (n-rich) and 
^^Ni Ni. 
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Fig. 4-4. Energy dependence of flows at hr^ed = 0.45 [145]: (full circles) Fe — Fe 
protons; (open circles) Ni — Ni protons; (squares) Fe — Fe neutrons - (a) asystiff; 
(b) asysoft; (c),(d) same for a^N = ^frn^ no isospin dep. - The full diamonds in (a) 
represent the proton balance energy data of ref . [142] for the Fe — Fe (right) and 
Ni — Ni systems 

In Figs. 4-4, 4-5 we clearly see the isospin effects on the tranverse flows, in par- 
ticular the sensitivity to the density dependence of the symmetry term of the 
nuclear EOS: an as?/siij(f behavior, more attractive for protons above normal 
density for the Fe asymmetric case, gives a clear shift in the balance energy 
as well as a larger (negative) flow at 55 AMeV. i.e. below the balance. Both 
effects are in agreement with the data and are disappearing in the asysoft 
choice. Of course also the isospin and density dependence of the NN cross 
sections is important (see the (c) plots) but we note that a noticeable sensi- 
tivity to the isovcctor part of the EOS is still present. In particular we see 
that the isospin dependence of the mean field is able to keep the transverse 
flow difference between protons in Fe — Fe and Ni — Ni. However a systematic 
study over different systems with more "exotic" isospin content is necessary 
to confirm this result. 



An important effect predicted by the simulations is the clear difference be- 
tween neutron and proton flows. Due to the difficulties in measuring neutrons 
this should be seen in a detailed study of light isobar flows. Moreover we 
like to recall that clusters arc better probing the higher density regions. This 
point is quantitatively shown in the Fig. (4-6) where we present the transverse 
momentum vs. rapidity distributions for ^He — Triton clusters in semicen- 
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Fig. 4-5. Centrality dependence of flows at beam energy 55 AMeV: (stars) Fe — Fe 
protons; (open circles) Ni — Ni protons; (open diamonds) Fe — Fe neutrons; (open 
squares) Ni — Ni neutrons - (a) asystiff; (b) asysoft; (c),(d) same for gnn = 
no isospin dep. - Experimental points: (full circles) Fe — Fe protons; (full squares) 
Ni — Ni protons. 

tral Fe — Fe collisions at 55 AMeV, i.e. below the balance energy [113]. We 
can estimate a 20% larger (negative) flow for the ^He ions, just opposite to 
what expected from Coulomb effects. This appears to be a clear indication of 
the contribution of a much reduced (negative) neutron flow in the case of an 
asystiff force, i.e. a more repulsive symmetry term just above po- The effect 
would disappear in an asysoft choice. Of course even for this observable the 
Isospin — MD will modify the symmetry energy dependence, as already noted 
before for the cluster yields. 



4.3.4 Differential Flows 

It has been suggested that in the study of symmetry energy by mean of the 
neutron and proton collective flow one can even exploit the difference in the 
so called isospin fractionation, i.e. the difference in the neutron-proton pre- 
equilibrium emission. In fact it has been pointed out that a soft/stiff symmetry 
energy leads to a different N/Z emission [112,46,55], see also previous section. 
This can be attributed to the density pressure gradient, which is quite different 
in the two cases. 

We note that the gradient of the pressure should not be confused with the 
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Fig. 4-6. Mean transverse momentum in the reaction plane vs. reduced rapidity for 
light 3He — 3H isobars in the Fe — Fe collisons at 55 AMeV beam energy (i.e. below 
the balance) for semicentral impact parameter, bred = 0.6. Asystiff parametrization. 

so called Kgym, see also the discussion in Sect. 2. Therefore the interpretation 
given in [46] which relies entirely of the different Ksym, around saturation, for 
the two used parametrization (from —69 to +61 MeV) should be taken with 
caution. In fact between the soft and stiff symmetry energy there is, at the 
same time, a change in the slope L, from 54 to 95 MeV, that indeed gives a 
larger contribution to the pressure gradient, according to : 



dPasy ^ Upd{E/A) 



dp 



dp 



sym . sym 



dp" 



— ( g-^' + g^sym 



P (4-4) 



where / = {N — Z)/A is the asymmetry parameter. 

Despite the possibile interpretation, in order to make the analisis of collective 
flow more sensitive to the symmetry potential, the neutron— proton differential 
flow, defined in Eq.4-3, has been introduced [127]. In such a way one combines 
constructively the difference in the neutron-proton collective flow and the 
difference in the number of protons and neutrons emitted. At the same time 
the influences of the isoscalar potential and the in-medium nucleon-nucleon 
cross sections are also reduced. However the measurement of such a differential 
flow demands not only for the measurement of neutron collective flow but also 
for a precise assessment of their number, which most likely is impossible. On 
the other hand the idea to combine more than one isospin contribution in one 
observable is certainly important for elusive effects as those coming from the 
symmetry energy. Moreover the usual problems caused by the neutrons can 
be overtaken looking at clusters. For example one can use the definition of 
differential collective flow and apply it to the He isospin doublet, see 

previous subsection. 
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4-4 Effective Mass Splitting and Collective Flows 

The problem of momentum dependence in the isospin channel is still open. 
Intermediate energies are important in order to have high momentum particles 
and to test regions of high baryon (isoscalar) and isospin (isovector) density 
during the reactions dynamics. Now we present some qualitative features of 
the dynamics in heavy ion collisions related to the splitting of nucleon effective 
masses. 

We report here on expected effects of the isospin MD, studied by means of 
the Boltzmann-Nordheim-Vlasov BNV transport code, Ref. [67], implemented 
with a BGBD — like mean field with a different {n,p) momentum dependence, 
as already introduced in Sect. 2, see ref. [65]. Free (isospin dependent) isotropic 
cross sections are used in the collision term. In order to enhance the role of 
the energy dependence of the mean field we focus on semi-peripheral Au + Au 
collisions at the energy of 250 AMeV and on transverse and elliptic flows of 
neutrons and protons in high rapidity regions, see [146]. 

4-4-i Tranverse flows 

The evaluation of the balance energy for the proton more inclusive Directed 
Flow, integrated over all transverse momenta, in Au + Au collisions represents 
a very good test for the used effective interactions at lower energies. The 




20 60 40 bO 60 
E,„, (MeV/A) 



Fig. 4-7. Balance energy in Au + Au collisions for intermediate impact parameters. 
Solid circles refer to the case m* > m*, empty circles to m* < m*. Triangles 
represent the data from ref. [147] 
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Fig. 4-8. Transverse flow of protons and neutrons for Au+Au reactions at 
250 AMeV, b/bmax = 0.5, in the rapidity interval 0.7 < < 0.9. 

results are shown in Fig.4-7 and compared with NSCL-MSU data. The proton 
directed flows are calculated for semicentral coUisions, i.e. bred = 0.5 [145]. 



Both paramctrizations arc well reproducing the experimental balance energy 
region. As expected we see a slightly earlier balance of the flow in the m* > m* 
case due to some more repulsion for fast protons but the effect is hardly 
appreciable at these energies, see also the previous subsection. We will now 
look more carefully at the mass-splitting effects with increasing beam energy 
and for more exclusive observables. 

n/p transverse flow results for semicentral ^^"^ Au+^^'^ Au coUisions at 250 AMeV 
beam energy, for particles in regions of relatively high rapidities are shown 
in Fig. 4-8 where the left and right panels illustrate the two m* < m* and 
m* > m* cases. It is useful to work with momenta and rapidities normalized to 
projectile ones in the center of mass system (cm), defined as y^^^ = {y/yproj)cm 
and ^ = (pt/Pproj) cm- 
High rapidity selections allow us to find effects of the mean field MD, which 
become more important with increasing momenta. Indeed the different be- 
haviour of mean field depending on mass splitting is evident all over the range 
of transverse momentum shown here. We emphasize that the average mo- 
mentum of these particles is generally beyond the Fermi momentum. Indeed a 
rough estimate in the center of mass system gives already for the beam parallel 
component 

p.^mo- 0.8{yproj)cm ^ 260 - 280 MeV/c 
In this momentum region a smaller effective mass determines greater repulsive 
interaction. When m* < m* neutrons feel greater repulsion than protons and 

their deflection in the reaction plane for pf'^ < 0.6 is close to protons, where 



62 



the Coulomb field is also acting. At variance in the case m* > m* the proton 
flow is greatly enhanced if compared with neutrons. 

We hke to note that in both plots of Fig. 4-8 the proton flow is above the 
neutron one, particularly at low transverse momenta. This is a clear effect of 
the Coulomb repulsion. This is also responsible for the non-symmetric behavior 
of neutron and proton flows in the two cases of opposite mass splitting. 

At high transverse momentum we must take into account the increase of 
the out-of-plane component, so in this case increased repulsion determines 
a greater reduction of transverse flow. It causes a crossing between proton and 
neutron flow in the case m* > m*. This behaviour is clearly seen if we look at 

the different slope of the curves in Fig. 4-8 (right panel) around p^^^ > 0.6. 

A quite suitable way to observe such effective mass splitting effect on the 
tranverse flows is to look directly at the difference between neutron and proton 
flows 

shown in Fig.4-9. 

The variation due to the mass sphtting choice is quite evident in the whole 
Pt range of emitted nucleons. Since the statistics is much smaller at high pts 
(large error bars in the evaluation) the effect should be mostly observed in the 
pf^ < 0.5 range. 




_ f) r\c: I \ \ \ \ I \ \ \ \ I \ \ \ \ 

0.5 1 1.5 
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Pt 

Fig. 4-9. Pt dependence of the difference between neutron and proton transverse 
flows in Au + Au collisions, same energy, centrality and rapidity selection as in 
Fig.4-8. Solid circles refer to the case m* > m*, empty circles to m* < m*. 
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4.4-2 Elliptic flows 



Also in this case a good check of our effective interaction choices is provided by 
some more inclusive data in the medium energy range. In Fig.4-10 we report a 
comparison with the FOPI — GSI data in the interesting beam energy region 
of the change of sign of the elliptic flow for protons, integrated over large 
transverse momentum contributions. The agreement is reasonable for both 
parametrizations. Also in this case we see a slightly larger overall repulsion 
(earher zero-crossing) for the m* > m* case. 

Results on the pt dependence of the neutron/proton elliptic flows in Au + 
Au reactions at 250 AMeV are shown in Fig. 4-11. We can see that the Py 
component rapidly grows at momenta pf'^ > 0.6, causing the flows to fall 
down towards negative values. At low momenta, emitted particles undergo a 
deflection in the reaction plane, more pronounced for nucleons with smaller 
effective mass. In fact, comparing the two panels of Fig.4-11 we see that around 
the maximum in-plane deflection 0.4 < pf^^ < 0.8 neutron and proton flows 
are inverted. We also find different slopes for protons and neutrons around 
Pt^^ > 0.7, depending on the relative sign of mass splitting, in agreement with 
the previous results in the transverse flows. Around pf''^ ~ 1, i.e. projectile 
momentum in the cm, the elliptic flow finally takes negative values, with 
slight differences for the two parametrizations. Again the proton and neutron 
behaviors are not perfectly inverted in the two cases because of the proton 
Coulomb repulsion, in particular at low pt- Also in this case the best observable 
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Fig. 4-10. Energy dependence of the elliptic flow in Au + Au semicentral collisions 



at mid-rapidity, |y' '| < 0.1, integrated over the > 0.8 range. FOPI data with 
similar selections, see ref. [148], are given by the full triangles. Solid circles refer to 
the case m* > m*, empty circles to m* < m*. 
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Fig. 4-11. Elliptic flow of protons and neutrons in Au+Au reaction at 250 AMeV , 
b/bmax = 0.5, in the rapidity interval 0.7 < \y^^^\ < 0.9. 

to look at is the pt dependence of the difference between neutron and proton 
elliptic flows 

Vr''{y,Pt) = Vi{y,pt) - V,^{y,pt) 
shown in the Fig.4-12. 

Again the difference is quite evident in the whole pt range of emitted nucleons. 
We note that also here the statistics is much smaller at high pts (large error 
bars in the evaluation). Anyway the mass splitting effect appears very clearly 
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Fig. 4-12. Pt dependence of the difference between neutron and proton transverse 
flows in Au + Au collisions, same energy, centrality and rapidity selection as in 
Fig.4-11. Solid circles refer to the case m* > m*, empty circles to m* < m*. 
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m the low pr < 0.5 ranee where the results are more reliable. 

As a general comment we stress the importance of more exclusive data, with 
a good pt selection. Indeed from Figs. 4-8, 4-9 (for Vi) and Figs. 4-11, 4-12 
(for V2) we clearly see that more inclusive Vi{y), V2{y) data, integrated over 
all pt's of the emitted nucleons at a given rapidity, will appear to be less 
sensitive to isospin — MD effects. As already discussed from general features 
of the mean field MD we expect a kind of compensation between low and 
high Pt contributions. As already remarked the Coulomb repulsion acting on 
protons in this heavy system tends to reduce the dynamical effects of nucleon 
mass splitting, in particular for low momentum particles. In spite of that the 
signals appear quite clear in the simulations. In any case this represents a 
further indication of the importance of momentum selections. 

4.4-3 Changing the stiffness of the symmetry term 

It is well known that the collective flows for asymmetric systems are sensi- 
tive to the density dependence of the symmetry term of the nuclear EOS, in 
particular at high transverse momentum [113,46,25]. For this reason it is im- 
portant to be sure that signals previously discussed are mainly due to the n/p 
effective mass splitting and not strictly depending on the choice of the stiffness 
of the symmetry energy. In order to check this point we show some calcula- 
tions using a symmetry energy with a density behaviour very different from 
the one used before, but keeping the same n/p-effective mass splittings. We 
refer to it as asy-soft, compared with the asy- stiff pYeviouslj used. For details 
see ref.[65]. Our results show that while the nucleon flows are depending on 
the stiffness of the symmetry term, like in other calculations [113,46,25], the 
transverse momentum behavior of the differences keeps the same sensitivity to 
the effective mass sphttings. Fig. 4-13 illustrates the n/p-flow differences as a 
function of transverse momentum for the same energy, impact parameter and 
high rapidity selections used before. Both transverse Vi~"'{y,pt) and elliptic 
V2~^{y,Pt) flow differences show very small variations from the asy — stiff 
case (i.e. Figs. 4-9 and 4-12) over the whole pf'^ range, in particular in the 
region pf*^ < 0.5, of interest for the better statistics. This result represents 
an important indication that we are selecting experimental observables that 
directly probe the isospin effects on the momentum dependent part of the 
nuclear EOS. 
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Fig. 4-13. pt dependence of the difference between neutron and proton transverse 
(left panel) and elliptic flow (right panel) in Au + Au collisions at 250 AMeV , 
b/bmax = 0.5, in the rapidity interval 0.7 < < 0.9. Calculations are performed 
using the asy — soft symmetry energy. Solid circles refer to the case m* > m*, 
empty circles to m* < m*. 
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5 Isospin Dynamics at the Fermi energies 



5. 1 Stochastic BNV transport theory 

We will show that Heavy Ion Collisions in the Fermi energy domain are rather 
sensitive to the different density /momentum dependences of the symmetry 
term in a quite transparent way. We present studies with stable and unstable 
ions in order to increase the asymmetry range of the nuclear systems. 

Since dynamical instabilities are playing an essential role in the reaction dy- 
namics at Fermi energies it is essential to employ a stochastic transport theory. 
An approach has been adopted based on microscopic transport equations of 

Boltzmann-Nordheim-Vlasov {BNV) type [86,149,150,151,141] where asym- 
metry effects are suitably accounted for [112,113] and the dynamics of fluctu- 
ations is included [52,154]. 

The transport eqiiations, with Pauli blocking consistently evaluated, are inte- 
grated following a test particle evolution on a lattice [141,156,67]. A parametriza- 
tion of free NN cross sections is used, with isospin, energy and angular de- 
pendence. The same symmetry term is utilized even in the initialization, i.e. 
in the ground state construction of two colliding nuclei. 

We consistently include isospin effects in the treatment of the stochastic term 
of the transport equation. Indeed the evolution under the influence of fluctua- 
tions is described by a transport equation with a stochastic fluctuating term, 
the so-called Boltzmann-Langevin equation (BLE) [157,158,159,79,160]. We 
follow two methods to include the fluctuations effects. The Stochastic Mean 
Field {SMF) approach, [52,154,155], and a simplified procedure, computa- 
tionally much easier [83,156,161,162], based on the introduction of density 
fluctuations by a random sampling of the phase space, (Phase Space Sam- 
phng, PSS). The amplitude of the noise is gauged to reproduce the dynamics 
of the most unstable modes [83] . 

The SMF method is based on a fully self-consistent dynamical treatment 
of the fluctuations during the time evolution. We notice that cr^ =< (/ — 
/)^ >, the variance around the mean trajectory of the system, given by the 
Boltzmann-Vlasov evolution in phase space /(r, p, t), in each phase space cell 
obeys the equation of motion [80,154]: 



with 2 D{t) the correlation function of the fluctuating term and the relaxation 



d 
di 



^y+2Dit). 



(5-1) 
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time r(t) = l/{w^ + w~), where and w~ are the transition probabihties 
into and out of the cell. The statistical value cTq = /o(teg)(l — /o(^eg)) at 
equilibrium suggests an ansatz for the correlation function of the fluctuation 
term of the form: 

2D{t) ^ {1 - f)w+ + fw- (5-2) 



i.e. the magnitude of the fluctuations is given by the total number of collisions 
(fluctuation-dissipation theorem) [159]. Then we obtain for the time evolution 
of the difference relative to local statistical fluctuations A = (7^ — /(I — /): 

= -^A. (5-3) 

dt r{t) ^ ' 



Since A = is a solution of Eq.(5-3), if the variance is initially locally set 
equal to its statistical value, it will always follow the evolution of the average 
distribution function according to the relation = /(I — /). Thus from our 
knowledge of the mean "trajectory" /(r, p, t), at any time step we can extract 
local statistical fluctuations in each phase space cell. 

From a projection over coordinate space we derive local density variances wich 
are implemented with a Montecarlo method at suitably chosen time intervals, 
see ref.[155]. In this way we can have also the trajectory branchings in regions 
of instability. The procedure is correct if we assume a local thermal equilibrium 
in each phase space cell. This is appropriate for the problems discussed here, 
namely fragment production in the expansion/separation phase, and it can be 
tested numerically. 

For each system we have checked the equivalence of the SMF approach with 
the sampling PSS method in the description of the collision dynamics, from 
fast particle emissions to the fragment production. The analysis of the results, 
prsented in the following, is based on events collected with both numerical 
procedures. We note that in this way the transport simulation will represent 
a kind of event generator. In general we have accumulated 400 — 500 events 
for each physical case. 

In order to simplify the analysis of the observables most sensitive to isospin 
effects we keep flxed the isoscalar part of the EOS, corresponding to a soft be- 
havior of symmetric nuclear matter (5" A'"M) with compressibility K = 201MeV. 
We will usually consider three choices for the density dependence of the sym- 
metry term. From Fig.2-3 of Sect. 2 these are asy-soft (dashed line), asy-stiff 
(solid line) and asy- super stiff {\oYig-dash.e(i line). In this way we try to disen- 
tangle dynamical symmetry term effects. 
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5.2 Main mechanisms of dissipative collisions: Fusion vs. Deep-Inelastic Com- 
petition 

We start the analysis at relatively low energies. Moreover we will also consider 
medium-light unstable beams eventually available in this energy range fron 
SPIRAL/ 1 — II and other new Radioactive Ion Beam (RIB) facilities. Here 
we investigate isospin effects on the competition between binary events and 
heavy residue formation which occurs at intermediate impact parameter b ~ 

In Fig. 5-1 we show plots, for two choices of the symmetry term, of the time- 
evolution of the density projection on the reaction plane for the mean tra- 
jectory (averaged over several events) in semicentral coUisions {b — Afm) of 
neutron rich {"^^AriN/Z = 1.56) Ni{N/Z = 1.29)) and N ^ Z, neutron 
poor (^^V +^ Ge) ions at ?,QAMeV [112]. 

The effect on the main reaction mechanism of the different density depen- 
dences of the symmetry term is quite evident in Figs. 5-1 (a, 6). It can be un- 
derstood in terms of the amount of repulsion present during the interaction 
of the two surfaces, i.e. in a region below normal density. For neutron rich 
colliding ions the asy-soft choice in general leads to a less repulsive dynamics 
since the symmetry field for protons is more attractive, see Fig. 3 of Sect. 2. 
We sec a stronger interaction between the two partners, a larger dissipation of 
the relative energy, and thus the system will more likely enter a fusion path. 
Fig. 5-1 (a). Experimentally this means a larger incomplete fusion cross section 
at medium energies. 

For neutron poor systems, Figs.5-l(c, d) we see just the opposite. Naively we 
would expect more repulsion since we have larger Coulomb and direct n — p 
collision contributions. Actually in the asy-soft case we have now a dominant 
binary deep-inelastic mechanism. Fig. 5- 1(c). At variance, in the asy-stiff choice 
Fig. 5-1 ((i) we get a clear fusion dominance. This result is amazing for two 
reasons: i) The dynamics of a symmetric N = Z system appears to be isospin- 
dependent; ii) The stiffer symmetry term leads to a larger attraction. A way to 
understand this is that the two exotic neutron poor systems develop a "proton 
skin" which is overlapping in the interaction zone. The protons in the asy-stiff 
case see a stronger repulsive field and arc more easily promptly emitted. As a 
consequence the two partners feel a smaller Coulomb repulsion and more likely 
fusion can be reached. Consistently a larger yield of fast proton emission is 
predicted with the asy-stiff parametrization [112]. We like to recall that some 
sensitivity to isospin effects in deep inelastic collisions was already shown in 
ref. [153] from Landau- Vlasov simulations of Ca-isotope coUisons at IbAMeV . 

These results suggest that in low energy dissipative coUisions an observable 
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t = Ofm/c t = 140fm/c t = 280fm/c t = 420fm/c 



Fig. 5-1. Density plots at different times in a reaction between neutron-ricfi ions 
(a, 6 Ar Ni) and neutron-poor ions {c,d V +^'^ Ge). (a,c): asysoft sym- 
metry term; {b,d): asystiff. See text 



sensitive to the stiffness of the symmetry term can be just the relative yield of 
incomplete fusion vs. deep-inelastic events. Moreover for n-rich systems in the 
asy-soft case we expect more interaction time available for charge equilibra- 
tion. This means that even the binary events will show a sensitivity through 
a larger isospin diffusion. At variance, in the asy-stiff case the two final frag- 
ments will keep more memory of the initial conditions. This point will be 
futher discussed later in the Section. 



The dependence of the interaction time on the stiffness of the symmetry term 
will also influence the competition between binary and neck — fragmentation 
events [163,53], where Intermediate Mass Fragments (IMF, in the range 3 < 
Z < 10) are formed in the overlapping region, roughly at mid-rapidity in 
semicentral reactions. This represents the new leading dissipative mechanism 
at the Fermi energies, as we will discuss in detail later on. Therefore the 
expected sensitivity to the isovector EOS appears very stimulating. 
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5.3 Isospin in Fragmentation Dynamics: Survey of experimental results 

In the last decade the experimental interest on the isospin dynamics in heavy 
ions collisions at Fermi energies and EOS symmetry term effects on various 
reaction mechanisms has been steadly rising, largely stimulated by the per- 
spective of the second generation radioactive beam facilities. 

A large variety of systems and phenomena have been investigated. This clearly 
shows the richness of the field but also reveals some uncertainty in focussing 

on the most relevant physics questions. One of the aims of this report is just 
to contribute on that direction, even trying to select the most sensitive ob- 
servables. In this respect we will first organize a discussion of the performed 
experiments, in a precise correspondence to theoretical analyses and sugges- 
tions presented in the following subsections. 



5.3.1 Isospin equilibration and fragment production 

The systems ^°Ca,^° Ar ,^°CZ-F ^^Fe,^^ Ni have been among the first devoted 
to study the role of the isospin degree of freedom at Fermi energies [164]. 
Let us note that the C/-projectile is a radioactive secondary beam. The total 
mass of composite system (CS) is the same in all cases allowing to explore 
the effects of the N/Z content on the IMF production. At 53AMeV, the 
isotopic ratios of ^^B/^^B, ^BepBe, '^Li/^Li, exhibit an increasing trend with 
the {N/Z)cs- Differences in the yield in central selections, formed with the 
two beam-target combinations, ^^Ar Ni and '^^Ca Fe, show that the 
IMF emission occurs before the isospin degree of freedom is fully equilibrated. 
Similar studies were performed at lower energies, at 33 and AbAMeV , in ref. 
[165], for n-rich vs. n-poor isobaric IMF ratios. Differences in the isotopic 
ratios were seen at AbAMeV but not at SSAMeV. The authors claim that 
with increasing beam energy we see a transition from a regime in which the 
isospin equilibration is reached before fragment production to one when IMF 
emission happens before charge equilibration. This could be a nice evidence 
of isospin transparency for these medium-light systems, but a more accurate 
analysis of the reaction mechanism would be required (centrality cuts) . 

An example of the importance of the centrality cuts can be seen from the 
results obtained at low energy, Ejygam = ^OAMeV, for the heavier A = 58 
systems ^^Fe{n — rich)/^^Ni{n — poor), in various projectile/target combina- 
tions. The isobaric and isotopic ratios ^^B/^^C, ^°Be/^°B, ^Lif^Be, ^^C/ ^^C, 
^^B/^^B, '^Li/^Li, measured at various laboratory angles scale aproximately 
with the {N/Z)cs, [166], clearly showing that both projectile and target nu- 
cleons contribute to the source composition. However, for the ^^Ni H-^* Fe 
and ^^Fe Ni mixed cases, different values of ^'^C/^^C or ^^B/^^C at 40° 
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seems to suggest an incomplete charge equilibration even at this energy for 
the source situated at midrapidity. We will see from the transport simulations 
that at these forward angles in semicentral collisions we can have fragments 
produced in the neck region but still driven from the projectile remnants, so 
the discussion on isospin equilibration is more involved. 

For the same Fe, Ni systems the correlation between the mean number of 
IMFs and charged particle multiplicity, a rough indication of centrality, has 
been investigated as a function of incident energy in [167]. Neutron rich combi- 
nations systematically produce more IMFs for the same value of charged par- 
ticle multiplicy. The effect is robust at the three lower energies, 45, 55, 75AMeV, 
as clearly evidenced when an appropriate normalization to the system total 
charge, accounting for the larger number of protons in the Ni case, is in- 
troduced. However at higher energy, lObAMeV, the isospin dependence of 
multifragment production appears much weakened. This could be an interest- 
ing indication of the onset of different mechanisms for fragment production at 
higher energies, well above the Fermi region discussed in this section. 



5.3.2 Isospin Distillation and Isoscaling 

In other collisions the formed composite systems differ not only in N/Z but 
also in mass. In these cases, for a proper interpretation of isospin effects, 
attention is payed even to the mass dependence of observables. The n-rich/n- 
poor ^'^'^Sn/^^'^Sn systems have been extensively investigated with outstanding 
results. 

The first evidence of the Isospin Distillation (i.e. different asymmetry in the 
liquid and gas phases) in central collisions was in fact observed in ^^^Sn+^'^'^Sn 
and ^^'^Sn+^^'^ Sn reactions at AOAMeV, [168], where the correlations between 
the average number < Nimf > of IMFs and the multiphcities of charged 
particles Nc, of light-charged particles Nlc {Z < 2), and neutrons iV„ were 
measured. It was shown, [168,169], that in the n-rich ^^"^Sn +^'^'^ Sn case the 
multiplicity of the emitted neutrons for the maximum value of < Njmf > 
was much larger than expected from the neutron excess of the 124-1-124 system. 
A clear Neutron Distillation effect was seen for the first time, associated to 
the largest yield of fragment formation: very n-rich gas phase coexisting with 
a more symmetric cluster phase, the so-called healthier fragments of ref.[169]. 

The Sn reactions in all possible n-rich/n-poor target/projectile combinations 
have been extensively studied at the Cyclotron Lab. of MSU at the beam en- 
ergy Eheam = 5(}AMeV , [74,170]. The Isospin Distillation in central collisions 
was observed and analysed using an approach aimed to reduce as much as 
possible the effects of sequential decays by normalizing the isotopic yield from 
different Sn + Sn systems to that of the more charge-symmetric case. It was 
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then realised that the yield ratio of a given isotope in two different reactions 
R2i{N, Z) = Y2{N, Z)/Yi{N, Z) obey an Isoscalinglaw. its logarithm is a lin- 
ear function on the number of neutrons and protons, with the slope given by 
the the isoscaling parameter a and /? respectively, [74,170]. 

The double ratio of the measured isotope and isotone yields can be used, in a 
grand-canonical statistical approach, to extract the isospin composition of the 
free nucleons, i.e. the proton/neutron density ratios.. Observing its dependence 
on the asymmetry of the composite system, {N/Z)cs, it was clearly established 
that a neutron enrichement of the free nucleonic gas takes place. Within the 
same thermodynamical model the absolute value of p„/ Pp in the nucleonic gas 
can be evaluated. The output is in fact significantly larger than the initial 
neutron to proton ratio: between 1.7 to 3.4 for ^^'^Sn +^^'^ Sn (initial 1.24) 
and 2.8 to 8.2 for ^^"^Sn +^^4 (initial 1.48). This effect is certainly present 
but it appears overestimated within the grand canonical picture. In the same 
scheme the isoscaling parameters, being related to the difference of neutron 
and proton chemical potentials, can provide an evaluation of the symmetry 
term of the EOS in dilute matter [170,120]. Very low values are deduced, 
with a large uncertainty due to the determination of the temperature of the 
equilibrated source. Alternative interpretations have been recently proposed, 
e.g. see ref. [171], were the isospin dynamical aspects play a greater role and 
no global equilibration is required. Information even on the density gradient 
of the symmetry energy can be reached, see the following discussions based 
on transport simulations. 

Reactions for similar systems, ^'^'^Sn,^'^^ Xe+ ^'^^Sn,^^^ Sn, have been recently 
studied at lower energies, 28AMeV [172]. Isotopic yield ratios measured at 40° 
laboratory angle show a reduced isospin distillation effect. This could be an 
indication of the quenching at lower energies of the liquid-gas transition for the 
fragment production. However the centrality selection of the events should be 
accurately checked since in semiperipheral collisions other mechanisms, with 
different isospin dynamics, can produce fragments in the same angular range, 
as discussed below. In any case we expect a rise and fall of the liquid-gas mech- 
anism with increasing beam energy for central collisions, and consequently of 
the related isoscaling signal. This trend has been recently observed for the 
58^g 58 reactions systematically studied at 30, 40 and 47AMeV beam en- 
ergies, [75,173,174]. 

Isoscaling behaviour and neutron enrichement of the free nucleon gas have 
been revealed also for ^^'^'^'^'^Sn+ ^^^^^]\fi central collisions at 35AMeV [175]. 
Finally we have to note that the isospin distillation was also noticed in mul- 
tifragmentation reactions induced by high energy protons [176]. 
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5.3.3 Isospin Dynamics for Neck Fragmentation in Semiperipheral Collisions 

It is now quite well established that the largest part of the reaction cross 
section for dissipative collisions at Fermi energies goes through the Neck Frag- 
mentation channel, with IMFs directly produced in the interacting zone in 
semiperipheral collisions on very short time scales. We can expect different 
isospin effects for this new fragment formation mechanism since cluster are 
formed still in a dilute asymmetric matter but always in contact with the re- 
gions of the Projectile-Like and Target-Like remnants almost at normal den- 
sities. 

A first evidence of this new dissipative mechanism was suggested at quite 
low energies, around 19AMeV, in semicentral ^°°Mo +^°° Mo, ^'^^Sn -t--*-^" Sn 
reactions [177,178]. A transition from binary, deep-inelastic, to ternary events 
was observed, with a fragment formed dynamically that influences the fission- 
like decay of the primary projectile- like {PLF) and target-like (TLF) partners. 
Prom in-plane fragment angular distribution a decrease of scission-to-scission 
lifetimes with the mass asymmetry of the PLF /TLF "fission-fragments" from 
3000/m/c down to 200/m/c was deduced. Similar conclusions were reached in 
[179]. Consistent with the dynamical scenario was the anisotropic azimuthal 
distribution of IMF's. In fact the IMF alignement vfith respect to the (PLF*) 
velocity direction has been one clear property of the "neck-fragments" first 
noticed by Montoya et al.[180] for ^^^Xe +^^'^^ Cu at b^AMeV. 

A rise and fall of the neck mechanism for mid-rapidity fragments with the 
centrality, with a maximum for intermediate impact parameters b ~ \hmaxi 
as observed in [181], suggests the special physical conditions required. The 
size of the participant zone is of course important but it also appears that a 
good time matching between the reaction and the neck instabilities time-scales 
is required, as suggested in refs. [163,53]. In fact a simultaneous presence, in 
noncentral collisions, of different IMF production mechanisms at midrapidity 
was inferred in several experiments [182,183,184,185,186,188]. 

An accurate analysis of charge, parallel velocity, and angular distributions has 
been extended to high fragment multiplicities in [186] by Colin et al. They 
have noticed a "hierarchy effect": the ranking in charge induces on average a 
ranking in the velocity component along the beam, fpar, and in the angular 
distribution. This means that the heaviest IMF formed in the mid-rapidity 
region is the fastest and the most forward peaked, consistently with the forma- 
tion and breakup of a neck structure or a strongly deformed quasiprojectile. 
A very precise and stimulating study of the time scales in neck fragmentation 
can be carried out using the new 47r detectors with improved performances on 
mass resolution and thresholds for fragment measurements. Such kind of data 
are now appearing from the CHIMERA collaboration [187,188]. 
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We have reported here the experimental evidences for "neck fragmentation" 
since we can immediately expect large isospin dynamical effects from the pres- 
ence of large density gradients and from the possibility of selecting various 
time-scales for the fragment formation. This will be widely discussed in the 
theory part of this Section. The first evidences of isospin effects in neck frag- 
mentation were suggested by Dempsey et al. in [189] from semiperipheral col- 
lisions of the systems ^'^^'^^'^Xe+ ^^"^'^"^"^Sn at 55AMeV, where correlations be- 
tween the average number of IMF's, Njmf, and neutron and charged particle 
multiphcities were measured. The variation of the relative yields of ^ife/^'^iJe, 
^He/Li with Vpar for several Zpi^p gates shows that the fragments produced 
in the midvelocity region are more neutron rich than are the fragments emit- 
ted by the PLF. Enhanced triton production an midrapidity was considered 
in ref.[181], and more recently in [190], as an indication of a neutron neck 
enrichement. 

P.M. Milazzo et al. [191,192] analyzed the the IMF parallel velocity distri- 
bution for ^^Ni+^^Ni semiperipheral collision at 30AMeV. The two-bumps 
structure for IMFs with 5 < Z < 12, located around the center of mass 
velocity and close to the quasiprojectile {PLF*) source respectively, was ex- 
plained assuming the simultaneous presence of two production mechanisms: 
the statistical disassembly of an equilibrated PLF* and the dynamical frag- 
mentation of the participant region. The separation of the two contributions 
allows for several interesting conclusions. The average elemental event mul- 
tiplicity N[Z) exhibits a different trend for the two processes: in particular, 
the fragments with 5 < Z < 11 are more copiously produced at the midveloc- 
ity region. This experiment has a particular importance since isospin effects 
were clearly observed, in spite of the very low inital asymmetry. The measured 
isotopic content of the fragments is clearly different in the two mechanisms. 
The experimental heavy isotope/light isostope yield ratios, ^''C/^^C, ^"^B/^^B, 
^^BepBe, ^Li/^Li, show a systematic decreasing trend as a function of par- 
allel velocity from cm. to PLF values. All these results indicate a neutron 
enrichement of the neck region, even when initially the system N/Z is close 
to unity. 

Plagnol et al. [193] have examined, for the system Xe -\- Sn between 25 and 
SOAMey, the competition between midvelocity dynamical emission and equi- 
librium evaporation as well as its evolution with incident energy. The onset 
of the neck emission takes place around 2^AMeV and rises with the energy 
while the evaporative part remains quite invariant for a selected centrality. 
Neck matter is found to be more charge asymmetric: more neutron rich iso- 
topes are favored at midvelocity in comparison to evaporation. Evidence of a 
neck-like structure and its neutron enrichement has been seen even in collisions 
with a rather light symmetric target, ^^Ni+ ^^C, '^'^Mg, at M.^AMeV [194]. 
The average N/Z ratio for isotopes with Z = 3,4 exhibits a clear increase 
from the PLF to the midrapidity zone. 
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We close with a couple of comments that indicate the need of new, possibly 
more exclusive, data. The reasons for a preponderance of neutron-rich isotopes 
emitted from the neck region are a matter of debate. Possible explanations be- 
ing, apart the density dependence of symmetry energy, also a fast light cluster 
production, especially of a-particles, which promptly leads to an amplifica- 
tion of neutron excess in the participant matter [195]. For completeness we 
have to mention that different analyses even give conflicting results on the 
n-enrichement of the clusters produced at mid-rapidity [196,197]. This shows 
that the reaction dynamics is in general very complicated, and we can have 
different isospin effects even in competition. The point will be clearified in the 
next subsections. 



5.3.4 Isospin Diffusion 

The isospin equilibration appears of large interest also for more peripheral 
collisions, where we have shorter interaction times, less overlap and a compe- 
tition between binary and neck-fragmentation processes. The specific feature 
at Fermi energies is that the interaction times are close to the specific time 
scales for isospin transport allowing a more detailed investigation of isospin 
diffusion and equilibration in reaction between nuclei with different N/Z asym- 
metries. The low density neck formation and the preequilibrium emission are 
adding essential differences with respect to what is happening in the lower 
energies regime. Tsang et al. [198] probed the isospin diffusion mechanism for 
the systems ^'^'^Sn+ ^^"^Sn at E = bOAMeV in a peripheral impact parame- 
ter range b/bmax > 0.8, observing the isoscaling features of the light isotopes 
Z = 3 — 8 emitted around the projectile rapidity. An incomplete equilibration 
was deduced. The value of the isoscaling parameter a = 0.42 ±0.02 for ^'^'^Sn+ 
^^^Sn differs substantialy from a = 0.16 ± 0.02 for ^^^Sn+ ^^"^Sn. The isospin 
imbalance ratio [199], defined as 

2x — a;124+124 _ ^112+112 
-^»(^) ^ ^124+124 _ ^112+112 (^"'^) 



{i — P,T refers to the projectile/target rapidity measurement, and x is an 
isospin dependent observable, here the isoscaling a parameter) was estimated 
to be around Rp{a) = 0.5 {vs.Rp{a) = 0.0 in full equilibration). The authors 
pointed out that this quantity can be sensitive to the density dependence of 
symmetry energy term since the isospin transfer takes place through the lower 
density neck region. 

The isoscaling analysis of heavy residues identified in peripheral reactions 
^^Kr-\- ^"^"^'^^"^Sn at E = 25AMeV was employed to trace the degree of N/Z 
equilibration as a function of the excitation energy of the primary residue, 
[200]. The evolution of isospin exchange with the centrahty or interaction 
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time, can be established. It was concluded that for the most damped binary 
events, corresponding to the largest excitation energies, charge equilibration 
could be attained. 

5.4 Theoretical interpretation: overview. 

Using the above described Stochastic Mean Field (SMF) approach dissipa- 
tive collisions of Sn, Sn and Sn, Ni isotopes at Fermi energies, in various 
asymmetry combinations, have investigated in great detail . In the transport 
simulations the isovector part of the effective forces is only modified, following 
the general philosophy of the research project which is behind this report. We 
have seen before the amazing richness of the data on isospin effects in reactions, 
still remaining in the stable beam sector. Here we will try to select the key 
observables in order to extract precise information on the isovector part of the 
nuclear EOS in regions far from the normal conditions. In particular we report 
on a study of the 50AMeV coUisions of the systems ^^'^Sn+^^'^Sn ^^^Sn+'^^^Sn 
and ^^^Sn +^^^ Sn, [55], where data are available from NSCL — MSU exps. 
for fragment production and isospin diffusion, [74,170,198]. 

We also analyse the neutron rich, mass asymmetric reaction ^'^'^Sn Ni, at 
SBAMeV bombarding energy, and the corresponding neutron poor reaction 
^^^Sn Ni, to look in detail at the neck fragmentation dynamics and its 
isoscaling properties, [56]. These reactions have been experimentally studied 
with the CHIMERA detector at LNS, Catania [187,188]. As already re- 
marked, the great advantage of the CHIMERA detector for semi central col- 
lision studies is the possibility of measuring target-rapidity data, due to the 
low threshold of the detection telescopes, togheter with a large precision in 
IMF mass separation. As we will see this will allow new correlation measure- 
ments of fundamental importance for understanding the reaction mechanism 
and the related isospin dynamics. 

5.5 From central to semi-central collisions: bulk fragmentation. 
5.5.1 Fragmentation dynamics 

At central impact parameters in this energy range the reaction mechanism 
manifests the features of bulk fragmentation, as obseved in various dynamical 
simulations, [201,54,202]. One can identify quite generally three main stages of 
the collision, as observed also from the density contour plot of a typical event at 
h = 2fm displayed in Fig. 5-2: (1.) In the early compression stage, during the 
first 40 — 50/m/c, the density in the central region can reach values around 
1.2 — 1.3 normal density; (2.) The expansion phase, up to 110 — 120/m/c, 
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brings the system to a low density state. The physical conditions of density 
and temperature reached during this stage correspond to an unstable nuclear 
matter phase; (3.) In the further expansion fragmentation is observed. 

According to stochastic mean field simulations, the fragmentation mechanism 
can be understood in terms of the growth of density fluctuations in the pres- 
ence of instabilities. The volume instabilities have time to develop through 
spinodal decomposition leading to the formation of a liquid phase in the frag- 
ments and a gas of nucleons and light clusters. As seen in the Figure, the frag- 
ment formation process typically takes place up to a freeze-out time (around 
260 — 280/m/c). This time is well defined in the simulations since as the time 
of saturation of the average number of excited primary fragments. The clusters 
are rather far apart with a negligible nuclear interaction left among them. 

An intermediate behaviour, between the fragmentation mechanism illustrated 
above and the formation of an elongated "neck" region between the two col- 
lisional partners, extensively discussed later, is observed at 6 = 4/m (see also 
Fig. 5-8). The "Freeze-Out time", when the nuclear interaction among clusters 
disappears, is decreasing with impact parameter. This gradual transition sug- 
gests that it may be unappropriate to think in terms of a unique fragmentation 
mechanism and more so to assign a fixed size or shape for a multifragmenting 
source, even passing from b — 1 — 2fm to 6 = 4 — 5/m. 

Guided by the density contour plots we can investigate the behaviour of some 
characteristic quantities which give information on the isospin dynamics in 
fragment formation. 

In many of the following figures (i.e. see Fig. 5-3) we report as a function of 
time: 

(a) Mass A in the liquid phase (solid line and dots) and gas phase (solid line 
and dots). 

(b) The asymmetry parameter I — {N — Z) / [N -\- Z) in the gas "central" (solid 
line and squares), gas total (dashed+squarcs), liquid "central" (solid+circles) 
and IMF's (clusters with 3 < Z < 15 , stars). The horizontal line indicates 
the initial average asymmetry. "Central" means a box of dimension 20 fm 
around the CM. of the total system. 

(c) Mean fragment multiplicity Z > 3 whose saturation defines the freeze-out 
time and configuration. 

Moreover we show properties of the "primary" fragments in the Freeze- Out 
Configuration: 

(d) Charge distribution probability P(Z), 

(e) Average asymmetry distribution lav(Z) and 

(f) Fragment multiplicity distribution P(N) (normalized to 1). 
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Fig. 5-2. Sn+ Sn collision at 50AMeV: time evolution of the nucleon density 
projected on the reaction plane. Central b = 2/m collision: approaching, compres- 
sion and separation phases. The times are written on each figure. The iso-density 
lines are plotted every 0.02/m~^ starting from 0.02/m~^. 

5.5.2 Isospin distillation 

We look first at the influence of the mass-to-charge ratio of the colliding ions on 
fragment production comparing, for a fixed asy-EOS (asy-stiff), the neutron- 
rich and neutron-poor Sn + Sn systems. The results are plotted in the Figures 
5-3 and 5-4. For ^'^'^Sn+ ^'^'^Sn we notice a neutron dominated pre-equihbrium 
particle emission during the first 50fm/c. The liquid phase becomes more 
symmetric during the compression and expansion. At variance, the system 
^^"^Sn +^^'^ Sn is on the p-rich side of the valley of stability and along the 
fast nucleon emission the asymmetry of the liquid increases. Hence the pre- 
equilibrium stage appears very important to assess the effective asymmetry of 
the system that will eventually fragment. 

From the beginning of the fragment formation phase of the evolution, between 
110 and 280/m/c, the trends of the liquid and gas phase asymmetry of the 
two systems become similar. Fig. 5-3, 5-4(b). In the "central region" the liquid 

asymmetry decreases while an isospin burst of the gas phase is observed. This 
behaviour is consistent with the kinetic spinodal mechanism in dilute asym- 
metric nuclear matter leading to the Isospin Distillation between the liquid 
and the gas phase. 
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Fig. 5-3. ^^^Sn+^^'^Sn b = 2fm 
collision:time evolution (left) 
and freeze-out properties 

(right), ASY-STIFF EOS. 



Fig. 5-4. ^^^Sn+^^^Sn b = 2fm 
collision:time evolution (left) 
and freeze-out properties (right), 
ASY-STIFF EOS. 



The effects of this process are clearly seen in the IMF isospin content, in 
both cases lower than at the beginning of the spinodal decomposition. Fig. 5- 
3,5-4(c). Opposite trends for fragments with charge above and below Z ^ lb 
can be observed. For heavier products the average asymmetry increases with 
the charge, a Coulomb related effect. However, the asymmetry rises again for 
hghter fragments. This can be a result of the differences in density and isospin 
between the regions in which the fragments grow, due to the fact that not all 
of them form simultaneously, as shown in the density contour plot. 

The charge distribution of primary fragments has a rapidly decreasing trend, 
typical of a multifragmentation process. Fig. 5-3, 5-4(d). The neutron rich 
system favors events with larger IMF's multiphcities, see Fig. 5-3, 5-4(c, f). 
This feature was evidenced also in experiments, as discussed before. 



5.5.3 Isoscaling analysis 

In recent years isospin effects on fragment production have been discussed 
in terms of an isoscaling behaviour, revealed first in multifragmentation data 
[74,170]. It was experimentally observed that comparing two different reac- 
tions, one neutron rich (label 2) and one neutron poor (label 1), the ratio 
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between the yields of a given N^Z isotope Rtl — Y2{N, Z)/Yi{N, Z) follows 
the relation: 

lni?2i = C + Na + Zp (5-5) 



with a, P isoscaling parameters and C a constant. 

Much work has been done on the interpretation of the isoscaling parameters 
in the framework of statistical approaches and the relation to the value of 
the symmetry energy coefficient. For instance, in a grand-canonical statisti- 
cal approach, these parameters are related to the neutron (proton) chemical 
potential differences in the nuclear environments, assumed to be at the same 
temperature, where the fragments are created in the two reactions. It can be 
also related to the ratio of the free neutron (proton) density in the two sources: 

, /?.^ , (5-6) 

In turn, these are related to the symmetry energy properties: 



where esym is the symmetry energy per nucleon (see Sect. 2). We note that 

limi+ - 

dp 



only in the limit ^tmn ~ we get the "liquid-drop" results used in 



refs. [203,204]: 



A//„ = ^esym[{-^) - (^) ] and App = 4esym[{j^) - (^ 



with \App\ > \Apn\ (5-7) 



Since we have always the relation: 



In(^) ^ a-P ^ ^e^yMih-h) (5-8) 



a large interest is rising on the possibility of a direct measurement of the 
symmetry energy in the fragment source from the isoscaling a, (5 parameters 
[203,204] . 

However , it appears that an isoscaling behavior does not necessarily require 
the existence of a fully equibrium, i.e. a grand canonical ensemble formula. 
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Fig. 5-5. Stochastic transport simulations of Sn + Sn central collisions at ^QAMeV. 
Isoscaling behaviors: primary fragments (upper curves); after sequential decays 
(lower curves). 



An isoscaling law is in fact observed also in dynamical simulations of the 
stochastic mean-field type [171]. In Fig. 5-5 we report the isoscaling behavior 

of primary fragments observed in stochastic transport simulations of central 
Sn + Sn collisions at 50 AMeV. In the lower part the isoscaling plots are 
shown after a secondary evaporation calculation. It is seen that the sequential 
neutron emission reduces the a slope considerably. 



Actually one can easily show that isoscaling is expected also in the case where 
the fragment isotopic distributions in the two reactions have a gaussian shape 
around a given average asymmetry. Then the isoscaling parameters are deter- 
mined by the difference of those averages divide by the widhts of the distri- 
butions. The isoscaling parameters are expected to be larger in the dynamical 
relatively to the statistical approach, since the widths of the isotopic distribu- 
tions are smaller than the ones obtained at equilibrium. 
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5.5.4 Symmetry Energy effects 



Let us now investigate how different assumptions on the density dependence 
of the symmetry energy may influence the observables discussed above. In 
Figs. 5-6, 5-7 we report the results obtained respectively with the asy-soft 
and the asy-superstiff symmetry term. We see that the average value of N/Z 
of the fast particles emitted during the expansion phase is affected by the 
symmetry term. In the asy-soft case, below po, neutrons are less bound than 
in the asy-superstiff case (opposite for protons). A more neutron- rich prompt 
particle emission with the soft asy-EOS leads to a lower asymmetry of the 
initial dilute matter undergoing spinodal decomposition. This is seen from 
Figs.5-6(b) (asy-soft) and 5-7(b) (asy-superstiff). As expected the results of 
Fig. 5-3 (b) (asy-stiff) are somewhat in between. 

The influence of the symmetry term is clearly seen in the isospin content of the 
primary intermediate mass fragments. The asy-soft choice is more effective for 
the isospin distillation effect, producing the most symmetric IMF's, as the 
comparison of Figs. 5-6 ,5-7(b),(e) shows. This can be explained energetically 
observing that this choice has a larger symmetry energy at very low densi- 
ties. The differences between the results observed using the two parametriza- 
tions can also be directly interpreted in terms of the corresponding density 
behaviours of the n,p-chemical potentials, see the comments on Fig.2-4 of 
Sect. 2. 

5. 6 From semicentral to peripheral collisions: neck fragmentation 

Summarizing the main experimental observations, we would like to stress the 
following peculiarities of a "dynamical" IMF production mechanism in semi- 
peripheral collisions: 

1. An enhanced emission is localized in the mid-rapidity region, intermediate 
between PLF and TLF sources, especially for IMF^s with charge Z from 3 
to 15 units. 

2. The JAfF's relative velocity distributions with respect to PLF (or TLF) 
cannot be explained in terms of a pure Coulomb repulsion following a statisti- 
cal decay. A high degree of decouphng from the PLF (TLF) is also invoked. 

3. Anisotropic IMF's angular distributions are indicating preferential emis- 
sion directions and an alignment tendency. 

4. For charge asymmetric systems the hght particles and IMF emissions keep 
track of a neutron enrichment process that takes place in the neck region. 

A fully consistent physical picture of the processes that can reproduce observed 
characteristics is still a matter of debate and several physical phenomena can 



84 




Fig. 5-6. ^^^Sn+^^'^Sn b = 2fm 
collision:time evolution (left) 
and freeze-out properties 

(right), ASY-SOFT EOS. 



Fig. 5-7. 124^^^124^^ 5 ^ 2/m 
collision:time evolution (left) 
and freeze-out properties (right), 
ASY-SUPERSTIFF EOS. 



be envisaged, ranging from the formation of a transient neck-like structure 
that would break-up due to Rayleigh instabilities or through a fission-like 
process, to the statistical decay of a hot source, triggered by the proximity with 
PLF and TLF [209,205,206]. Dynamical transport models suggest since long 
time the possibility of observing neck emission [163,207,53,208], [189,181,195]. 
We show in Fig. 5-8 the density contour plots of a neck fragmetation event 
at 6 = 6/m, for the reaction ^'^^Sn+^'^^Sn at 50AMeV, obtained from the 
numerical calculations based on the stochastic transport approach described 
before, [56]. 



5. 6. 1 Neck dynamics and IMF properties: the iso-migration 

In the following we will review results regarding the fragmentation mechanism 
in semiperipheral collisions, as well as the related isospin dynamics and effects 
of the symmetry term density dependence. First we examine again the Sn+Sn 
collisions at 50AMeV, but at impact parameter b — 6/m. We discuss the same 
observables as in central coUisions. 

The development of a neck stucture in the overlap region of the two colliding 
nuclei is evidenced in Fig. 5-8. During the interaction time this zone heats and 
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Fig. 5-8. Sn + Sn collision at 50AMeV: time evolution of the nucleon density 
projected on the reaction plane. Left column: b = Afm. Right column: h = 6/m. 

expands but remains in contact with the denser and colder regions of PLF 
and/or TLF. The surf ace/ volume instabilities of a cyhndrically shaped neck 
region and the fast leading motion of the PLF and TLF will play an important 
role in the fragmentation dynamics. At the freeze-out time, with the neck 
rupture at about 140/m/c, intermediate mass fragments are produced in the 
mid rapidity zone. In some events fragments form very early while, in others, 
they can remain for a longer time attached to the leading PLF's or TLF's. In 
the same figure we note an intermediate behavior between multifragmentation 
and neck fragmentation at b — Afm. 

In the Figs. 5-9, 5-10 we plot the same quantities as before (Sect. 5. 5. 2) for n-rich 
and n-poor Sn + Sn collisions for impact parameter b — 6/m (asy-stiff EOS). 
In the charge distribution probability, Figs. (5-9, 5-10) (d), one can clearly sep- 
arate the PLF /TLF residues and the IMF^s with upper limit charge Z ^ 15. 
We call the IMF's produced by such a mechanism Neck Intermediate Mass 
Fragments (iV/MF's). For neutron-rich systems events with two or even three 
NIMF^s appear, are slightly more hkely. Figs. (5-9, 5-10) (f). The isospin con- 
tent reveals new distinctive features as seen in Figs. (5-9, 5-10) (e). The IMF's 
formed in the neck region are much more neutron rich than the corresponding 
fragments produced in semicentral collisions. Moreover the PLF and TLF 
residues have definitely a lower asymmetry than NIMF^s. To interpret this 
behaviour we have to keep in mind that the neck region is always in contact 
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Fig. 5-9. 124^^+124^^ ^ ^ Qj^ 

collision:time evolution (left) 
and freeze-out properties (right), 
ASY-STIFF EOS. 



Fig. 5-10. ^^"^Sn+^^^Sn b = 6fm 

collision: time evolution (left) 

and freeze-out properties (right), 
ASY-STIFF EOS. 



with a normal-density phase (the spectators) during the fragment formation 
and this determines a different isospin dynamics. Moreover the clusters are 
formed in a nuclear matter not very dilute relative to saturation density. We 
also notice that the neck breaking leads, in some events, to fragments closer to 
one of the spectators with increased interaction between the two. Thus there 
should be a smooth transition to PLF/TLF fast-fission type of events. 

In neck fragmentation we are testing the symmetry energy in different regions 
of nucleon density. It is energetically more favorable to have protons (neutrons) 
migration from (to) the neck region to (from) the more dense spectators, leav- 
ing the nuclear matter in the neck neutron-rich at the time of breaking. Clearly 
the isospin dynamics is ruled by the same energetic arguments as in the case 
of isospin distillation in central collisions. The main difference relative to bulk 
fragmentation is related to the density gradient between the lower density neck 
and the spectators, which can trigger proton and neutron flows in opposite 
directions (isovector mode). Therefore, we propose to call this phenomenon 
isospin migration in contradistinction to the distillation phenomenon driven 
by isoscalar-like unstable fluctuations. 
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Fig. 5-11. ■^'^'^Sn+^'^Ni at 35AMeV. Typical evolution of the density contour plot 
for a neck fragmentation event at 6 = 6/m. 





t=ldo.OJmi'c 
1 1 




1 J 


: <^ 


' v-^ 1 


: \~ 


~ 1 1 


1 


1 1 





Fig. 5-12. Density contour plots as in Fig.5-11. Early stage of fragment formation 
(top) and same events close to the freeze-out (bottom) for four ternary cases a), b), 
c), d) in neck fragmentation. 



5.6.2 Neck dynamics at lower energy 



It is also interesting to study the features of the neck dynamics at lower energy. 
Here we consider the ^'^^Sn+ '^^'^Ni semiperipheral collisions at incident energy 
around 35AMeV, where nice data are appearing from the CHIMERA col- 
laboration at the LNS [187,188]. The reaction evolution at impact parameter 
b — 6/m is shown in Fig.5-11. Calculations are performed using the asy-stiff 
parameterization. The dynamics corresponds to a mainly two-center system. 
We notice the superimposed motion of the PL and TL pre-fragments linked to 
the formation of a neck-like structure with a fast changing geometry between 
40/m/c and 140 - 160/m/c. 

Neck instabilities favour the appearance of NIMF^s, after 150/m/c, in a 
variety of places and ways as can be seen by looking at Fig. 5-12. For four 
events, we illustrate two characteristic stages, the early phase of the fragment 
formation process and the configuration close to freeze-out. 



88 




b(fm) 



Fig. 5-13. Sn + Ni reactions: Impact parameter dependence of the neck fragmenta- 
tion probability (ternary events). White circles: neutron rich reaction. Black circles: 
neutron poor reaction. Asystiff EOS. 



Apart from the events corresponding to the fast NIMF production we have a 
large fraction of binary events, with a presence in the exit channel of only pri- 
mary excited PLF's and TLF's. Induced deformations suggest that they can 
split, also asymmetrically, on time scales long with respect to neck fragmenta- 
tion, but much shorter that in an equilibrium fission process. In Fig. 5-13 we 
plot the dependence on the impact parameter of the ternary event probability. 
The maximum value, about 25%, is attained around mid-centrality, between 
b = 6 — 7/m, decreasing on both sides to about 10% for b = 5fm and b = 8fm. 
It is even smaller at 6 = 4/m, in spite of a stronger dissipation. The longer 
interaction time favors a reabsortion of the neck matter at this energy. This 
is a different behaviour in comparison to the dynamics at bOAMeV, where, 
as we have already emphasized, the higher available energy makes possible 
the transition to multifragmentation with increasing centrality. At greater im- 
pact parameters, b = 9fm, a smaller overlap and a faster separation are also 
suppressing this mechanism. 

Prom the simulations we can extract an interesting information on the time 

scale of the Neck-IMF production. In Fig. 5-14 we show, for different impact 
parameters, the probability distribution of the time interval between the in- 
stant of the first separation of the dinuclear system and the moment when a 
Neck-7MF is identified (scission-to-scission time). A large part of the NIMFs 
are formed in short time intervals, within 50/m/c. 

Finally we would like to remark that the neck fragmentation shows a depen- 
dence on the nucleon-nucleon cross sections and the i?0 S'-compressibility. The 
latter point is particularly interesting since it seems to indicate the relevance 
of volume instabilities even for the dynamics of neck. This appears consistent 
with the short time scales shown before, see also the discussion in ref . [56] . 



89 



0.4 

0.3 

a. 0.2 

0.1 

0.0 



In 


b=5tm 


-, , , , 1 , , , , 1 , , , , 1 





25 50 75 100 125 150 



0.5 
0.4 

0.3 t- 
0.2 
0.1 
0.0 



b=6fm 




25 50 75 100 125 150 



0.5 






0.4 






b=7fm 


0.3 








0.2 










0.1 






0.0 




,i,,,,i:Tm,,,^, 



25 50 75 100 125 150 

t(fm/c) 



0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



b=8fm 




25 50 75 100 125 150 

t(fm/c) 



Fig. 5-14. The probability distribution of scission-to-scission time in the neck frag- 
mentation for impact parameters from 5 to 8fm. Asystiff EOS. 

5.6.3 Analysis of kinematical observables 

The nonstatistical features of the NIMF production are revealed in various 
kinematic correlations. The corresponding observables can be measured in 
exclusive experiments. 

An interesting kinematical observable is the asymptotic relative velocitity of 
the neck-produced IMFs with respect to the PLF (TLF), Vrei{PLF, TLF) = 
|vpLF,TLF — vimfI- This is compared with the relative velocity reached in a 
pure Coulomb-driven separation, signature of a statistical fission process of 
a compound PLF* or TLF* system, as provided by the Viola systematics 
[210,211]: 



\ 



-(0.755 



Z1Z2 



red 



7.3) 



(5-9) 



where Ai,A2,Zi, Z2 are the mass and charge numbers of the fission products 
and Mred is the corresponding reduced mass. 

For each Neck-/MF we can evaluate the ratios r = Vrei{PLF) / Vyioia{PLF) , 
(rl = Vrei{TLF)/vvioia{TLF)). In Figure 5-15 we plot rl against r for each 
NIMF. We call such a representation a Wilczynski — 2 plot [212]. The solid 
lines represent the loci of the PL-{r = 1) and TL-{rl = 1) fission events 
respectively. The values (r, rl) appear simultaneously larger than 1 suggesting 
a weak NIMF correlation with both PLF and TLF, in contrast to a statistical 
fission mechanism. The process has some similarities with the participant- 
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Fig. 5-15. Wilczynski-2 Plot: correlation between deviations from Viola systematics, 
see text. Results are shown for two asy — EOS. 

spectator scenario. However the dynamics appear much richer than in the 
simple sudden abrasion model, where the locus of the r — rl correlation should 
be on the bisectrix, apart the Goldhaber widths, see ref.[205]. Here the wide 
distributions of Fig.5-15 reveal a broad range of fragment velocities, typical of 
the instability evolution in the neck region that will lead to large dynamical 
fluctuations on NIMF properties. 

Another very selective correlation is provided by the anisotropy in the frag- 
ment emission, measured by the in-plane azimuthal angle, ^piane, defined as the 
angle between the projection of the PL — IMF scission axis onto the reaction 
plane and the separation axis between PLF and TLF, [178]. A |$p;ane| — 
collects events corresponding to asymmetric "fissions" of the PL-system very 
aligned along the outgoing PL — TL separation axis. At variance a statisti- 
cal fission dominance would correspond to a fiat ^piane behavior. The ^piane 
distributions, corresponding to the same NIMF events analysed before, are 
covering a quite limited angular window, close to the full alignment configura- 
tion, ^piane = 0°. The distribution becomes wider when we approach the two 
r, rl = 1 fines of the Fig.5-15. 

We can transform in time — scales the correlations discussed before. In fact an 
induced asymmetric, fast fission can also manifest some deviations from Viola 
systematics. From our simulations, this mechanism takes place on longer time 
scales compared to neck fragmentation. If a light fragment escapes later from 
its PL/TL partner, friction certainly will attenuate the dynamical effects and 
relative velocities will deviate less from Coulomb values. In the plane r — rl 
such events are located closer to the line r = lorrl = l.A prolonged contact 
with the PLF or TLF will induce also different angular distributions in com- 
parison to the neck fragments (less alignement, wider ^ plane distribution). As 
a limiting case we recover the statistical fission mechanism for a IMF pro- 
duction from an equilibrated PL/TL excited residue, with no Viola-deviations 
and no alignement. We conclude that an analysis based on the Wilczynski-2 
plot, Fig.5-15, of the IMF produced in semi central collisions will be able to 
select and study fragments formed on a large variety of interaction time-scales. 
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Fig. 5-16. NIMF isospin content for asysoft (circles), a.sy stiff (rombs) and 
superasystif f (squares) EOS as a function of r-deviation from Viola systemat- 
ics, at impact parameters from b = 5fm to 8/m. The two solid lines represent the 
mean asymmetries of the projectile (top) and target (bottom). 



of particular interest for the mass and isospin dynamics. 



5.6.4 Symmetry term effects 

Prom the simulations of ref.[56] for the ^^'^Sn+^^Ni system it appears that the 
symmetry term of the EOS does not seem to influence sensitively the main 
features of the neck mechanism, i.e. probability of ternary events, deviations 
from Viola systematics and NIMF angular distributions. However, one has 
to consider that with the studied system we cannot reach high charge asym- 
metries, the asymmetry parameter I = (N — Z)/A ranging from 0.193 for the 
projectile to 0.125 for the target, with an average / = 0.17. 

A more promising observable seems to be the isotopic content of the Neck- 
IMF. For the three asy — EOSs introduced in Sect. 2 we plot in Fig. 5-16 
the average isotopic composition / of the NIMF^s as a function of the PLF 
r-deviation from Viola systematics. At all impact parameters clear differences 
are evident. The average asymmetry does not depend strongly on r. We no- 
tice however that it increases with the stiffness of the symmetry potential 
around and below saturation. The superasystiff parameinzation, i.e. with an 
almost parabolic increasing behavior around po, produces systematically more 
neutron rich NIMFs. This effect is clearly due to a different neutron/proton 
migration at the interface between PL/TL "spectator" zone of normal density 
and the dilute neck region where the NIMFs are formed, as it was already 
observed for the reactions at 50 AMeV. 
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Fig. 5-17. Isoscaling in neck fragmenta- Fig. 5-18. Isoscaling in neck fragmenta- 
tion: lnR2i dependence on N. tion: lnR2i dependence on Z. 

From the Fig. 5- 16 we see that the isospin content of NIMFs carries important 
information on the isovector part of the effective interaction and the related 
isospin dynamics in the early stages of the reaction. We expect a different 
isospin pattern for the IMF's produced at later times from induced and/or 
statistical fission of the more charge symmetric PLF*'s and TLF*'s. The 
Wilczynski-2 plot will help to make the selections. We have to remind that 
all the results presented here refer to properties of primary (excited) frag- 
ments. The neutron excess signal appears likely to be washed out from later 
evaporation decays. A reconstruction of the primary fragments with neutron 
coincidence measurements would be very important. 



5.6.5 Isoscaling analysis 

Relative to the earlier discussion it is of interest to test whether the isoscaling 
behavior can manifest itself even in the neck-fragmentation, which is a clear 
dynamical process, related to short characteristic time scales. 

In Figs. 5-17, (5-18) is shown the (Z) dependence of /ni?2i, for Z = 1 
to Z = 9 light fragments produced in the neck region, as obtained in the 
simulations with an asystiff — EOS parametrization. R^i is the yield ratio 
of the n-rich ^'^'^Sn Ni vs. the n-poor ^^'^Sn Ni system, see before 
Sect.5.5.3. 

We do see the isoscaling signal, an exponential N- and Z- dependence of 
the yield ratios with quite well defined a, (5 slopes. Although we cannot use 
explicit equilibrium relations, like Eq.5-8, we still expect a symmetry energy 
dependence of the isoscaling parameters. Since the fragment formation takes 
place in the neck region, we predict that its isospin content will dictate the 
values of the isoscaling parameters. Indeed, we may assume that for a neutron 
poor system, closer to the symmetric case, the differences between the various 
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asysoft 


asystiff 


superasystif f 


a 


0.69 


0.95 


1.05 




-0.67 


-1.07 


-1.18 



Table 6 

The isoscaling parameters a and /3 in neck fragmentation for three asy — EOS. 

asy — EOS on the isotopic and isotonic distributions are reduced, and in a 
first approximation, identical. At variance, for the neutron rich system, we 
have seen that passing from asysoft- to superasystif f -EOS more neutron- 
rich Neck-7MFs are formed. Therefore, the corresponding distributions have 
to be steeper. As reported in Table 6, a nice increase (in the modulus) of the 
isoscaling parameters with the increasing stiffness of the symmetry energy is 
observed. 



5. 7 Charge Equilibration in Peripheral Collisions 

5.7.1 Isospin diffusion in presence of inhomogenous density distributions 

We turn now to discuss the isospin equilibration mechanism in even more 
peripheral collisions where only two primary fragments {PLF — TLF) are 
observed in the exit channel. When the ions in the entrance channel have 
a different N/Z ratio, we expect that isospin diffusion will lead the system 
towards a more uniform asymmetry distribution. The degree of equilibration, 
correlated to the interaction time, may provide insights on transport properties 
of fermionic systems [213,214], in particular on the diffusion coefficient of 
asymmetric nuclear matter, [215,216]. 

Here we are focus on a charge asymmetric collision ^'^'^Sn-\-^^'^ Sn, at 50AMeV 
bombarding energy, to which we refer as the mixed system, (M), where some 
data also exist, [198]. The simulations are performed for peripheral collisions 
at impact parameters 6 = 8,9, 10/m. Binary events are selected. We define the 
average interaction time, tc, as the time elapsed between the initial touching 
and the moment when the PLF and TLF reseparate. From our simulations 
we obtain tc ~ 120, 100, 80/m/c respectively for the three impact parame- 
ters. Typical density contour plots, at 6 = 8fm and b = 10/m, are shown 
in Fig. 5-19. In the overlap region, after around 40/m/c, the formation of a 
lower density interface is evidenced. The isospin migration takes place during 
the transient contact between the two spectator regions with density close 
to the normal one, separated by a dilute neck region. Thus we have con- 
centration and density gradients ruling the isospin diffusion. In contrast, in 
deep-inelastic collisions at lower energies, the isospin equilibration is driven 
by the N/Z difference between two interacting nuclei with an uniform density 
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Fig. 5-19. ^^^Sn+^^^Sn collision at 50 AMeV, b = 8fm (up) and b = 10 fm (down): 
density contour plots. 

profile until the separation, i.e. only the concentration gradient is active. 

In binary events the charge asymmetry of primary projectile (target)- like 
fragments, PLF and TLF provide the essential information about the isospin 
equihbration rate. At 6 = 8/m, however, around 25% of the events were 
ternary and as shown in the previous section for even more central events this 
mechanism becomes dominant. The IMF's formation in the overlap region 
will influence the final isospin distribution, rending more difficult the inter- 
pretation of the results. Thus we select here only binary events. We quantify 
the degree of equilibration through the isospin imbalance ratio Ri Eq.(5-4), 
[198,199]. The measured isospin dependent quantities are directly the isospin 
content of the fragments at separation (in the i — P,T rapidity regions) , 
as they result from the mixed reaction 124/112, from the reactions between 
neutron rich nuclei, (124/124), and between neutron poor nuclei, (112/112), 
respectively. We report its dependence on tc for asysoft (squares) and asy- 
superstiff (circles) EOS in Fig. 5-20. A good degree of isospin equilibration 
corresponds to a "convergence" to of both Rp (upper curves) and Rt (lower 
curves). Prom our results we conclude that an asystiff-like EOS provides a 
better agreement to the experimental observations, shown as arrows in Fig.5- 
20, [198,216]. In fact in the MSU experiment there is no particular selection 
on binary events. We expect that the presence of events with production of 
Neck-7MFs will induce an apparent larger isospin equilibration since more 
neutrons will migrate to the neck region (and viceversa for protons), as dis- 
cussed before. 

A clear difference between the two equations of state is evident especially for 
the longer interacting time, h — Sfm case. Smaller values of the isospin imbal- 
ance ratios for asysoft EOS point towards a faster equilibration rate. In Refs. 
[198,216] a possible explanation was proposed related to the observation that 
below normal density the asysoft EOS has a larger value of symmetry energy. 
Therefore an enhanced isospin equilibration has to be triggered if the diffusion 
takes place at lower density. In fact the mechanism of charge equilibration is 
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Fig. 5-20. ^^^Sn +^^^ Sn 50 AMeV collision: interaction time evolution of the pro- 
jectile (upper) and target (lower) isospin imbalance ratios. Squares: Asysoft EOS. 
Circles: Superasystiff. The arrows correspond to the data of Ref.[198]. 

more complicated at these energies due to reaction dynamics (fast particle 
emissions, density gradients, etc.), with interesting compensation effects, as 
shown in the following. 



5.7.2 Isospin diffusion and EOS dependence 

Let us focus on the "mixed" ^"^"^Sn +^^^ Sn case. The isospin content of the 
two residues at separation is determined by the interplay between the nucleon 
emission from each ion and by the nucleon transfer through the neck: 



Ip = -^{hp - -J-igP - -T—ipT + -7— hp) (5-10) 
^OP ^OP ^OP 

It — -j— (-/or - -j — igT + -j — iPT — Itp) (o-J-J-j 

At Aqt ^ot ^ot 

Here /j, A represent asymmetries and masses for: i — P,T, PLF,TLF at 
separation; i = OP, OT, initial Projectile/Target; i = gP,gT, fast emitted 
nucleons from Projectile/Target (the label g stands for "gas"); i = PT,TP, 
nucleons transferred from projectile (target) to target (projectile). 

In Fig. 5-21 we plot the time evolution of the quantities Igp, Igx and Agp, Agx 
for asysoft and superasystiff EOS. We remark that Igp it is much larger than 
Iqp ^ 0.19. The same is true for the target but the difference is smaller. Con- 
sequently the pre-equilibrium emission reduces the N/Z difference between 
the two nuclei thus competing with the transfer process. 

When the two equations of state are compared a more neutron rich composi- 
tion of pre-equilibrium is seen in the asysoft case since below normal density, 
from where most of the emitted nucleons originate, the neutrons (protons) are 
less (more) bound for this asy-EOS. The differences between the two asy-EOS 
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Fig. 5-21. ^^^Sra +^^'^ Sn b = 8fm (left) and b = 10/m (right) collision: time 
evolution of isospin content and mass of projectile (target) pre-equilibrium emitted 
particles Igp, Agp {IgT, -^qt)- The squares (circles) indicate the asymmetry directly 
transferred from projectile (target) to target (projectile) IpT (Itp), at separation 
time. Asysoft: full symbols, solid lines. Superasystiff: empty symbols, dashed lines. 
Horizontal lines give the initial target and projectile asymmetries. 

are reduced at larger impact parameters, as seen in the results for b = 10/m 
in Fig. 5-21. The neutron rich nuclear skin seems to play a more important 
role for such more peripheral collisions. For the projectile both preequilibrium 
emission and nucleon transfer drive the system towards a more symmetric 
configuration. The two processes tend to compensate for the target. Therefore 
the projectile asymmetry has a more pronounced deviation from the corre- 
sponding initial value in comparison to the target. 

From the Fig.5-21 we realize that asysoft interactions are very efficient in trans- 
ferring isospin to the "gas" while the isospin diffusion in the nuclear medium 
appears more reduced. The opposite is happening with the stiff symmetry 
term. The final result on isospin equilibration comes out from a kind of com- 
pensation, with the asysoft case showing a larger equilibration, see Fig.5-20. 
The evidence that the isospin content of the free nucleon emission is largely 
driving the isospin equilibration is extremely important. We can then look for 
a consistent picture of several obscrvables ranging from isospin effects on fast 
nucleon emissions and collective flows, see previous Sect. 4, to the imbalance 
ratios of N/Z or other isospin dependent quantities. 

Before closing this discussion we would like to elaborate a little more on the 
physics of isospin diffusion in the nuclear medium at the Fermi energies, and in 
particular on the conclusion from the simulations of a much stronger effective- 
ness of a stiff symmetry term. From the Fig.5-21 we see a clear dependence 
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on the cLsy-EOS of the asymmetry transferred between the two nuclei. In 
particular, for b — 8fm we observe that 



jisupasystif f) ^ j{asysoft) ^ = 192 

j^supasystiff) ^ ^ q^q^ ^ /(f^^^/*) (5-12) 

The origin of these inequalities lies in the presence of density gradients due 
to the more dilute interface. For the superasystiff EOS the neutron migration 
towards the neck region is favored from both participants. This explains simul- 
taneously 7(f«f«-f-*^//) > i^p and 7^-«f«-f-*^//) > i^^. For the asysoft EOS this 
effect is weakened because the symmetry part of the mean- field for protons 
and neutrons does not change significantly with density below saturation. 

This argument can be made more explicit if we observe that the proton and 
neutron migration is dictated by the spatial gradients of the corresponding 
chemical potentials ^„(pp,p„,T) and fip{pp, pn,T) [217]. The currents of the 
two species can be expressed as follows: 



jn = -CtVHn{Pp,Pn,T) 
jp = -CWpp{pp, Pn, T) 

where C is a constant. 
In terms of Landau parameters, i.e.: 

= V + ^o'', <l-n,p q'^n,p (5-13) 

apg> 



-Ct[\ 



-Ct[\ 



\dpn, 

'dpp^ 

dpnj 



Pp,T 



'dpn\ 

' dpp \ 



we obtain: 



^^2 7V„ 2 TVj ^ 2^ 7V„ ivj ^ 

= -DinVp - D^nVI (5-14) 

■'^ 2 Np 2 Np ^ ^ 2^ Np Np ^ ' 

= -D^pVp - D2pVI (5-15) 

In the last expressions we have changed the variables Pn, Pp '^o p = Pn + Pp and 

I ^{pn- Pp)/p- 
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Hence it is clear that not only concentration gradients arc important in the dy- 
namics of isospin migration and equilibration, but also density gradients play 
an essential role. Actually the final result comes out from a delicate balance 
between the two effects. The stronger density dependence of the stiff sym- 
metry term below normal density enhances the density gradient contribution, 
finally leading to the larger isospin diffusion of Eqs.(5-12). As already shown 
the influence of the symmetry term on fast nucleon emissions can partially 
compensate the effect. 
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6 Effective interactions in the isovector channel: relativistic ap- 
proach and the role of the 5 meson 

The QHD effective field model represents a very successful attempt to de- 
scribe, in a fully relativistic picture, equilibrium and dynamical properties of 
nuclear systems at the hadronic level [27,28,29]. Consistent results have been 
obtained for the nuclear structure of finite nuclei [218,219,220], for the NM 
Equation of State and liquid- gas phase transitions [19] and for the dynam- 
ics of nuclear collisions [146,221]. Relativistic Random-Phase- Approximation 
(RRPA) theories have been developed to study the nuclear collective response 
[222,223,77,224,225,226]. 

This report focusses on the dynamical response and static (equilibrium) prop- 
erties of Asymmetric Nuclear Matter (ANM). We use a relativistic kinetic 
theory with the aim of a transparent connection between the collective and 
reaction dynamics and the coupling to various channels of the nucleon-nucleon 
interaction. One of the main points of discussion is the relevance of the cou- 
phng to a scalar isovector channel, the effective 5[ao(980)] meson. In fact the 
scalar-isovector coupling is not necessarily connected to the exchange of a real 
5-meson, but such a channel appears automatically in relativistic Hartree-Fock 
and Dirac-Brueckner-Hartree-Fock treatments from exchange and correlation 
contributions [33,36]. 

The introduction of the isovector-scalar channel in covariant approaches should 
play a key role in the effective interaction in asymmetric matter, see ref.[12,13]. 
This point has not received great attention before for two main reasons: 

i) The (5— channel has not been considered a priori, just on the basis of the 
weak contribution to the free Nucleon-Nucleon interaction, [227,228]. But in 
the spirit of the Effective Field Theory as a relativistic Density Functional 
Theory, (the EFT/DFT framework [229]), the relevance of this channel could 
be completely different in nuclear matter, due to medium and many-body 
effects, as noted before. In particular, we expect a large contribution from 
exchange terms of the strongly coupled isoscalar channels, [36,12,13]. 

ii) The extension is not well supported by the existing data on exotic nu- 
clei, as remarked in the refs. [34,220,230,231]. Clearly these conclusions are 
mainly derived from the lack of information on observables more sensitive to 
the density dependence of the symmetry term. In particular finite nuclei stud- 
ies cannot easily disentangle the effects originating from different isovector 
mesons, since the densities are mainly below saturation, see the discussion in 
refs. [220,34,232]. 

We note that recently , see the conclusions of refs. [30,31], the 5-channel has 
been reconsidered as an interesting improvement of covariant approaches in 
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the framework of the EFT/DFT philosophy. One of the main tasks of oTir 
work is then to try to select the dynamical observables more sensitive to it, 
[12,13]. 

In this respect the results reported here on the collective response and reaction 
dynamics can be useful in order to solve the open problem of the determina- 
tion of the scalar-isovector coupling. Since contributions to this channel are 
mainly coming from correlation effects [36], the correct microscopic approach 
should be to derive the relative coupling constant, in a Q-f^-D-Relativistic Mean 
Field (RMF) framework, from Dirac-Brueckner-Hartree-Fock (DBHF) cal- 
culations. Several attempts have been recently performed, see [33,34,233,234], 
but the results are up to now not fully model independent. 

An important outcome of our work is to show that the two effective couplings, 
vector and scalar, in the isovector channel influence in a different way the static 

(symmetry energy) and dynamic (collective response, reaction observables) 
properties of asymmetric nuclear matter. This will open new possibilities for 
a phenomenological determination of these fundamental quantities. 

We will often derive transparent analytical results. In order to show also some 
quantitative effects of the dynamical contribution of the 5-channel we have to 
fix in some way the corresponding coupling. Wc have used a constant value (sec 
Table 7) extracted from the DBHF analysis of refs. [33,34], where it actually 
appears not strongly density dependent in a wide range of baryon densities. 
Some results are also presented with Fock correlations explicitly accounted for 
{NLHF case, see following). 



6. 1 QHD effective field theory 

We start from the QHD effective field picture of the hadronic phase of nuclear 
matter [27,28,29]. To include the main dynamical degrees of freedom of the 
system we will consider the nucleons coupled to the isoscalar scalar a and 
vector cu mesons and to the isovector scalar 5 and vector p mesons. 

The Lagrangian density for this model, including non-linear isoscalar/scalar 
cr-terms [235], is given by: 

L = 'ipi^^iid'^ - g^V - gpB'' • t) - {M - g,<t> - ggr ■ 5 )]V + 

\{d,5- d^5- m's5') - Ig.^G^'' + \mlB,W (6-1) 
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where W^^ix) = d^'V'ix) - ^''V^'{x) and G^'''{x) = a^B'^(x) - d'^B^'ix) 



Here iIj{x) is the nucleon fcrmionic field, (f){x) and Vlx) represent neutral 
scalar and vector boson fields, respectively. 6{x) and B^(x) are the charged 
scalar and vector fields and r denotes the isospin matrices . 

From the Lagrangian, Eq.(6-1), with the Euler procedure a set of coupled 
equations of motion for the meson and nucleon fields can be derived. The 
basic approximation in nuclear matter applications consists in neglecting all 

the terms containing derivatives of the meson fields with respect to the mass 
contributions. Then the meson fields arc simply connected to the operators of 
the nucleon scalar and current densities by the following equations: 

%/U + A%^ + B%^ ^i){x)'il^{x) =p~s (6-2) 

W{x) = g^V^ = Ui^ixh^ipix) = U% , 
B^{x) = g,B^ = f,i;{x)Ym^) , 

6{x) = gsS = fsi^{x)Tip{x) (6-3) 

where $ = g^(j), = {ga/niaf, A = a/gl, B = b/gj, ^ {gu>/m^f, 
fp = {9p/mpf, fs = {gs/msf. 

For the nucleon fields we get a Dirac-like equation. Indeed after substituting 
Eqs. (6-2,6-3) for the meson field operators, we obtain an equation which con- 
tains only nucleon field operators. The equations can be consistently solved 
in a Mean Field Approximation (RMF), where most applications have been 
performed, i.e. in a self-consistent Hartree scheme [28,29]. 

The inclusion of Fock terms is conceptually important [236,237] since it auto- 
matically leads to contributions to various meson exchange channels, also in 
absence of explicit coupling terms. A thorough study of the Fock contributions 
in a QHD approach with non-linear self-interacting terms has been recently 
performed [238], in particular for asymmetric matter [36]. 

6.1.1 Relativistic transport equations with Fock terms 

We now discuss a kinetic approach consistent with the previous approximation. 
We are concerned with a semiclassical description of nuclear dynamics, so that 
the nuclear medium is supposed to be in states for which the nucleon scalar and 
current densities are smooth functions of the space-time coordinates. Within 
the QHD model we focus our analysis on a description of the many-body 
nuclear system in terms of one-body dynamics. Correlation effects can be 
effectively included at the level of coupling constants, as noted in the discussion 
of the results. 
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We move to the quantum phase-space introducing the Wigner transform of 
the one-body density matrix for the fermion field [239,240]. The one-particle 
Wigner function is defined as: 

[F{^,P)U = / d'Re-^'-''{: + f )^«(^ " f ) , 

where a and /3 are double indices for spin and isospin. The brackets denote 
statistical averaging and the colons denote normal ordering. The equation of 
motion is derived from the Dirac field equation by using standard procedures 
(see e.g. [239,240]), as: 



'-d^b^F{x,p)U+p^[-rF{x,p)U - MF^p{x,p) 



9u 



R 

^ J e-^-^ <: M^M,M^-)r-B,{x.) :> 
9l_ [ <. ^^^x+)ilj^{x.)5{x_) :> 

= (6-4) 



with x+ = X + ^ and x_ = x — ^. When we insert the Eqs. (6-2,6-3) for the 
meson field operators we clearly see the appearance of Fock contributions (at 
the lowest order in density matrices). 

The Wigner function is a matrix in spin and isospin space; in the case of asym- 
metric NM it is useful to decompose it into neutron and proton components. 
Following the treatment of the Fock terms in non-linear QHD introduced in 
Refs. [238,36], we obtain the kinetic equation: 



2 
2 



d^j''F^'\x,p)+Yp*iF^'\x,p) - M*F^\x,p) + 



L3,{x)r ± /pJ3m(^)7^ - Lps{x) T hPs^{x)\ F^^\x,p) = 0, (6-5) 



where i — n,p, and here and in the following upper and lower signs correspond 

to protons and neutrons, respectively. Here A = - dp, with acting only on 
the first term of the products and pss = psp ~ Psn and j3^{x) = j^{x) — i^'(a;) 
are the isovector scalar density and the isovector baryon current, respectively. 
We have defined the kinetic momenta and effective masses as: 
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M*^M- Ups{x) ± {x) , (6-6) 
with the effective couphng functions given by: 



~^_ld^_J_^f^xll 3 _3 
PS 8 dpsix) 2p/ ^"^^ dpUx) ^ + 2^^ 8^^ ' 



lrf$(x)_ 5 3 3 
8dps{x)^ A^'^^ A^"^ 8- 
~ __ld^x)_ 1 _ 1 9 

8dps{x)^2^^ 2^p^^^'^ 
1 d<^{x) 1 3 

(6-7) 



where 8 TrF'^{x) = p| + J^j'' + pis + j^.^]^^- We recall that wc are dealing 
with a transport equation so the currents and densities, in general, are varying 
functions of the space-time, at variance with the case of nuclear matter at 
equilibrium. Thus also the effective couplings are space, i.e. density, dependent. 

The expression of Eq.(6-6) for the effective mass embodies an isospin contri- 
bution from Fock terms even without a direct inclusion of the 5 meson in the 
Lagrangian. As seen from Eqs.6-7 the usual RMF approximation (Hartree 
level) is recovered from the Hartree-Fock results, by changing the effective 
coupling functions fi{i — a,u!,p,5), Eqs.(6-7), to the the explicit coupling 
constants /j. 



6.1.2 Equilibrium properties: the nuclear Equation of State 

The energy density and pressure for symmetric and asymmetric nuclear matter 
and the n, p effective masses are self-consistently calculated in terms of the four 
boson coupling constants, = {:^), i = o-,uj,p,6, and the two parameters 

of the a self-interacting terms, A = and B = [12,13]. Here we will 
present results at the Hartree level. Tfie extension to Fock contributions is 
easily performed following the scheme discussed before. 

From the Lagrangian, Eq.(6-1), the energy-momentum tensor is derived in the 
mean field approximation as 
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-^mluxu;^ - ^mlBxB^]g^,. (6-8) 

with C/(0), the nonhnear potential of the a meson C/(0) = |a0^ + By 
using the field equations for mesons, the equation of state for thermal matter 
is derived, as in [12], with energy density 



i=n,p 



(27r) 



+W) + ^pi + 5^/e3 + ^pJs. (6-9) 



and pressure 



where = y k"^ + M*"^ and the nucleon effective masses are given in Eq.6- 
6. The ni{k) and nj(/c) in Eqs. 6-9, 6-10 are the p,n fermion and antifermion 
distribution functions: 



ni{k) = 



l + exp{(E,*(fc)-K)/T}' 



The effective chemical potentials /x* arc given in terms of the vector meson 
mean fields 



l^t = K - fo:PB =F fppB^ , (6-12) 

where are the thermodynamical chemical potentials /ij = de/dpi^ which at 



zero temperature reduce to the Fermi energies Epi = ^kli + Mf. The baryon 
densities pb and the scalar densities ps are given by (7 is the spin/isospin 
degeneracy) 
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Table 7 

Parameter sets of the QHD models discussed here 



parameter 


NLp 


NLpS 


NLHF 


fa ifm^) 


11.3 


same 


9.15 




6.5 


same 


3.22 


fp ifm') 


1.1 


3.15 


1.9 


fs ifm') 


0.00 


2.4 


1.4 


A (/m-i) 


0.02 


same 


0.098 


B 


-0.004 


same 


-0.021 



r r\^k M* 

Ps-lJj^s^Hk)+n{k))- (6-13) 

At the temperatures of interest here the antybaryon contributions are actually 
negligible. 

The isoscalar meson couplings are fixed from symmetric nuclear matter prop- 
erties at T = 0: saturation density po = 0.15/m~^, binding energy E/A — 
-16MeV, nucleon effective mass M* = 0.7Mn {M^ = 939MeV) and incom- 
prcssibility Ky = 2A0MeV at po- The fitted f^. f^.A.B parameters arc re- 
ported in Tabic 7. They are quite standard for these minimal non-linear RMF 
models, with values compatible with microscopic DBHF estimations [33,34] in 
a wide range of densities. In order to isolate the effects in the isovector channel 
we will review results obtained using the same isoscalar interaction and with 
the p{NLp)— or the p+6{NLp6)— couplings in the isovector part (refs. [12,13]. 
The symmetry energy at saturation is fixed to the Bethe-Weiszacker fourth 
parameter 04 = 30.5 MeV. NLHF stands for the non-linear Hartree-Fock 
scheme described before. 



6.1.3 Symmetry energy 

The symmetry energy in ANM is defined from the expansion of the energy 
per nucleon E{pB,I) in terms of the asymmetry parameter I defined as / = 

PB PB A 

E{pbJ) = '-^^ = E{pB) + E,yMl^ + 0{l') + ... (6-14) 
Pb 
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and so in general 



_ l d^E{pBj) , _ 1 dh 



In the Hartree case an explicit expression for the symmetry energy is easily 
derived [11,12,13]: 



^kin _|_ Tppot 
sym ' sym 



^ E^Z + E^:ln , (6-16) 



where kp is the nuclcon Fermi momentum corresponding to ps, Ep = \J (kp + M*^) 
and M* is the effective nucleon mass in symmetric NM, M* — Mjv — 
A{kF, M*) represents the integral 

Ait,, M-) . j ^3.^^,-^ = 3 - ^) (6-17) 



We remark that A{kFi M*) is very small at low densities, and actually can still 
be neglected up to a baryon density ps — 3po, [12,13]. Then in the density 
range of interest here we can use, to leading order, a much simpler form of the 
symmetry energy, with more transparent 5-meson effects: 



Esym{PB) 



QEl 2 



fp - fs 



M* 



Pb 



(6-18) 



We see that, when the 5 is included, the empirical 04 value actually corresponds 
to the combination [fp — fsi^Y] of the (p, S) couphng constants. Therefore 
if /a 7^ we have to increase correspondingly the p-couphng (see Fig.l of 
ref.[ll]). 

In Table 7 the NLp set corresponds to fs — 0. In the NLpS interaction fs 
is chosen as 2.4/m^, roughly derived from the analysis of ref.[34]. As already 
noted before this choice is not essential for our discussion: the aim of our work 
is to show the new dynamical effects of the 5-meson coupling and to select the 
corresponding most sensitive observables. 

In order to have the same 04 at saturation wc must increase the p-coupling 
constant by a factor three, up to fp = 3.15/m^. Now the symmetry energy 
at saturation density is actually built from the balance of scalar (attractive) 
and vector (repulsive) contributions, with the scalar channel becoming weaker 
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with increasing baryon density. This is clearly shown in Fig. 6-1. This is in- 
deed the isovector counterpart of the saturation mechanism occurring in the 
isoscalar channel for symmetric nuclear matter. From this consideration we 
get a further support to the introduction of the 5-coupling in the symmetry 
energy evaluation. 




Fig. 6-1. p- (open circles) and 5- (crosses) contributions to the potential symmetry 
energy, second and third terms of Eq.(6-16). The dashed line is the approximate 
(5-contribution of Eq.(6-18). The solid line is a linear extrapolation of the low density 
behaviour. 

In Fig. 6-2 we show the total symmetry energy for the different models. At 
subnuclear densities, Pb < Po, in both cases, NLp and NLpd, from Eq.(6-18) 
we have an almost linear dependence of Esym on the baryon density, since 
M* ~ Ep as a good approximation . Around and above po we see a steeper 
increase in the (p -|- 6) case since M*/Ep is decreasing. 













135 










100 






/ / 
/ / 










y 

y 




75 




■ / / 
•/ / 






50 










25 












C 




1 , , , , 1 , , 


. . 1 






1 3 


3 





Fig. 6-2. Total (kinetic -|- potential) symmetry energy as a function of the baryon 
density Dashed line (iVLp). Dotted line {NLp5). Solid line NLHF. 

In conclusion when the 5-channel is included the behaviour of the symmetry 
energy is stiffer at high baryon density from the relativistic mechanism dis- 
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cussed before. This is in fact due to a larger contribution from the p relative 
to the 5 meson. We expect to see these effects more clearly in the relativistic 
reaction dynamics at intermediate energies, where higher densities are reached 
(see Sect. 8. 
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Fig. 6-3. Baryon density variation of the isovector effective coupUng when the Fock 
terms are included, model NLHF. 

When the Fock terms are evaluated the new "effective" couplings Eqs.(6-7) 
naturally acquire a density dependence. This is shown in Fig. 6-3 for the isovec- 
tor terms. The decrease of the p coupling at high density accounts for the slight 
softening of the symmetry energy. Fig. 6-2, in agreement with DBHF expec- 
tations, see [34,232]. Details of the calculation can be found in Refs. [36,12]. 

In Fig. 6-4 we show the Equation of State (energy per nucleon) for pure neutron 
matter (/ = 1) obtained with the two parameter Sets, NLp and NLpS. The 
values, in particular for the NLpS case, are in good agreement with recent 
non-relativistic Quantum-Monte-Carlo variational calculations with realistic 
2— and 3— body forces [47]. The inclusion of the (5-coupling leads to a larger 
repulsion at baryon densities roughly above l.Spo- This would be of interest 
for the structure of neutron stars. It is also relevant for the possibility of a 
transition to new forms of deconfined nuclear matter [241]. 



6.1.4 Symmetry Pressure and Symmetry Incompressibility 

From the previous analysis we expect to see also interesting 5-effects on the 
slope {symmetry pressure) and curvature {symmetry incompressibility) of the 
symmetry energy around po- These quantities and their relevant physical 
meaning were already discussed in Sect. 2 in a non-relativistic framework. 

From Eq.(6-18) we get the potential contribution to the density variation of 
Esym (after some algebra): 
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Pb/Po 



Fig. 6-4. EOS for pure neutron matter. Dashed line: Nip. Solid line: NLp6). 



Ipot 



U^^y-- - i^^i (6-19) 

Around normal density po the first term is fixed by the 04 value (cfr. Eq.(6-18)). 
The second term (always positive since [DM* / dps) < 0) gives a net increase 
of the slope, due to the 5-field introduction, as seen also in Fig.6-2. To be 
more quantitative, with our parametrization we get a potential contribution 
to the slope L of ASMeV from the first term and a genuine 5-contribution of 
about 20MeV from the second one. When we include also the kinetic part 
(from the first term of Eq.(6-18)) we have a total slope parameter going from 
L{p) = +UMeV to L{p + 5) = +103Me1/. 

We note again that the slope parameter, or equivalently the Symmetry Pres- 
sure Psym = PoL/3, is of great importance for structure properties, being 
linked to the thickness of the neutron skin in n-rich (stable and/or unstable) 
nuclei [242,243,244,30] , and to the position of the drip-line. Moreover the 
same parameter gives an estimate of the shift of the saturation density with 
asymmetry (at the lowest order in /^), easily obtained from a linear expansion 
around the symmetric value po{I = 0) as shown in Sect. 2, [52]. 

It is instructive to perform a similar analysis for the curvature parameter 
Kgym- Now the potential contribution is exclusively given by the 5-meson, with 
a definite positive sign, as we can see from the previous discussion. It is a large 
effect on the total since the kinetic part of the symmetry incompressibility is 
quite small [245]. Compared to non-relativistic effective parametrizations, in 
Sect. 2, when we add the 5-meson, we move from a linear to a roughly parabolic 
Ps-dependence of the symmetry energy. 

With our parameters we pass from a Kgymip) = +7MeV (only kinetic) to a 
Ksym{p + 5) = -\-120MeV [245]. Thus this quantity appears extremely inter- 
esting to look at experimentally. The problem is that the effect on the total 
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incompressibility of asymmetric matter, that likely could be easier to measure, 
is not trivial. What really matters for the total incompressibility is the combi- 
nation {Ksym — QL), as shown in Sect. 2, with the possibility of a compensation 
between the two terms. Just by chance this is actually what is happening in 
our calculations since for the above combination we get —A97MeV in the case 
of only p-coupling, and —bOAMeV when we add also the 5-field. 

6.1.5 Finite Temperature Effects 

In the temperature range of interest in this paper, below the critical tem- 
perature Tc of the liquid-gas phase transition of the order of 15 — 16 MeV, 
temperature effects on the symmetry properties are not expected to be large. 
Indeed the contributions of antifermions, that could modify all the terms with 
the scalar densities are still very reduced. For both cases, p and {p + 5), the 
temperature variation of symmetry energy is quite small. We have a reduc- 
tion mainly coming from the kinetic contribution due to the smoothing of the 
n/p Fermi distributions. We note that this result is in full agreement with 
relativistic Brueckner-Hartree-Fock calculations [247]. 

The temperature effect on the nucleon mass splitting, given by a difference of 
n/p scalar densities Eq.(6-6), is even smaller. In the following subsection we 
will study in detail the phase diagram of heated asymmetric nuclear matter, 
focussing in particular on the instability regions. 

6.2 Mechanical and chemical instabilities 

Heavy-ion collisions can provide the possibility of studying equilibrated nu- 
clear matter far away from normal conditions, i.e. to sample new regions of 
the NM phase diagram. In particular the process of multifragmentation al- 
lows to probe dihite nuclear matter at finite temperatures. In the symmetric 
case we expect to see a phase transition of first order of liquid-gas type, as 
suggested from the very first equations of state built with effective interactions 
[248,249,250]. 

As discussed extensively in Sect. 3 in a non- relativistic framework, for asym- 
metric nuclear matter a qualitatively new feature in the liquid-gas phase tran- 
sition is expected, the onset of a coupling to chemical instabilities (compo- 
nent separation). This will show up in a novel nature of the unstable modes, 
the mixture of density and charge fluctuations leading to an Isospin Distilla- 
tion. Indeed, equilibrium thermodynamics as well as non-equilibrium kinetics 
both predict that an asymmetric system will separate into more symmetric 
larger fragments ("liquid" phase) and into neutron- rich light fragments ("gas" 
phase) . Since the effect is driven by the isospin dependent part of the nuclear 
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equation of state, here we will look at the influence of the 5-coupling on this 
new liquid-gas phase transition. 

As discussed earlier the instability condition of a two-component, n/p, ther- 
modynamical system is given by 



where P is the pressure, Hp is the proton chemical potential and y the proton 
fraction Z/A, related to the asymmetry parameter I — 1 — 2y. Eq.(6-20) is 
equivalent to constrain the free energy to be a convex function in the space 
of the n.p density oscillations, 6pn,Spp. In charge symmetric matter isoscalar 
(total density) 6pn + Spp and isovector (concentration) 6pn — Spp oscillations 
are not coupled and we have two separate conditions for instability: 



—— < 0, mechanical, i.e. vs. density oscillations, and 

-^^] < 0, chemical, i.e. vs. concentration oscillations. (6-21) 

Jtp 



In asymmetric matter the isoscalar and isovector modes are coupled and the 

two separate inequalities do not anymore mantain a physical meaning, in the 
sense that they do not select the nature of the instability. Inside the general 
condition Eq.(6-20) the corresponding unstable modes are a mixing of density 
and concentration oscillations, very sensitive to the charge dependent part of 
the nuclear interaction in the various instability regions [20]. 

In dilute asymmetric NM (n-rich) the normal unstable modes for all realistic 
effective interactions are still isoscalar-like, i.e. in phase n — p oscillations 
but with a larger proton component. This leads to a more symmetric high 
density (liquid) phase everywhere under the instability line defined by Eq.(6- 
20) and conscquentely to a more neutron-rich gas {Isospin Distillation) . Such 
"chemical effect" is driven by the increasing symmetry repulsion going from 
low to roughly the saturation density and so it appears rather sensitive to the 
symmetry energy of the used effective interaction at subnuclear densities. 

In this section we study the effect of the ^-coupling on the instability region 
given by Eq.(6-20) in dilute asymmetric nuclear matter and on the structure 
of the corresponding unstable modes. We start from an identity valid for any 
binary thermodynamical system 
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Fig. 6-5. Limits of the instability region in the T,pB plane for various asymmetries. 
Dashed lines: NLp). Solid lines: NLpS) [a being the asymmetry parameter]. 



where pg, pq {q — n,p) are respectively neutron/proton chemical potentials 
and densities. Thus from the chemical potential on each isotherm, Eq.(6-12), 
we can easily compute the limits of the instability region in the T, ps plane 
for dilute asymmetric NM in the two choices, without and with the 5-meson. 

The results are shown in Fig.6-5. The inclusion of a 5-field (solid hues) appears 
not to affect much the instability limits even at relatively large asymmetry 
/ = 0.8 (A^ ~ 9Z). We notice just a small reduction and a shift to the left 
(lower densities) of the whole region: this can be related to a slightly larger 
symmetry repulsion. We can understand the relatively small 5-effect on the 
stability border by remembering that for low densities, well below po, the 
symmetry term has roughly the same linear behaviour in both (p) and (p + 5) 
schemes, fixed by the 04 parameter (see the discussion after Eq.(6-18)). Only 
for very large asymmetries it appears relatively easier in the 6 case to be in 
the stable liquid phase. 

A larger difference can be seen in the behavior of the quantity Eq.(6-20) inside 
the instability region. This is plotted in Fig. 6-6 for various asymmetries at zero 
temperature. The solid curves (with 5-couphng) are systematically above the 
dashed ones, possible signature of a weaker instability. In fact this is not the 




(6-22) 
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Fig. 6-6. The quantity Eq.(6-20) inside the instability region at T = and various 
asymmetries. Dashed hues: NLp. Solid lines: NLpS. 

case since, as shown in the following, the 5-meson is not affecting the low 
density unstable modes. In order to better understand the origin of this effect 
we study the structure of the corresponding unstable modes. We follow the 
Landau dispersion relation approach to small amplitude oscillations in Fermi 
liquids, as in Sect. 3. For a two component {n,p) matter the interaction is 
characterized by the Landau parameters Fq'^ , {q, q') = {n,p) defined as 

N,{T)^^S,,,, + Fr' (6-23) 

where Nq{T) represents the single particle level density at the Fermi energy. At 
zero temperature it has the simple form (H — c — 1) Ng = -^^^^r^, q — n,p. 

In the symmetric case (Fg"" = F^^, Fq^ = Fq^"), the Eqs.(6-21) correspond to 
the two Pomeranchuk instability conditions 



= F^"" + Fq"^ < -1 mechanical 
= Fq™ - Fo"P < -1 chemical. (6-24) 

From the dispersion relations Fq will give the properties of the density (isoscalar) 
modes while F^ is related to the concentration (isovector) modes. For asym- 
metric NM we have corresponding generalized Landau parameters Fq^, Fq^ 
which characterize the asymmetric collective response. They can be expressed 
as fixed combinations of the Fq '"^ for each baryon density, asymmetry and 
temperature. This transformation reduces Eq.(6-22) to a "diagonal" form [20] 
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Fig. 6-7. Behaviour of the generahzed Landau parameters Fq^ inside the instabihty 
region at zero temperature and asymmetry I = 0.5 {N = 3Z). Dashed hues: NLp. 
SoUd hues: NLpS. 



l + Foy(l + F»J 



OgJ 




(6-25) 



As already discussed, in the unstable region of dilute asymmetric NM we have 
isoscalar-like unstable modes with 1 -|- F^g < 0, while the combination 1 -|- F^g 
will always stay positive. In Fig. 6-7 we report a full calculation of these two 
quantities in the unstable region at zero temperature, for asymmetry I = 0.5, 
with and without the (5-coupling. The 5-meson almost does not affect at all the 
unstable mode, given by the F^g parameter. Hence the limits of the instability 
region, for I — 0.5, are not changed (see the Figs. 6-5, 6-6). We have a larger 
effect on the Fq^ parameter which describes "stable" is ovector-like modes, that 
actually can propagate as good zero-sound collective motions since Fgg > 0. 
This can be expected from the isovector nature of the 5-meson. From here we 
get the main differences seen in Fig. 6-6 for the Eq.(6-20) inside the instability 
region, just the product of the two quantities plotted in Fig.6-7 (from Eq.(6- 
25)). 

Another interesting aspect of the comparison between Fig. 6-6 and Fig.6-7 is 
the shift of the "maximum instability" density region. From the thermody- 
namical condition reported in Fig. 6-6 it seems that the largest instability (the 
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Fig. 6-8. Neutron and proton effective masses vs. the baryon density for / = 0.5 
(A^ = 3Z), NLpS parameters. The dashed hne corresponds to symmetric nuclear 
matter. 

most negative value) is around pB = In fact from Fig.6-7 we see 

that the fastest unstable mode, corresponding to the most negative Pomer- 
anchuk condition for I + Fq^, is actually present for more dilute matter, around 
Pb — 0.02/m~^. This shows the relevance of the hnear response analysis. 

Finally the fact that the 5-coupling mostly affects the stable, isovector-like, 
modes is of great interest for possible effects on the isovector Giant Dipole Res- 
onances studied around normal density within the RMF approach in asym- 
metric systems. This will be shown in the next Section on the relativistic linear 
response in asymmetric matter. 



6.3 Nucleon Effective Mass Splitting 

An important qualitatively new result of the (5-meson coupling is the n/p- 
effective mass splitting in asymmetric matter [11,12], see Eq.(6-6). In Fig.6- 
8 we report the baryon density dependence of the n/p effective masses for 
/ = 0.5 (A^ = 3Z) asymmetry, calculated with our NLp6 parameters, com- 
pared to the symmetric case [246]. We see a splitting of the order of 15% at 
normal density po, increasing with baryon density. Unfortunately, from the 
present nuclear data we have a very limited knowledge of this effect, due to 
the low asymmetries available. However we can expect important effects on 
transport properties ( fast particle emission, collective flows) of the dense and 
asymmetric NM that will be reached in Radioactive Beam collisions at inter- 
mediate energies. 




116 



The sign itself of the sphtting would be very instructive, since it is a quite 
controversial point (see the comments at the end of Sect. 2). As we can see 
from Eq.(6-6), in a relativistic approach in n-rich systems we have a neutron 
effective mass always smaller than the proton one. We note again that a de- 
creasing neutron effective mass in n-rich matter is a direct consequence of 
the relativistic mechanism for the symmetry energy, i.e. the balance of scalar 
(attractive) and vector (repulsive) contributions in the isovector channel. 



6.3.1 Dirac and Schrddinger Nucleon Effective Masses in Asymmetric Mat- 
ter 

The prediction of a definite m* < m* effective mass splitting in RMF ap- 
proaches, when a scalar 5 — like meson is included, is an important result 
that requires some further analysis, in particular relative to the controversial 
predictions of non-relativistic approaches, see Sect. 2. 

Here we are actually discussing the Dirac effective masses mj)(n,p), i.e. the 
effective mass of a nucleon in the in- medium Dirac equation with all the meson 
couplings. The relation to the Schrddinger eSective masses m*g{n,p), i.e. the 
"k-mass" due to the momentum dependence of the mean field in the non- 
relativistic in-medium Schrodinger equation is not trivial, see [237,28,252,253]. 
We will extend the argument of the refs. [252,253] to the case of asymmetric 
matter. 

We start from the simpler symmetric case without self-interacting terms. The 
nucleon Dirac equation in the medium contains the scalar self-energy = 
~faPs s-nd the vector self-energy (fourth component) Eq = fuiPs and thus the 
corresponding energy-momentum relation reads: 

(e + m - Eo)' =p' + {m + E,)^ = / + ^g_2g) 

i.e. a dispersion relation 

e = -m + Eo + \ff^^^ (6-27) 



Prom the total single particle energy E — e + m expressed in the form E — 
\Jk1^-\- m'^, where koo is the relativistic asymptotic momentum, using Eq.(6- 
26) we can get the relation 



2m 



e + 



2m 



|^ + S. + Eo + -^(E^-E^) + ^6 
2m 2m m 



^ + Ueff{pB,Ps,e) (6-28) 
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i.e. a Schrodinger-type equation with a momentum dependent mean field that 
with the dispersion relation Eq.(6-27) is written as 



eff 



m 2m 2m m 



E 
m 

-i2 I v'S 



(6-29) 



D 



The relation between Schrodinger and Dirac nucleon effective masses is then 



m ^ m 



= , s„ = ^*D _ , ^ , ^ (6-30) 



Since at saturation the two self-energies are roughly compensating each other, 
Es + Eo — —50MeV the two effective masses are not much different, with the 
5" — mass slightly larger than the D — mass. 

In the case of asymmetric matter, neutron-rich as always considered here, we 
can have two cases: 



• Only p meson coupling 

Now the scalar part is not modified, we have the same scalar self energies 
Eg for neutrons and protons and so the same Dirac masses. The vector self 
energies will show an isospin dependence with a new term TfpPBs, repulsive 
for neutrons (— sign, since we use the definition p(B,s)3 = P{B,s)p ~ P{B,s)n)- 
As a consequence we see a splitting at the level of the Schrodinger masses 
since Eq.(6-30) becomes 

TTl 

m*s{n,p) = (6-31) 

m + [hs + ho)sym T JpPB3 

in the direction of m* < m* (here and in the following upper signs are 
for neutrons). 

• p + 5 coupling 

The above splitting is further enhanced by the direct effect of the scalar 
isovector couphng, see Eq.(6-6). Then the Schrodinger masses are 

TTl 

The new term in the denominator will further contribute to the m*„ < 
m*p splitting since we must have fp > fs in order to get a correct symmetry 
parameter 04. The effect is larger at higher baryon densities because of the 
decrease of the scalar ^53, due to the faster reduction of psn- 
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Actually in a non-relativistic limit wc can approximate in Eq.(6-28) directly 
the energy e with the asymptotic kinetic energy leading to the much simpler 
relation: 

TJi 

m*s= ~ m - So = - (S^ + Sq). (6-33) 

In the case of asymmetric matter this leads to the more transparent relations 

• Only p meson coupling 

m*s{n,p) = m*jj- (S^ + Soj^ym ± /pPsa (6-34) 

• p-\- 5 coupling 

= ^*Dsyrn " (^^ + sym ± (/pPB3 " f&PSz) (6-35) 

In conclusion any relativistic effective field model will predict the definite 
isospin splitting of the nucleon effective masses m* < m* in the non-relativistic 
limit. We expect an increase of the difference of the neutron/proton mean field 
at high momenta, with important dynamical contributions that will enhance 
the transport effects of the symmetry energy. 



6.3.2 The Dirac-Lane Potential 

In the non-relativistic limit of the Eq.(6-28) we can easily extract the neutron- 
proton mean optical potential using the isospin dependence of the self-energies 

=F fppB3 , 

^s,q = ^s,sym ± fsPS3 , (6-36) 

With some algebra we get a compact form of the Dirac — Lane Potential 

Un-Up 



U Dirac— Lane 

Po 



/,(l-5^)-/5^(l + 5^) +/p^e (6-37) 



21 

Po 
— ( 

m 



It is interesting to compare with the related discussion presented in Sect. 2 for 

the non-relativistic effective forces, mainly of Skyrme-like form. First of all 
wc predict a definite positive E — slope given by the quantity fp—, which is 
actually not large within the simple RMF picture described here. This is a 
obvious consequence of the fact that the "relativistic" mass sphtting is always 
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Fig. 6-9. Energy dependence of the Dirac-Lane potential in the RMF picture (solid: 
NLp; dashed: NLp6) and in the phenomenologic Dirac Optical Model of Madland 
et al. [254,255], see text. 



in the direction m* < m*. The realistic magnitude of the effect could be larger 
if we take into account that some explicit momentum dependence should be 
included in the scalar and vector self energies, as discussed in ref. [146]. An 
interesting point in this direction comes from the phenomenological Dirac 
Optical Potential {Madland — potential) constructed in the refs. [254,255], 
fitting simultaneously proton and neutron (mostly total cross sections) data 
for collisions with a wide range of nuclei at energies up to 100 MeV. Recently 
this Dirac optical potential has been proven to reproduce very well the new 
neutron scattering data on ^osp^ gg ]\^qY [256] measured at the Svendberg 
Laboratory in Uppsala. 

The phenomenological Madland — potential has different implicit momen- 
tum dependences (exp/log) in the self-energies, [254,255]. In Fig. 6-9 we show 
the corresponding energy dependences for the Dirac — Lane potentials, com- 
pared to our NLp and NLp6 estimations. When we add the (5-field we have 
a larger slope since the p — coupling should be increased. The slopes of the 
phenomenological potentials are systematically larger, interestingly similar to 
the ones of the Skyrme — Lyon forces of the Fig. 2-6 of Sect. 2. We note that 
in the Madland — potential the Coulomb interaction is included, i.e. an extra 
repulsive vector contribution for the protons. The first, not energy dependent, 
term of Eq.(6-37) is directly related to the symmetry energy at saturation, 
exactly like in the non-relativistic case, see Eq.(2-22) of Sect. 2. 
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The conclusion is that a good systematic measurement of the Lane potential 
in a wide range of energies, and particularly around/above 100 MeV, would 
answer many fundamental questions in isospin physics. 
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7 Collective modes of asymmetric nuclear matter in the relativistic 
approach 



7.1 Linear Response Equations 



We discuss here collective oscillations that propagate in cold nuclear matter 
under the influence of the mean field dynamics. These studies can be con- 
sidered a relativistic extension of the method introduced by Landau to study 
liquid-^ife [73,257,71] and recently applied to investigate stable and unstable 
modes in nuclear matter [72,50,51]. The starting point is the relativistic kinetic 
transport equation of Sect. 6, Eq.(6-5). We look for solutions corresponding to 
small oscillations of the Wigner function F{x,p) around the equilibrium value. 
Therefore we put 

F{x,p) = H{p) + G{x,p) (7-1) 



where H{p) is the distribution function at equilibrium and G{x,p) represents 
its fluctuations, ref.[13]. 

In the linear approximation, i.e. neglecting terms of second order in G{x,p), 
the equations for the Wigner functions become 



(1-'-A){F{x)tFs{x))H^M > (7-2) 

for neutrons and protons (i = n,p, upper/lower signs), where M* — M — 
fcrPs^fs PS3- The quantities F{x) and F3{x) are the isoscalar and the isovector 
components of the self-consistent field: 



F{x) = -^LG{x) + ^L^,G>^{x) - 8^psG{x) - 8^psG^{x) 

ops OJn 

-S^psGsix) - S^psGt^ix) + 8^^,fGix) , (7-3) 

OpS3 Oj3n Ops 



-8fsGs(x) + 8/,7^G^(a;) - sf^ pss G{x) + 8^^,j^G{x) .(7-4) 

ops ops 



Following the notation of Sect. 6 the coupling functions fi are in general den- 
sity dependent due to the Fock contributions. The Hartree approximation is 
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recovered by vanishing all the derivative terms in the quantities F{x) and 
F^ix) (except since still = ^ps)/ps)- 

We obtain equations for the collective oscillations by multiplying Eqs. (7-2) by 
7a. After performing the traces, we equate to zero both the real and imaginary 
parts of the result [222,77]. Furthermore, by Fourier transforming and inte- 
grating over four-momentum, we get the set of equations for the scalar and 
vector fluctuation of each species, [13]. We note that the formalism developed 
here in Hartrcc-Fock approximation includes a linear response theory valid also 
for any other approach to QHD beyond RMF. Namely when we consider in 
general a dependence on the baryon density of all meson-nucleon couplings, 
hke in the Density Dependent Hadronic (DDH) model [258,259,220,260], in- 
spired by the Dirac-Brueckner-Hartree-Fock (DBHF) approach to account 
for many-body correlations. 

The set of equations developed in the Hartree-Fock approximation includes 
the ones corresponding to the usual Hartree approximation (RMF). As al- 
ready mentioned, it is easily obtained by considering the coupling /j to each 
channel equal to the coupling constant of the corresponding meson. The result 
is of appreciable plainer structure due to the constant value of all couplings, 
except fa [261]. Since the physics results become more transparent, in the fol- 
lowing we will stay in the Hartree scheme, keeping well in mind that the Fock 
contributions can be easily included. In fact they amount to having some extra 
contributions in the various interaction channels without qualitative modifi- 
cations of the physical response. 

The normal collective modes are plane waves, characterized by the wave vector 
{k^ ~ (A;°, 0, 0, 1 k |)). They are determined by solving a set of homogeneous 
linear equations. The solutions correspond only to longitudinal waves and do 
not depend on and |k| separately, but only on the ratio 




The zero-sound velocities are given by those values of Vg for which the rel- 
evant determinant of the set of equations vanishes, i.e. by the dispersion 
relations. Correspondingly the neutron/proton structure of the eigenvectors 
(normal modes) can be derived. Again it should be noted that in asymmetric 
nuclear matter isoscalar and isovector components are mixed in the normal 
modes. Here this can be argued by the fact that in each linearized equation 
(7-2) both proton/neutron densities and currents are appearing, [13]. However 
we also recall that one can still identify isovcctor-like excitations as the modes 
where neutrons and protons move out of phase, while isoscalar-like modes are 
characterized by neutrons and protons moving in phase [50,20]. 
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1.2 The Role of Scalar/ Vector Fields in the Dynamical Response 



Before showing numerical results for the dynamical response of asymmetric 
nuclear matter in various baryon density regions and using the different effec- 
tive interactions, we would like to analyse in more detail the structure of the 
relativistic linear response theory in order to clearly pin down the role of each 
meson coupling. 



7.2.1 Isovector Response 

One may perhaps expect that once the asymmetry parameter 04 is fixed, the 
velocity of sound is also fixed [263] . On the other hand our results will clearly 
show a different dynamical response, e.g. with or without the (5-meson channel, 
for interactions which give exactly the same 04 parameter. Essentially this 
is due to the fact that the response depends on the slope of the symmetry 
energy for isovector oscillations around the saturation point. In order to get 
a complete quantitative understanding of this effect we will first consider the 
case of symmetric nuclear matter in the Hartree scheme, where the dispersion 
relations assume a particularly transparent analytical form. 

In this case it is also possible to decouple the collective modes into pure 
isoscalar and isovector oscillations, [50,20] as also seen in Sect. 3. After a 
straightforward rearrangement we have a dispersion relation for the isovec- 
tor modes. Now in order to find the zero-sound velocity one has to evaluate 
the determinant of a 2 x 2 matrix (and not a 4 x 4, as in the asymmetric NM). 
After a straightforward rearrangement we have a dispersion relation [13] for 
the isovector modes 



1 + Nf 



M*' , , PS 



/.(I - vt) - fs^ [1 - fsA{kF, M*) -fpj^v 



,2 



^(s)=0 .(7- 



Here Np — ^^^fi is the density of states at the Fermi surface and s = Vg/vp. 
ip{s) is the usual Lindhard function of the Landau Fermi Liquid response 
theory: 



s-1 



The quantity A{kF, M*) is the same integral discussed in the Sect. 6, Eq.(6-17). 



In fact the structure of the dispersion relation is the same for the isoscalar 
excitations, of course one has to replace the isovector fluctuations {SpBs, Spss) 
with the isoscalar ones {Sps, Sps), and the coupling constants of isovector 
mesons with those of the isoscalar mesons. 
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At this point we can make the following approximation 



2 2 



to evaluate the expression inside the square brackets in Eq.(7-5). Looking at 
numerical results shown later (Fig. 7-2) this is a good approximation to within 
3% for all effective interactions considered. As in Sect. 6 we also neglect the 
integral A{kp, M*). Then Eq.(7-5) assumes the transparent form: 



1 + 



6^ 



ppot 
sym 



fp kp 
2 Ef 



M* 



PS PB 



(7-6) 



The potential part of the symmetry energy explicitly appears in the dispersion 
relations, but together with an important correction term which exhibits a 

different fp, fs structure with respect to that of -Bf"^, Eq.( ) of Sect. 6. We 

can easily have interactions with the same value at normal density but 
with very different isovector response. E.g. when including the S channel we 
have to increase the /p coupling in order to have the same a^, as discussed in 
Sect. 6, but now the "restoring force" (coefficient of the Lindhard function in 
the Eq.(7-6)) will be reduced. 

A similar effect has been pointed out in a detailed non-relativistic Skyrme — 
RPA study of the Giant Dipole Resonance in heavy nuclei (^°^P6) using effec- 
tive interactions with various isovector terms [262]. A separate sensitivity of 
the average resonance frequencies on the symmetry energy 04 and on its slope 
has been found. In a covariant scheme we can see from Eq.(7-6) that such be- 
haviour can be achieved only by using two isovector fields, at the lowest order. 
This result shows more generally that a dynamical observable can be more 
sensitive to the microscopic structure of the isovector interaction than static 
properties. For instance in a careful study of the neutron distributions, [30], 
it is clearly shown that these "equilibrium" observables are almost equally 
correlated to value, slope and curvature of the symmetry term. 



7.2.2 Isoscalar Response 

We already noted that for symmetric NM there is a close analogy between the 
isoscalar and the isovector response in the RMF approach. In the isoscalar 
degree of freedom the compressibility will play the same role of the symmetry 
energy in the dispersion relation equations. Also in this case we will have an 
important correction term coming from the interplay of the scalar and vector 
channel. 
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The Eq.(7-5) now becomes [261]: 

M*2 



1 + Nf 



Ivf / n \ 

fM - v's) - (l - UA{kF, M*) - ^ vl) 



V{s) = 0(7-7) 



that can be reduced with the same approximations to the isoscalar equivalent 
of the Eq.(7-6): 



1 + 



E*p 
3 k'p 



T^pOt 



9/a 



Ef 




(^(s)=0, 



(7-8) 



Here the K^l^ is the potential part of the nuclear matter compressibility that 
in the Hartree scheme has the simple structure [261] (see Eq.(16) of ref. [263]) 



Knm{pi 



3 kp 
~El 



+ 9 



foj fa 



M* 



PB — -f^ATM + ^NM 



(7-9) 



With this analogy, the previous discussion can be transferred to isoscalar os- 
cillations with the role of Egym now "played" by the compressibility. In fact in 
the isoscalar sector one always takes into account both the scalar and vector 
fields in any RMF model. However, the coupling constant f^, can assume very 
different values depending on the chosen value for effective masses Mq. This 
is easy to understand since in the RMF limit the saturation binding energy 
has the simple form 

E/A{po)^E*F + f^po-MN 

where is the bare nucleon mass. Thus in order to have the same satu- 
ration values of po, E/A{po), when we decrease Mq we have to increase f^. 
We then come to the natural conclusion that two EOS with different effec- 
tive masses, even if the compressibilities are the same, are expected to have 
different dynamical behaviour. This is a very general feature present also in 
non-relativistic approaches. 

From studies in RMF on monopole resonances in finite nuclei it seems that a 
higher value of compressibility is required with respect to non-relativistic cal- 
culations. Many authors state that this certainly demands for a clarification 
[225,226,264,265]. Even if the monopole resonance is not directly connected 
to the isoscalar collective mode in nuclear matter, our discussion nicely sug- 
gests to look at the interplay between effective mass and compressibility. For 
example we can estimate by means of the dispersion relation Eq.(7-8) that 
we can have a shift between the compressibility and the "effective compress- 
ibility" of the order of ~ 100 MeV. Therefore an effective interaction e.g. 
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with K ~ 300 AAl eV can reproduce the same frequencies of another one with 
K ~ 200 AM eV (and a shghtly larger M^). 



7.2.3 Landau Parameters 



Next wc discuss the rclativistic equations for collective modes in terms of the 
Landau parameters. Interesting features will appear from the comparison to 
the analogous non-relativistic case of Sect. 3. We will focus first on the isovector 
response, but, as shown before, the structure of the results will be absolutely 
similar in the isoscalar channel. 

The general non-relativistic expression for dispersion relation of isovector 
modes can be found in ref.[257]: 
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where Fq is the "isovector" combination of the Landau Fq parameters for 
neutrons and protons Fq = F^" — Fq^, already introduced in Sect. 2. F^ are 
the equivalent quantities for the momentum dependent part of the mean field. 
In the relativistic approach, for symmetric nuclear matter, we have: 



M*2 

F" — F — Fi 



Ef l + fsA{kF,M*) 

1 + -^^pVp 

where F^ = Npfi {i = p, 6) with Np = 2iV„,p. Note that the F^ contribu- 
tion results only from the vector coupling. By using the expression of F^°^ 
(Eq.(6-18) of Sect. 6), we can write Eq.(7-10) in the same form of Eqs.(7-5, 
7-6). The result is a similar expression but without the term in fs inside the 
brackets in Eq.(7-6). As said before this is not the leading term, however, 
around saturation density it amounts to about 10% of the total correction. 



Moreover turning to the analogy with isoscalar channel, the corresponding 
term is now the coupling of the a field, which is much larger. This purely rel- 
ativistic contribution could be up to a 20%. We underline this point because 
generally the linear response in RMF is discussed calculating the Landau pa- 
rameters and then using these values directly in the non-relativistic expression 
for collective modes [263,266]. 

From the analysis in terms of the Landau parameters, we can describe the 
effect of the scalar-vector coupling competition previously discussed in the 
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following way. The symmetry energy fixes the Fq , in fact: 
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(7-12) 



but in the dynamical response also the F" enters, linked to the momentum 
dependence of the mean field, mostly given by the vector meson coupling. The 
results are completely analogous in the isoscalar channel, with the compress- 
ibility given by 



with the "isoscalar" combination Fq = F^'^ + Fq'^ . The relativistic forms of 
the isoscalar Landau parameters are exactly the same as in Eq.(7-ll), just 
substituting the 5, p coupling constants for the cr, u) ones [261] . 

7.3 Isovector Collective Modes in Asymmetric Nuclear Matter 

We now turn to Asymmetric Nuclear Matter. We discuss first results for the 

isovector collective oscillations, which are driven by the symmetry energy 
terms of the nuclear EOS. The aim is mainly to investigate the effect of the 
scalar-isovector channel. This is normally not included in studying the isovec- 
tor modes, or generally the properties of asymmetric matter in a relativistic 
approach, while it should be naturally present on the basis of the analysis pre- 
sented in Sect. 6. Moreover we emphasize again that a Hartree-Fock scheme 
leads to the presence of a scalar-isovector channel in any case, even without 
the inclusion of the 5 -meson field [238,36]. 

We will first discuss results obtained in the Non-Linear Hartree scheme (NLH) 
including either both the isovector (p + 6) mesons {NLp5 parametrization) or 
only the p meson, (NLp). Even though the Hartree approximation has a sim- 
pler structure, it contains all the physical effects we want to point out. Finally 
from the complete Hartree-Fock (NLHF) calculations we will confirm the 
dynamical contribution of the scalar isovector channel. For NLH calculations 
we use the paramctrizations of Table 7, Set{NLp) and Sct{NLpS), of Sect. 6. 
In the Hartree-Fock case the coupling constant fs is adjusted to the value 
fs{po) = 2.0 fm^ of the NLHF model, Eqs.(6-7) of Sect.6. 

In Fig.7-la we show the sound velocities in the Hartree approximation as 
a function of the asymmetry parameter / for different baryon densities. We 
actually plot the sound velocities in units of the neutron Fermi velocities. This 
is physically convenient: when the ratio is approaching 1.00 we can expect 
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that this "zero" sound will not propagate due to the strong coupling to the 
"chaotic" single particle motions ("Landau damping"). This quantity then 
will also directly give a measure of the "robustness" of the collective mode. 
Dotted lines refer to calculations including (p + 5) mesons, long-dashed lines 
correspond to the case with only the p meson. Calculations are performed for 
Pb = Po and = 2 po- We note that the results of the two models differ 
already at zero asymmetry, / = 0. At normal density (po curves), in spite of 
the fact that the symmetry energy coefficient, = Esym{Po), is exactly the 
same in the two cases, significant differences are observed in the response of 
the system. Prom Fig. 7-1 (a) we expect a reduction of the frequency for the 
bulk isovector dipole mode in stable nuclei when the scalar isovector channel 
{6 — like) is present. Moreover we note that, in the NLp case, the excitation of 
isovector modes persists up to higher asymmetries at saturation density. These 
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Fig. 7-1. Isovector-like modes: (a) Ratio of zero sound velocities to the neutron 
Fermi velocity Vpn as a function of the asymmetry parameter I for two values 
of baryon density. Long dashed line: NLp. Dotted line: NLpS. (b) Corresponding 
ratios of proton and neutron amplitudes. All lines are labelled with the baryon 
density, po = 0.16/m~^. The full circles in panel (b) represent the trivial behaviour 
of —(pp/pn) vs. asymmetry parameter. 



are non-trivial features, related to the different way scalar and vector fields 
are entering in the dynamical response of the nuclear system. Such behaviours 
are therefore present in both collective responses, isoscalar and isovector, as 
seen in the previous subsection. 
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Differences are observed even at = 2po, where however also the symmetry 
energy is different. At higher densities a larger Egym is obtained in the case 
including the 6 meson and this leads to a compensation of the effect observed 
at normal nuclear density. In particular, at higher asymmetries / the collective 
excitation becomes more robust for NLp6. Differences are observed also in the 
"chemical" structure of the mode, represented by the ratio Spp/Spn, which is 
plotted in Fig.7-l(b). The ratio of the out of phase n/p oscillations is not fol- 
lowing the ratio of the n/p densities for a given asymmetry, as shown by the 
full circles in the figure. We systematically see a larger amphtude of the neu- 
tron oscillations. The effect is more pronounced when the 6 (scalar-isovector) 
channel is present (dotted lines). 



7.3.1 Disappearance of the Isovector Modes 

For asymmetric matter it is found that for all interactions with increasing 
baryon density the isovector modes disappear: we call the density, where this 
occurs, pg"°**. E.g. from Fig.7-l(b) we see that the ratio 5pp/Spn tends very 
quickly to zero with increasing baryon density, almost for all asymmetries. 
Around this transition density we expect to have an almost pure neutron wave 
propagation of the sound. In Figs. 7-2 and 7-3 we show the results of the NLp 
case, but the effect is clearly present in all the models. 

In Fig. 7-2 we give the zero-sound velocities as a function of density for a 
definite asymmetry. For symmetric matter we have a real crossing of the two 
sound velocities, isoscalar and isovector, as shown in Fig.7-2(a). Above pg""** 
the isoscalar mode is the most robust. 

For asymmetric matter we observe a transition in the structure of the propa- 
gating normal mode, from isovector-like to isoscalar-like, Fig.7-2(b,c). Similar 
effects have been seen in a non-relativistic picture [50]. This mechanism is 
analysed in Fig. 7-3 for the ratio of proton to neutron amplitudes (including 
the sign) . In fact we see that the proton component of the propagating sound 
is quite small in a relatively wide region around the "transition" baryon den- 
sity, a feature becoming more relevant with increasing asymmetry, see the 
open circle line. This is quite interesting since it could open the possibility of 
observing some experimental signatures of the neutron wave effect. 

For a given asymmetry I the value of is different for the models con- 

sidered, as can be argued by the behaviour of 5pp/5pn at 2po in Fig-7-l(b). 
E.g. for J = 0.1 NLpS has the lower value (p^°^* — 2.4po), while NLp has 
the higher one (p^°** — 3.0po)- This is again related to the reduction of the 
isovector restoring force when the scalar-isovector channel (5-like) is present, 
as discussed before. 
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Fig. 7-2. Sound phase velocities of the propagating cohective mode vs. the baryon 
density (NLp case). Crosses: isovector-hke. Open circles: isoscalar-like. (a): sym- 
metric matter, (b): asymmetric matter, I = 0.1. (c): asymmetric matter, / = 0.5. 
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Fig. 7-3. Ratio of protons and neutron amplitudes in the propagating mode, for 
different asymmetries, as a function of the baryon density around the p^°** [NLp 
case). Crosses: / = 0.1, Fig.7-2b. Open circles: / = 0.5, Fig.7-2c. 



131 



l.Jf. Isoscalar Collective Modes in Asymmetric Nuclear Matter 
7.4-1 Exotic high baryon density modes 

From the previous analysis we have seen the isoscalar-hke excitations to be- 
come dominant at high baryon density, above the p^"*'* introduced before. 

Some results are shown in Fig. 7-4. It should be noticed that the frequency 
of the isoscalar-like modes is essentially related to the compressibility of the 
system at the considered density. In Fig. 7-4 (a) wc display the sound velocity 
obtained in Hartree and Hartree-Fock calculations at = 3.5 po; as a function 
of the asymmetry /. The differences observed among calculations performed 
within the Hartree or Hartree-Fock scheme for symmetric matter are due to 
a different behaviour of the associated equation of state at high density. At 
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Fig. 7-4. Sound velocity (a) and chemical composition (b) of isoscalar-like modes 
at pB = 3.5/90. Solid line: NLHF. Long dashed line: NLp. Dotted line: NLp5. The 
full circles in panel (b) represent the behaviour of Pp/ Pn vs. /. 

7 = the two Hartree models have exactly the same isoscalar mean fields, but 

for asymmetric nuclear matter the different behaviour of the symmetry energy 
leads to a different compressibility. The case NLp6, which has the stiffer Egy^ 
(resulting in a greater incompressibility for / > 0) with respect to NLp, shows 
also a larger increase of Vg/vpn with density. 
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Differences are also observed in the chemical composition of the mode (Fig.7- 
4(b)). The black circles show the behaviour of PjJ Pn vs. a. Note the almost pure 
neutron wave structure of the propagating sound, since the oscillations of pro- 
tons appear strongly damped {Spp/Spn -C Pp/Pn), see the previous comments. 

Before closing this discussion we have to remark that the isoscalar-like modes 
at high baryon density are vanishing if the nuclear EOS becomes softer. This 
is indeed the results of two recent models, [266,267], where the nuclear com- 
pressibility is decreasing at high baryon density because of a reduction of the 
isoscalar vector channel contribution. In ref. [266] this is due to self-interacting 
high order terms for the u meson, while in ref. [267] to a decreasing density- 
dependent coupling. 

Finally we note that all causality violation problems (superluminal sound ve- 
locities) observed in the non relativistic results at high baryon density, see 
ref. [263] and Fig. 3c in ref. [50], are completely absent in the relativistic ap- 
proach, as seen in the high density trends in Fig.7-2. 



7.4-2 Isospin Distillation in Dilute Matter 

As a good check of the relativistic approach we have also investigated the 
response of the system in the region of spinodal instability associated with 
the liquid-gas phase transition, which occurs at low densities. In this region 
an isoscalar unstable mode can be found, with imaginary sound velocity, that 
gives rise to an exponential growth of the fluctuations. The latter can rep- 
resent a dynamical mechanism for the multi-fragmentation process observed 
in heavy-ion collisions, see the discussion in Sects. 3 and 5. We find this kind 
of solution. In Fig. 7-5 we show the imaginary sound velocity and the ratio 
Spp/Spn as a function of the initial asymmetry for this collective mode. For 
all the interactions this ratio is different from the corresponding pp/ p„ of the 
initial asymmetry. This is exactly the chemical effect associated with the new 
instabilities in dilute asymmetric matter [19,20]. In particular it is seen that, 
when isoscalar-like modes become unstable, the ratio Spp/6pn becomes larger 
than the ratio Pp/ Pn (at variance with the stable modes at high densities, see 
Fig. 7-4). Hence proton oscillations are relatively larger than neutron oscilla- 
tions leading to a more symmetric liquid phase and to a more neutron rich gas 
phase, during the disassembly of the system. This is the Isospin Distillation 
effect in fragmentation, as already studied in detail in the previous Sects. 3 
and 5 in a non-relativistic frame. 

In conclusion we note that in dilute asymmetric NM we can distinguish two 
regions of instability, mechanical (cluster formation) and chemical (component 
separation). There is however no discontinuity in the structure of the unstable 
modes which are developing. For all realistic effective nuclear interactions (rel- 
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Fig. 7-5. Isoscalar-like unstable modes at pB = 0.4 po '■ Imaginary sound velocity 
(a) , in c-units, and ratio of proton and neutron amplitudes (b) as a function of the 
asymmetry. Solid line: NLHF. Dotted line: NLpS. Long dashed line: NLp. The full 
circles in panel (b) represent the behaviour of Pp/pn vs. the asymmetry parameter. 

ativistic and not) the nature of the unstable normal modes at low densities 
is always isoscalar-like, i.e. with neutrons and protons oscillating in phase, 
although with the neutron distillation effect discussed before, refs. [20,268]. 

7.5 General comments from the ANM collective response 

We have shown that from a detailed study of the collective response of asym- 
metric nuclear systems it would be possible to obtain information on the 
Lorentz structure of the in-medium interaction in the isovector channel. In fact 
we have singled out some qualitative new effects of the 5-meson-like channel 
on the dynamical response of ANM. Essentially, our investigation indicates 
that even if the symmetry energy is fixed, the dynamical response is affected 
by its internal Lorentz structure, i.e. the presence or not of an isovector- 
scalar field. This is implemented by the explicit introduction of an effective 
5-mcson and/or by the Fock term contributions. Both mechanisms are absent 
in the present rclativistic RPA calculations for finite nuclei. In the spirit of 
the EFT / DFT approach [229] it would be interesting to see the effect of an 
isovector scalar field extension, at the lowest order, on the existing covariant 
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RPA results. In general we see a close analogy in the structure of the linear 
response equations for isoscalar/isovector modes: 

• Same form of the dispersion relations. 

• Parallel role of E^"^ and Kf^lj in the determination of the restoring force. 

• Parallel structure of the corrections due to the scalar-vector meson compe- 
tition. 

This appears to be a beautiful "mirror" structure of the relativistic approach 
that seems to nicely support the introduction of a (5-meson-like coupling in the 
isovector channel, at least from a formal point of view. We like to remind that 
the same "mirror" structure of the relativistic picture is remarked in Sect. 6 
for equilibrium properties, saturation binding and symmetry energy, the ai 
and a4 parameters of the Bethe-Weiszacker mass formula. 
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8 Relativistic Heavy Ion Collisions: the coveiriant structure of the 
symmetry term 

Intermediate energy Heavy Ion Collisions, HIC, open the unique possibility 
to explore the Equation of State (EOS) of nuclear matter far from satura- 
tion, in particular the density dependence of the symmetry energy Esym{pB) 
. Within a relativistic transport model it is possible to see that the isovector- 
scalar 5-meson, which affects the high density behavior of the symmetry term, 
influences the isospin dynamics. The effect is greatly enhanced by a relativis- 
tic mechanism related to the covariant nature of the fields contributing to the 
isovector channel. 

An increasing -Esym(ps) leads to a more proton-rich neutron star whereas a de- 
creasing one would make it more pure in neutron content. As a consequence the 
chemical composition and cooling mechanism of protoncutron stars [273,274], 
mass-radius correlations [275,276], critical densities for kaon condensation in 
dense stellar matter [277,278] as well as the possibility of a mixed quark-hadron 
phase [279] in neutrons stars will all be rather different. It has recently been 
argued by means of simple thermodynamics considerations that even the onset 
of a quark-deconfined phase at high baryon density could present a sensitivity 
to the behaviour of Esym{pB) even for not very large asymmetries [241]. 

In the previous Sections we have seen how the search for EsymiPs) around 
saturation density has driven a lot of theoretical and experimental efforts. We 
have seen how HICs at Fermi energies can give the possibihty to extract some 
information on the symmetry term of the nuclear Equation of State (EOS) in 
region below and/or shgthly above the normal density. Here we will focus our 
attention on Relativistic HICs trying to select reaction observables particu- 
larly sensitive to the symmetry energy at higher density, where furthermore 
we cannot have complementary investigations from niiclcar structure like in 
the case of the low density behaviour. We stress again that HICs provide the 
unique way to create asymmetric matter at high density in terrestrial labo- 
ratories. Moreover effective interactions for high momentum nucleons can be 
probed. Calculations within transport approaches show that HICs around 
lAGeV allow to reach a transient state of matter with more than twice the 
normal baryon density. Moreover, although the data are mostly of inclusive 
type (and the colliding nuclei not very neutron rich), quite clearly a depen- 
dence of some observables on charge asymmetry is emerging. 

Collective flow, particle production and isospin equilibration results for reac- 
tions induced by stable and radioactive beams are discussed. The elliptic flows 
of nucleons and light isobars appear to be quite sensitive to the microscopic 
structure of the symmetry term, in particular for particles with large trans- 
verse momenta, since they represent an earlier emission from a compressed 
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source. Thus future, more exclusive, experiments with relativistic radioactive 
beams should be able to set stringent constraints on the density dependence 
of the symmetry energy far from ground state nuclear matter. 



Moreover we show that in fact a relativistic description of the nuclear mean 
field can account for an enhancement of isospin effects during the dynam- 
ics of heavy-ion collisions. The isospin dependence of collective flows and 
pion production has been already discussed in a non-relativistic framework 
[46,113,280,281] using very different EOS with opposite behaviours of the 
symmetry term at high densities, increasing repulsion (asy-stiff) vs. increas- 
ing attraction (asy-soft). The main new result shown here, in a fully relativis- 
tic scheme, is the importance at higher energies of the microscopic covariant 
structure of the effective interaction in the isovector channel: effective forces 
with very similar symmetry terms can give rise to very different isospin effects 
in relativistic heavy ion collisions. 

In recent years some efforts have been devoted to the effects of the scalar- 
isovector channel in finite nuclei, [220,230,231]. Such investigations have not 
shown a clear evidence for the 5-field and this can be understood considering 
that in finite nuclei one can test the interaction properties mainly below the 
normal density, where the effect of the 5— channel on symmetry energy and on 
the effective masses is indeed small [12,13] and eventually could be absorbed 
into non linear terms of the p field. In Sect. 7 we have also seen that the 
dynamical collective response is not much affected, in particular the spinodal 
instabilities. We can expect just some weak effects on the isovector-hke dipole 
modes. 

Here we show that heavy-ion collisions around lAGeV with radioactive beams 
can provide instead a unique opportunity to spot the presence of the scalar 
isovector channel. In fact, due to the large counterstreaming nuclear currents 
one may even exploit the different Lorentz nature of a scalar and a vector 



When both p-like and 5-\ike channels are considered EsymipB) can be written 
as (Sect. 6): 



with E* = sJk'^F + M*'^, M* the effective Dirac mass and /p,^ = {gp^s/nip^s)'^ 

are the coupling constants of the isovector channels. We see that, when the 
5-meson is included, the observed 04 = Esym{po) value actually assigns the 
combination [fp — fsi^)"^] of the (p, 5) coupling constants, [12,13]. In Fig.8-1 
we show the density dependence of the symmetry energy for three different 
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Fig. 8-1. Total (kinetic+potential) symmetry energy (in MeV) as a function of the 
baryon density. Solid: NLpS. Dashed: NLp. Short Dashed: NL — Dp. In the insert 
the density behaviour of the p coupling, fp (in frn^), for the three models is shown. 

Non — Linear — RMF models: two of these, one including only the p field 
(NLp), the other with p and S fields (NLpS), were already introduced in 
Sect. 6 (see Tabic 7). The last, {NLDp), includes only a p field but with a 
covariant density dependence of fp, [220,258]. This is tuned to give at high 
density the same Esym{pB) of the the {NLp6) case. As shown in the following, 
the latter model is useful for disentagling in the reaction dynamics the effects 
due to a difference in EgymiPB) from those directly linked to the strength of 
the p vector field [135]. Thus these models parametrize the isovector mean field 
either by only the vector field with fp = 1.1/m^, or with a balance between 
a vector field with fp = 3.3/m^ and a scalar one with fs = 2.4/m^, or finally 
by a normal density coupling fp{po) = 1.1/m^, at saturation, but with an 
increasing density dependence as shown in Fig. 8-1 (insert). We stress again 
that in NLpS the symmetry energy results from a balance between a scalar 
attraction, {6 — like), and a vector repulsion, (p — like), where the p-coupling 
is now roughly three times larger than in the NLp case. 

The strength of fs is fixed relatively well by DBHF [34,35] and DHF [36] 
calculations. Therefore the effects described in the following are not artificially 
enhanced, but based on a reliable estimate available at the moment. In any 
case the aim of our work is to present some qualitative new features expected 
in the reaction dynamics, in particular for collective flows, particle production 
and isospin stopping power, which result from the introduction of a scalar 
effective field in the isovector channel. 

We finally note that the use of a RMF approach with fixed couplings in high 
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density regions generally requires some caution due to the density dependence 
of the virtual meson couphngs because of correlations beyond the mean field 
scheme. In particular a softer EOS for the isoscalar part is favoured, as sug- 
gested from experimental analyses [129], from DBHF predictions and from 
reaction simulations for collisions of charge symmetric ions at intermediate 
energies [131]. In the calculations presented here the (a, a;) coupling constants 
have been modified accordingly, see details in ref. [232]. 



8.1 Relativistic transport simulations 



For the theoretical description of heavy ion collisions we solve the covari- 
ant transport equation of the Boltzmann type [283,146] within the Relativis- 
tic Landau Vlasov (RLV) method [259] (for the Vlasov part) and apply- 
ing a Monte-Carlo procedure for the collision term. RLV is a test particle 
method using covariant Gaussians in phase space for the test particles. The 
collision term includes elastic and inelastic processes involving the produc- 
tion/absorption of the A(1232Mey) and A'"*(1440Mey) resonances as well as 
their decays into one- and two-pion channels. Details about the cross sections 
for all the possible channels can be found in ref. [284]. An explicit isospin- 
dependent Pauli blocking term for the fermions is employed. 

For the following discussion it is useful to recall how a relativistic Vlasov 
equation can be obtained from Wigner Function dynamics of Sect. 6. The neu- 
tron/proton Wigner functions are expanded in terms of components with def- 
inite transformation properties. Consistently with the effective fields included 
in the minimal model one can limit the expansion to scalar and vector parts: 

From the kinetic equation (6-5) of Sect. 6 we obtain after some algebra a rela- 
tion between the vector and scalar components 

~ M* 

and thus a transport equation of Vlasov type for the scalar part fi{x,p*^) = 

{p^d'^ +\pUFr + M:{d^M:)]d/}Mx,p*^) = o (8-2) 



with the field tensors 

Fr = 9>r + d^pT- 

The trajectories of test particles obey to the following equation of motion: 
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M*{ 

. Pi 
' M* 



^x) 
(x) 



Fr(x,{r))+d'^M:(x) . 



(8-3) 



In order to have an idea of the dynamical effects of the covariant structure of 
the interactive fields, we write down, with some approximations, the "force" 
acting on a particle. Since we are interested in isospin contributions we will 
take into account only the isovector part of the interaction [135]: 



dr 

E* - 



Pi^ 
M* 



(8-4) 



In the second line we have neglected the contribution coming from the gradient 
of the current in the transverse direction, i.e. terms like (V^Js^) {I ^ m — 
X, y, z) and the derivative of the current with respect to time. This form will 
be used below to interpret collective flow. 



8.2 Collective flows 



Collective flows in heavy ion collisions give important information on the dy- 
namic response of dense, hot and asymmetric nuclear matter [121,123,124,125]. 
In particular the proton-neutron differential flow F^"(?/) [127] has been found 
to be a very useful probe of the isovector part of the EOS since it appears 
rather insensitive to the isoscalar potential and to the in medium nuclear cross 
section. The deflnition of the FP^{y) is 

1 Niy) 



where N{y) is the total number of free nucleons at the rapidity y, p^i is the 
transverse momentum of particle i in the reaction plane, and Tj is and -1 
for protons and neutrons, respectively. 

A typical result for the very asymmetric ^^"^871 -|-^^^ Sn reaction at l.bAGeV 
(semicentral collisions) is shown in Fig. 8-2. The error bars represent the statis- 
tical fluctuations due to the Monte-Carlo nature of the simulations. We notice 
that the differential flow in case of the NLp6 presents a stiffer behaviour rela- 
tive to the NLp model, as expected from the more repulsive symmetry energy 
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Esym{pB) at high baryon densities, in Fig.8-1. On the other hand it may be 
surprising that a relatively small difference at 2po results in a such different 
collective flows. 



The calculation can be repeated using the NLDp interaction, i.e. with only 
a p contribution but tuned to reproduce the same EOS of the NLpS case. 
The results, short-dashed curve of Fig. 8-2, are very similar to the ones of the 
NLp interaction. Therefore wc understand the large flow effect as mainly due 
to the different strengths of the vcctor-isovcctor field in the NLpS and the 
NLp, NLDp models. In fact if a source is moving the vector field is enhanced 
(essentially by the local 7 = -^^ Lorentz factor) relative to the scalar one. 

In the second line of Eq.(8-4) we have written the isovcctor contribution to the 
effective force. We are actually interested in the difference between the forces 
acting on a neutron and on a proton. Oversimplifying the HIC dynamics 
we consider locally neutrons and protons with the same 7 factor (i.e. with the 
same velocity). Then Eq.(8-4) can be expressed approximately in the following 
transparent form (^53 ~ ^pa): 



dr dr 



7 



Vp3 = 



-El± + 2(7 - l)/p 



PB 



Vp3 (8-6) 



where 7 = ^ is the local Lorentz factor of the collective motion and Eq.(8-1) 
has been used. 

Keeping in mind that NLp5 has a three times larger p-coupling it is clear 
that dynamically the vector-isovector mean field acting during the HIC is 
much larger than the one of the NLp, NLDp cases. Then the isospin effect 
is mostly caused by the different Lorentz structure of the "interaction" which 
results in a dynamical breaking of the balance between the p vector and S 
scalar fields, present in nuclear matter at equilibrium, [135]. This effect is anal- 
ogous to the interplay between the isoscalar vector- and scalar-fields which 
is seen in the magnitude and energy dependence of the real part of the op- 
tical potential, [131]. To characterize the effect on differential collective flows 
we have calculated the slope dF^"'{y) / d{y / Uproj) at mid-rapidity. Its value is 
46.7Mey/c for NLp6 and 23AMeV/c for NLp, i.e. a factor two difference. 
Calculations performed at lower beam energies show that below 500AMeV 
there is no essential difference in the differential flow predictions among the 
models discussed here. The effect from the strength of p fleld starts to become 
important around lAGeV, as expected from a relativistic mechanism. 

Another interesting observable is the elliptic flow V2{y,Pt), the second coeffi- 
cient in a Fourier series expansion of the azimuthal emission distribution, see 
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y/ yproj 

Fig. 8-2. ^^'^Sn+^^'^Sn reaction at 1.5 AGeV {b = 6fm): Proton-neutron differential 
transverse collective flow (in MeV/c), vs. rapidity, for the three different models for 
the isovector mean fields. Full circles and solid line: NLp6. Open circles and dashed 
line: NLp. Stars and short dashed line : NLDp. 
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Fig. 8-3. Difference between neutron and proton elliptic flow as a function of the 
transverse momentum in the ^^'^Sn -j-^^^ reaction at 1.5 AGeV (b=6fm) in the 
rapidity range —0.3 < y/ yproj ^ 0.3. Full circles and solid line: NLp5. Open circles 
and dashed line: NLp. Stars and short dashed line: NLDp. The statistical error 
bars of the NLDp curve are similar to the other cases and not shown. 

Sect. 4, 



Pt 

where pt = ^pl + Py is the transverse momentum [128,129]. A negative value 
of V2 corresponds to the predominant emission of matter perpendicular to 
the reaction plane (the so-called squeeze — out flow). The p^-dependence of 

V2, which has been recently investigated by various groups [129,124,130,131] 
is very sensitive to the high density behavior of the EOS since highly ener- 
getic particles {pt > 0.5) originate from the initial compressed phase of the 
collision, [131]. 
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In Fig. 8-3 we present thep^ dependence of the proton-ncTitron difference of tfie 
elliptic flow in the same very exotic ^^'^Sn+^^'^ Sn reaction at 1.5AGeV (semi- 
central collisions) for mid-rapidity emissions. The larger error bars correspond 
to a reduced statistics when a selection on different pt bins is introduced, in- 
creasing for larger pt values due to the smaller number of contributions. Prom 
Fig. 8-3 we see that with the (p -|- S) dynamics the high-pt neutrons show a 
much larger squeeze — out. This is fully consistent with an early emission 
(more spectator shadowing) due to the larger repulsive p-field. We understand 
this large effect since the relativistic enhancement discussed above is relevant 
expecially at the first stage of the coUision. The V2 observable, which is a good 
chronometer of the reaction dynamics, appears to be particularly sensitive to 
the Lorentz structure of the effective interaction. 

The same set of simulations has been repeated in ref.[135] for the more realistic 
^^^Sn +^2^ Sn reaction at 1.5 AGeV {b = 6/m), that hkely could be studied 
with the new planned radioactive beam facilities at intermediate energies. 
The results are shown in Fig. 8-4. The effect of the different structure of the 
isovector channel is still quite clear, of course with a reduction due to the 
smaller isospin density in the interaction region. Particularly evident is again 
the splitting in the high pt region of the elliptic flow. 



8.3 7r"/7r+ Ratios 



Observable effects originating from the high density symmetry energy are re- 
lated to differences in neutron and proton densities. Thus, we consider in the 
following the ratio of neutrons to protons as a function of time and space. One 
expects isospin effects on the different isospin channels of pions, since they are 
produced in nn, pp and np collisions via the decay of A and N* resonances. 
Using the same relativistic transport code the 7r~ vs. 7r+ production for central 
Au + Au collisions at different energies can be evaluated, sec [232]. Results are 
shown in Fig. 8-5. For the NLpd model the neutrons are emitted much earlier 
than protons from the high density phase due to a more repulsive mean held 
together with a lower n-effective mass. This effect is responsible for the reduc- 
tion of the n/p-ratio in the residual system and, particularly, it influences the 
particle production, [280,281]. 

The four isospin states of the A-resonance are produced in different scattering 
channels, when the threshold energy is available, e.g. nn — > pA~,nA°, pp 
pA"^, nA''"+, • • •. Thus A°'~ (A+'+"^) resonances are mainly formed in energetic 
nn- {pp-) collisions. Therefore, the n/p-ratio is indirectly related to that of the 
particle production. 

The n/p effective mass splitting in asymmetric matter can also directly affect 
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Fig. 8-4. ^^'^Sn +^^^ Sn reaction at 1.5 AGeV {b = 6/m) from the three different 
models for the isovector mean fields. Top: as in Fig. 8-2. Bottom: as in Fig. 8-3. Full 
circles and solid line: NLp5. Open circles and dashed line: NLp. Stars and short 
dashed line: NL — Dp. Error bars: see the text and the previous caption. 

the resonance production because of threshold effects. We note that resonances 
will also have isospin dependent in-medium effective masses that can be related 
to the nucleon effective masses through the isospin coupling coefficients in the 
process A <-> nN, [134,232]. In terms of the self-energies E we have: 



S,(A-) = S,(n) 

E,(AO) = h,{n) + 1e,(p) 

E,(A+) = lE,(n) + ^S,(p) 
Si(A++) = Si(p) and 

E,(7V*(+'°)) = E,(p,n). (8-7) 
where i — scalar and vector. 

Thus the effective masses of the four isospin states of the A-resonance will be 
different when a 5-meson is included in the calculation. We have seen from 
general considerations (see Sect. 6), that, in n-rich systems, a isovector scalar 
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Fig. 8-5. Central Au + Au collisions at O.GAGeV (upper) and l.OAGeV (bottom). 
Time evolution of n/p ratio and A resonance production in high density regions 
(p/po ^ 2.0) (first two columns) and of the total 7r~/7r+ ratio (right). Solid lines: 
NLp. Dashed lines: NLpS. 

meson field is leading to a neutron effective mass smaller than that for protons, 
in particular at high baryon densities. Consequently in the process nn — > pA~ 
less energy will be available for the A~ production in the NLpS case. We hke 
to note that this mechanism is not strictly linked to the density behavior of 
the symmetry energy. E.g., in density dependent coupling models we can have 
at high densities a lower Esym{p) because the p-meson coupling is decreasing 
but still a large m* — m* splitting if the 5-meson coupling stays constant or 
it is slightly increasing, as e.g. in the DBHF calculations of ref.[232]. 



In the A-resonance multiplicities (middle figures) one thus sees a decrease of 
the A°'^ isospin states due to the effect of the isovector-scalar (5-mcson. The 
pions are mainly produced from resonance decays, e.g. as A~ — > nvr", A° — > 
p7r~, A+ n7r+, • • •. Therefore the decrease of A'^'" resonances reduces 
the production of negative encharged pions which consequently decreases the 
7r~/7r^-ratio (right part in Fig. 8-5). This ratio then appears to be sensitive to 
the isospin term of equation of state at high densities. One should note that 
both models used here exhibit a similar asy-stiff behavior at high densities (see 
Fig. 8-1). We deduce that the n/p effective mass splitting mechanism described 
above is very important in inducing isospin effects on the particle production 
[285]. 

A comparison of total tt^/ tt" ratio with preliminary data, at different energies, 
is presented in Fig.8-6, in central Au -\- Au collisions. We note that the 5 
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Fig. 8-6. Energy dependence of the (7r"/7r^)-ratio in central Au + Au collisions. Re- 
sults of NLp (solid) and NLpd (dashed) parametrizations. The grey bars correspond 
to some very preliminary data from the FOPI — GSI collaboration [286]. 

influence (difference between solid and dashed curves) is slightly decreasing 
with the beam energy. This is an indication of the n/p effective mass splitting 
mechanism in A~ production, expected to be more important at the threshold. 

Generally isospin effects on pion production are decreasing with incident en- 
ergy, as seen in Fig.8-6, due to secondary collisions, i.e. pion absorption and 
A-rescattering, which are essentially not isospin-dependent. Thus the isospin 
dependence of pion production will be on average moderated. This trend was 
already found in earlier studies using a non-relativistic Quantum-Molecular 
Dynamics (QMD) approach [280]. 

For the same systems new very preliminary data on pion production have 
recently been reported from the FOPJ-coUaboration at GSI [286], shown as 
grey bars in the Fig.8-6. Generally our simulations provide a qualitative good 
description of the pion ratio with respect to beam energy, except perhaps 
at the lower energies. The inclusion of a 5 meson in the iso-vector part of 
the equation of state improves the comparison, at least at higher energies. In 
any case it appears that 7r~/7r """-ratios might be a probe of the high density 
symmetry energy and the role of the virtual 5-meson. It would be helpful to 
look at more exclusive data on pion production in asymmetric systems, [232]. 

The effect of the (5-field as discussed above is enhanced in the momentum 
spectra, i.e. the rapidity and transverse momentum distributions, due to the 
interplay of the isospin effects with the Coulomb interaction which affects 
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differently the tt"*" and tt^ distributions. Thus a clearer and more observable 
effect of the 5-field is predicted in the flow observables of the (7r~/7r"'")-ratios, 
as a function of energy. 



8.4 Isospin Transparency 

Isospin equilibration is seen in the collisions of systems with different asym- 
metries. It has been extensively investigated by the FOPI collaboration as a 
signature of trasparency [199,289]. Here we investigate it also with respect to 
the sensitivity to the iso — EOS. The idea is to study colliding systems with 
the same mass number but different N/Z ratio In particular, a combination 
of 44-Rii, N/Z — 1.18 and %Zr, N/Z — 1.4 has been used as projectile/target 
in experiments at intermediate energies of 0.4 and 1.528 AGeV [199,289]. The 
degree of stopping or transparency has been determined from the rapidity de- 
pendence of the imbalance ratio for the mixed reactions Ru{Zr) + Zr{Ru): 
i?(l/(o)) = ^RuZT^y{o)y^ZrRu^y{o)^^^ ^^^^^ N\y^^^) is the particle yield at a 

given rapidity for Ru -\- Zr and Zr -\- Ru with i — RuZr, ZrRu. The observ- 
able R is measured for different particle species, like protons, neutrons, light 
fragments such as t and ^He and produced particles such as pions {t:^'^), etc. 
It characterizes different stopping scenarios. 

In the proton case, moving from target to cm rapidity, R{p) rises (positive 
slope) for partial transparency, falls (negative slope) for full rebound scenarios 
and it is fiat when full stopping and total isospin mixing is achieved in the 
collision. An opposite behavior will appear for neutrons. Indeed we recall that 
in the full transparency limit R should have the value of R{p) = Z'^^'/Z^" = 
40/44 = 0.91 and R{n) = N^^/N^"" = 56/52 = 1.077 for protons and neutrons 
at target rapidity, respectively. Therefore, R{p) can be regarded as a sensitive 
observable with respect to isospin diffusion, i.e. to properties of the symmetry 
term. 

For the investigation of isospin equilibration we have analyzed transport re- 
sults obtained in the same way as carried out in the FOPI experiment, 
refs. [287,288]. In particular, central events are selected through the observ- 
able ERAT, the ratio of the mean transverse to the mean longitudinal kinetic 
energy. For this a phase space coalescence algorithm for fragment production 
in the final state is used, and the observable ERAT is calculated for each frag- 
ment. A detailed description of this analysis can be found in ref.[131], where 
it was shown that the charged particle multiplicity and ERAT distributions 
fit well the experimental data, as an important check of the phenomenolog- 
ical phase space coalescence model for fragment production. Thus the same 
centrahty cuts as in the FOPI experiments [289] can be used. 
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as function of the normalized 



Fig. 8-7. The imbalance ratio R{y^^'>) - j^zrRu^yW^ 

rapidity y^'^^ = y/yproj of free protons (left panels) and free neutrons (right panels) 
for the NLp and NLp6 models, for central (6 < 2 fm) Ru{Zr) + Zr(iiii)-collisions 
at 0.4 (top) and 1.528 AGeV (bottom) beam energies. The experimental data are 
taken from the FOP I collaboration [289]. The ratio in the initial system is shown 
by the horizontal dotted line. 




Fig. 8-8. Proton imbalance ratio as function of the normalized rapidity for central 
ib < 2 fm) Ru{Zr) + Zr(i2ti)-collisions at 0.4 AGeV as in the top-left part of 
Fig. 8-7. A curve is added (empty squares) corresponding to a NLp calculation with 

2 '-''/ree) • 



reduced nucleon-nucleon cross sections (half the free values, aNN{E) — ^■ 
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In Fig. 8-7 wc show the rapidity dependence of the imbalance ratio for free pro- 
tons {R{p)) and free neutrons {R{n)) at the two energies 0.4 and 1.528 AGeV. 
The imbalance ratio approaches unity at mid-rapidity for all particle types due 
to symmetry in the mid-rapidity region, as expected in central collisions. 

Going from target- to mid-rapidity the ratio nicely rises for protons, and de- 
creases for neutrons, a good signature of isospin transparency. The effect is 
more evident for the NLpS interaction. The observed difference between the 
two models can be understood since within the NLp6 picture neutrons expe- 
rience a more repulsive iso- vector mean field, particularly at high densities, 
than protons, and consequently much less nucleon stopping in the colliding 
system. 

However the influence of the iso-EOS on the imbalance ratio of protons and 
neutrons is not very large. At low intermediate energies (0.4 AGeV) one deals 
with moderate compressions of < 2 • psat where the differences in the iso- 
vector EOS arising from the S meson are small. At higher incident energies 
(1.528 AGeV, Fig.8-7-bottom) where a larger effect is expected, we actually 
see a slightly higher isospin effect on the imbalance ratios, at least for protons. 
In fact with increasing beam energy the opening of inelastic channels via 
the production/decay of A resonances through pions, as discussed in the last 
section, also contribute to the final result. This interpretation is confirmed by 
other studies [281]. 

Reduced in-medium Nucleon-Nucleon (NN) cross sections, in particular (T„p, 
will also increase the isospin transparency. This possibility can be investi- 
gated considering a factor of two reduction, a = ^(Tfree, of the free A^A^-cross 
section values used before, [288]. We note that such a reduction represents 
rather an upper limit of in-medium effects as compared to recent microscopic 
Dirac-Brueckner estimations [290]. In Fig.8-8 we show the results for the pro- 
ton imbalance ratios at 0.4 AGeV with 0.5(7/ree in the NLp case. We see an 
overall slightly increased transparency but not enough to reproduce the trend 
of the experimental values in the target rapidity region. On the other hand 
the reduction of the A^A^ cross sections, and in particular of cr(„p), leads too 
large a transparency in the proton rapidity distributions for central collisions 
of the charge symmetric Ru + Ru case; as is seen in our simulations and 
already remarked in previous IQMD calculations, [289]. Thus can exclude 
further reductions of the A^A^ cross sections. 

Since the calculations are performed with the same EOS for the symmet- 
ric nuclear matter, the same compressibility and momentum dependence, the 

observed transparency appears to be uniquely related to isovector — EOS ef- 
fects, i.e. to the isospin dependence of the nucleon self-energies at high baryon 
densities. The fact that protons and neutrons exhibit an opposite behavior 
for the imbalance ratios at target rapidity suggests that the detection of the 
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Fig. 8-9. Imbalance ratio of t to ^He for the same collision as in Fig. 8-7 for 0.4 AGeV 
beam energy. 
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Fig. 8-10. Imbalace ratio for the ratio of negative to positive charged pions for the 
same collision as in Fig. 8-7 for 0.4 AGeV beam energy. 

imbalance observable R{t/^He), i.e. the double ratio of the t/^He yield should 
reveal a larger sensitivity. The correct and practical method how to properly 
describe light fragment formation is still controversial. We report here on re- 
sults obtained in ref. [288] with simplest algorithm, namely a phase space 
coalescence model [131]. 

Fig. 8-9 shows the imbalance ratio R{t/^He) for central collisions at 0.4 AGeV 
incident energy. The isospin effect from the inclusion of the iso-vector, scalar 
6 meson in the NLp6 model is found to be very large near target rapidities. 
We note that the effect indeed can be hardly seen from the separate imbalance 
ratios for protons and neutrons at the same rapidity (see Fig. 8-7 top panels) 
apart from the difficulties of neutron detections. It would be therefore of great 
interest to experimentally measure directly this quantity. 




Finally, another sensitive observable should be, from the discussion in the 
last Section, the imbalance ratio of charged pions R{'k~ /ti^) (Fig. 8-10). At 
variance with the previous results for neutrons and light isobars, this ratio is 
reduced at target rapidity with the NLp6 model. This effect is consistent with 
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our understanding of the /t:~- ratios. Pions are produced from the decay 
of A resonances formed during the high density phase, see Fig.8-5. The tt^ 
abundance is then hnked to the neutron-excess of the high density matter, 
as discussed in the previous Section. We recall that the contribution of the 5 
meson leads to a more repulsive field for neutrons at supra-normal densities 
and consequently to less neutron collisions and finally to a smaller t:~ /t:^ 
ratio. 

We have analyzed the isospin transparency in relativistic collisions. We have 
observed that this observable is sensitive to the microscopic Lorentz structure 
of the symmetry term. Effective interactions with symmetry energies which are 
not much different at 2-3 times normal density Psat predict large differences 
in the isospin — transparencies, depending on the relative contribution of 
the various charged vector and scalar fields. Intermediate energy heavy-ion 
collisions with radioactive beams can give information on the symmetry energy 
at high baryon density and on its detailed microscopic structure. We have 
shown that such experiments provide a unique tool to investigate the strength 
of the 5 — like field. The sensitivity is enhanced relative to the static property 
Esym^ps) because of the covariant nature of the fields involved in collision 
dynamics. 



151 



9 Conclusion and outlook 



Nuclear reactions with neutron-rich (or radioactive) nuclei have opened the 
possibility to learn about the behaviour of the nuclear interaction and, in par- 
ticular, of the symmetry energy in a wide spectrum of conditions of density 
and temperature. This study appears extremely important: indeed a deeper 
understanding of the behaviour of neutrons and protons in a charge asymmet- 
ric nuclear medium is essential to test and to extend our present knowledge 
of the nuclear interaction and is of highest importance for the modelling of 
astrophysical processes, like supernova explosion or neutron stars. 

In the context of nuclear reactions, one has to identify the most sensitive 
observables that may provide information on the behaviour of the symmetry 
energy in several conditions of density and temperature. This is the main line 
followed along this Report. 

In Section 2 we discuss the symmetry energy dependence around normal den- 
sity and how this can affect important properties of neutron-rich nuclei, such 
as compressibility (and monopole frequency), saturation density and neutron 
skin. 

The behaviour of asymmetric matter at low density has been investigated in 
Section 3. We discuss in particular the phase diagram of asymmetric matter 
and the relevant features of instabilities. This subject appears important in 
connection to the possibility to observe a liquid-gas phase transition in vio- 
lent heavy ion reactions, where, after the initial coUisional shock, low density 
regions can be easily reached during the expansion phase. New important 
features, such as the isospin distillation effect, arc predicted and some exper- 
imental evidences have already appeared along this direction, though a more 
careful analysis is still needed, to disantangle among other possible contribu- 
tions to the distillation, such as pre-equilibrium effects, and to really prove 
the mechanism driving the fragmentation process. 

In Section 4 we have focused on features of the early stage of the reaction dy- 
namics between neutron-rich nuclei. Pre-equilibrium emission and collective 
flows appear particularly sensitive also to the the momentum-dependent part 
of the interaction. In asymmetric matter a splitting of neutron and proton 
effective masses is observed. The sign of the splitting is quite controversial, 
since the behaviour of neutron and proton optical potential at large energy 
has not been experimentally measured yet. Hence it appears very important 
to try to extract information on this fundamental question from nuclear colli- 
sions, where one can use probes, such as pre-equilibrium particles, particularly 
sensitive to the high density phase, where also high momenta are reached. 

In Section 5, we have explored several fragmentation mechanisms, occurring 
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at the Fermi energies, in the framework of a stochastic mean-field approach. 
We discuss the features of muhifragmentation in neutron-rich systems, and in 
particular the isotopic content of fragments. This can be connected to the be- 
haviour (the slope) of the symmetry energy at low density. For semi-peripheral 
reactions an interesting neutron enrichment of the overlap ( "neck" ) region ap- 
pears, due to the neutron migration from higher (spectator regions) to lower 
(neck) density regions. Also this effect is nicely connected to the slope of the 
symmetry energy. A careful comparison with experimental data would give 
important indications on the fragmentation mechanism and on the behaviour 
of the symmetry energy. 

In sections 6-7 we have discussed static properties and dynamical mechanisms 
in the context of relativistic approaches including isovector, both vector and 
scalar, channels (the p and 5 mesons) . The contribution of these two channels 
to static or dynamical properties has a different weight, leading to interesting 
effects on the frequency of the stable isovector modes, that is not simply 
related to the value of the symmetry energy (see Sections 6, 7). In the low 
density region, instabilities and the distillation mechanims are observed, in 
close parallehsm with the results of the non-relativistic treatment discussed 
above. 

Finally in Section 8 we discuss isospin effects in relativistic heavy ion collisions. 
The observable consequences of the inclusion of the 5 meson are enhanced by 
the Lorentz structure of the effective nuclear interaction in the isovector chan- 
nel. In particular, effects on particle production, such as pions or kaons, light 
particle collective flows and isospin diffusion are investigated. Once again, from 
the comparison with experimental data, one could get important information 
on the structure of the isovector interaction. 



9.1 Outlook: The Eleven Observables 

In conclusion, the study of isospin effects on static and dynamical nuclear 
properties appears as a very rich and stimulating field. Several probes can be 
used to get an insight on the behaviour of the symmetry energy in different 
conditions of density and temperature. A joint effort, from the experimental 
and theoretical side, should allow to extract relevant information on funda- 
mental properties of the nuclear interaction. We like to suggest a selection of 
Eleven Observables, from low to relativistic energies, that we expect partic- 
ularly sensitive to the microscopic structure of the in medium interaction in 
the isovector channel, i.e. to the symmetry energy and its "fine structure" : 

1. Competition of Reaction Mechanisms. Interplay of low-energy dissipative 
mechanisms, e.g. fusion (incomplete) vs. deep-inelastic vs. neck fragmentation: 
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a stiff symmetry term leads to a more repulsive dynamics. (Sect. 5) 

2. Energetic particle production. N/Z of fast nucleon emission: symmetry re- 
pulsion of the neutron/proton mean field in various density regions. Moreover 
at the Fermi energies we expect to see also effects from the n/p splitting of 
the effective masses. Even the spectra and yields of hard photons produced 
via (n,p) bremmstrahlung should be sensitive to the density and momentum 
dependence of the symmetry fields. (Sect. 4) 

3. Neutron/Proton correlation Junctions. Time-space structure of the fast par- 
ticle emission and its relation to the baryon density of the source. Again com- 
bined effects of density and momentum dependence of the symmetry term are 
expected (Sect. 4). 

4- E- slope of the Lane Potential. A systematic study of the energy dependence 
of the {n/p) optical potentials on asymmetric nuclei will shed lights on the 
effective mass splitting, at least around normal density. (Sects. 2 and 6) 

5. Isospin Distillation (Fractionation). Isospin content of the Intermediate 
Mass Fragments in central collisions. Test of the symmetry term in dilute mat- 
ter and connection to the possibility to observe a liquid-gas phase transition. (Sect. 5) 

6. Properties of Neck- Fragments. Mid-rapidity IMF produced in semicentral 
collisions: correlations between N/Z, alignement and size. Isospin effects on 
the reaction dynamics and "Isospin Migration" .(Sect. 5) 

7. Isospin Diffusion. Measure of charge equilibration in the "spectator" re- 
gion in semicentral collisions. Test of the interplay between concentration and 
density gradients in the isospin dynamics. (Sect. 5) 

8. Neutron- Proton Collective Flows. Together with light isobar flows. Check 
of symmetry transport effects. Test of the momentum dependence (relativistic 
structure) of the interaction in the isovector channel. Measurements also for 
different pt selections. (Sects. 4 and 8) 

9. Isospin Transparency. Measure of isospin properties in Projectile/Target 
rapidity regions in central collisions of "mirror" ions at intermediate energies. 
Similar effects, but due to the nucleon-nucleon cross section are expected to 
be smaller (Sect. 8). 

10. tt^/tt^ Yields. Since tt" are mostly produced in nn collisions we can expect 
a reduction for highly repulsive symmetry terms at high baryon density. Im- 
portance of a Pt selection. Similar studies for mesons with smaller rescattering 
effects would be of great interest. (Sect. 8) 

11. Deconfinement Precursors. Signals of a mixed phase formation (quark- 
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bubbles) in high baryon density regions reached with asymmetric HIC at in- 
termediate energies, versus the properties of the interaction considered (Sect. 8). 

We stress again the richness of the phenomenology and nice opportunities of 
getting several cross-checks from completely different experiments. 

For the points 3, 5, 6, 7, 8, 9, 10 from the transport simulations discussed 
here we presently get some indications of asy-stiff behaviors, i.e. increasing 
repulsive density dependence of the symmetry term, but not more fundamental 
details. Moreover all the available data are obtained with stable beams, i.e. 
within low asymmetries. 
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